Corrections with respect to the version submitted as Habihtationsschrift: 

• typographical errors corrected; 

• in section 2.5, the order of 2 sentences changed and the word 'rather' 
added; 

• in section 3.6, one citation for as updated; 

• in section 6.3, mistakenly introduced 'In general' substituted by 'For 
Wilson fermions'; 

• remark about the nucleon mass in quenched chiral perturbation theory 
moved from section 4.1.1 to section 5.5; 

• in section 4.2.2 modified statement about use of non-relativistic effec- 
tive field theories in lattice calculations of the T; 

• added a sentence introducing the decay constant data in the beginning 
of section 8.1.2; 

• new unquenched result for fs added in section 8.4. 

• typos in caption to Fig. 12 and in Fig. 30 corrected, orthographical typos 
corrected. 

• 1/Af fit to decay constants added. 
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1 INTRODUCTION 

One of the most appealing ideas is that the fundamental forces of Nature can be explained 
from mathematical symmetries. Indeed the known conservation laws result from space- 
time and internal symmetries, and the hope is that it is possible to find a consistent 
theory describing the interactions between all elementary particles and describing how 
the world is kept together. A large step forward on this way was the development of the 
Standard Model of elementary particle physics which includes Quantum Chromodynamics 
(QCD), the theory of strong interactions, based on a local SU{3) symmetry, and the 
Glashow-Salam- Weinberg model of the electromagnetic and weak forces based on a local 
f/(l) X SU{2) symmetry. It gives a surprisingly accurate description of the presently 
known elementary particle interactions, in experiments with particle collisions at energies 
up to several 100 GeV. One now tries to understand as much of the implications of the 
Standard Model as possible, and to test the model by enhancing the precision of theoretical 
predictions and experimental measurements and by searching for contradictions, to be able 
to find indications for a broader and perhaps simpler and more symmetric physical theory. 

An obvious expectation from the Standard Model, and QCD in particular, is to be 
able to explain the properties of mesons and baryons. For example, hadron masses and 
hadronic reactions should be predictable from QCD, and weak and electromagnetic matrix 
elements of hadrons experience contributions from the strong interactions whose magni- 
tude has to be known to be able to interpret the matrix elements theoretically. Of great 
interest are weak decays of hadrons which depend on the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements which are fundamental parameters of the Standard Model de- 
scribing the flavor changing weak currents of quarks and play a special role by providing 
a mechanism for CP violation within hadron physics. 

Among the CKM matrix elements which are at present intensively studied in exper- 
iments and theory and which are relevant to CP violation are those related to heavy 
flavored quarks, in particular the b quark. 

At energies less than 1 GeV, the interaction strength increases and non-perturbative 
are required. Placing space and time onto a four-dimensional lattice [1] it is possible to 
calculate matrix elements numerically on a computer within a path integral formalism. 
In the past few years the lattice has become a standard method in QCD calculations. 

To avoid finite size effects, the lattice extent L should be much larger than the wave- 
function size or the Compton wavelength of the particles one would like to describe, 
and the inverse lattice spacing should be much larger than the relevant masses and 
momenta in order to avoid cutoff effects from the finite lattice spacing a. The lightest 
hadrons, the pions, have a mass of around 140 MeV, whereas the B meson has a mass 
of 5.28 GeV and contains a heavy quark with a mass of 5 GeV. Lattice calculations of 
hadrons with heavy and light quarks have to be able to describe physics from this large 



4 



range of scales. 

An remarkable property of hadrons with heavy quarks is that the energy level splittings 
of b hadrons are much smaller than their masses: of the order of Aqcd = 200 — 500 MeV 
or smaller, where Aqcd is the energy scale where QCD becomes non-perturbative. The 
dynamics of the heavy quarks can be accounted for in an expansion in powers of the inverse 
heavy quark mass. This is a basis for effective field theories such as Heavy Quark Effective 
Theory (HQET) (for reviews see [2]) and non-relativistic QCD (NRQCD) [3,4], which can 
be used to simulate heavy quarks directly on the lattice while avoiding large discretization 
errors proportional to powers of the large mass. Besides, heavy quark effective field theory 
predictions can be used to guide extrapolations and interpolations of lattice simulation 
data as a function of the heavy quark mass. 

Simulations with light quarks are computationally expensive and sensitive to the finite 
lattice volume. Therefore one often uses quark masses much heavier than u and d quark 
masses and extrapolates the results to the physical values. A formalism for this can be 
derived using chiral perturbation theory (xPT), an expansion around the chiral (zero 
quark mass) limit describing the low-energy degrees of freedom of QCD, i.e. pions and 
other light mesons and their interactions with other hadrons at low energy. 

If one assumes that chiral symmetry of QCD is spontaneously broken, one expects 
massless Goldstone bosons in the particle spectrum. The physical light mesons can be 
identified with the corresponding Goldstone bosons if explicit breaking of chiral symmetry 
is added by including quark mass terms in the Lagrangian. Hadron interactions give 
corrections to the relations between quark and hadron masses. At low momenta the 
interactions e.g. in vr — tt scattering and n — N scattering are weak and can be treated 
using perturbative methods. 

The aim of this article is to summarize lattice calculations of the light and heavy-light 
hadron spectrum and decay constants the author performed within various collaborations, 
to discuss the systematics of these calculations in some detail, and to compare the results 
to recent other work. To make the article more self-contained, a short introduction to 
some formalisms relevant to the calculations is also given. 

In Section |21 we give a review of the basic theory and phenomena of QCD. A brief 
introduction to the lattice formalism is given in Section El The concepts of effective 
theories which are relevant in the context of QCD calculations, namely NRQCD, HQET 
and Chiral Perturbation Theory, are discussed in Section |3 Having made acquaintance 
with the lattice and the chiral effective field theory for the analysis of the light quark 
mass dependence, we discuss results on the light hadron spectrum and decay constants 
in Section In Section IHl we discuss lattice formalisms for heavy quarks using non- 
relativistic effective theories. In Sections [3 and |H1 lattice results on the 6-light spectrum 
and decay constants are presented. Section IHl contains some conclusions. 
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2 QCD AND THE STANDARD MODEL 



In this section we give a brief overview of QCD emerging from nuclear physics and the 
Standard Model of Elementary Particle Physics. For a recent review about hadron physics 
and methods in QCD see e.g. [5]. 

2.1 Hadrons and the Quark Model 

In the early 20th century it was discovered that atoms possess nuclei which are much 
smaller than the atom and which consist of smaller particles, protons and neutrons, which 
in turn are tied together by the strong nuclear force. We are interested in studying the 
strong force to understand what keeps all matter together. 

Yukawa postulated in 1935 that the interactions between nucleons are due to exchange 
of mesons [6]. From the short range of the nuclear force (interaction radius r ~ 1 fm), he 
concluded that the mesons should have a mass of ~ 1/r ^ 100 MeV. In 1947, the vr meson 
was discovered experimentally with a mass of ~ 140 MeV. K mesons were discovered soon 
afterwards. 

Since the 1970's indications began to accumulate that hadrons, i.e. nuclear particles 
and other particles interacting with the strong force consist of even smaller particles: 
baryons such as protons (p), neutrons (n), hyperons and their like consist of three quarks 
and mesons such as pions, kaons, the J/ip and the T of a quark and an antiquark. 

The masses of proton and neutron are very closely degenerate, and they can be re- 
garded as transforming as an approximately degenerate doublet of 'nucleons' in the fun- 
damental representation of a global SU (2) symmetry group. The corresponding quantum 
number is isospin (J). Within the quark model, the nucleons and pions are composed 
of u and d quarks, both almost massless, which sit in the fundamental representation of 
SU{2). 

The discovery of hadrons with a new quantum number, the strangeness led to 
the introduction of the strange (s) quark. The spectra of strange and lighter hadrons 
can be described using an approximate symmetry between hadrons with strangeness and 
those without. In the 'Eightfold Way' model of Gell-Mann and Ne'eman [7], the three 
quarks can be identified with the generators of an approximate SU{?>) flavor symmetry 
in the fundamental representation. The symmetry group contains isospin as well as the 
quantum number hypercharge Y = B + S, where B is the baryon number. The low- 
lying mesons composed of a quark and an antiquark can be grouped into irreducible 
representations contained in the 33* product representation of SU{?>), the singlet and the 
octet, and baryons consisting of three quarks can be grouped into the octet or decuplet 
representations contained in the product 333. The light baryon states, for which lattice 
results will be discussed later, and their quark content are shown in Table ^ 

The discovery of the A"*'+ gave an indication for the existence of an additional degree 
of freedom, the color. According to the quark model it consists of three u quarks with 
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Quark content 


Octet 


Decuplet 


uuu 




A++ 


uud 


P 


A+ 


udd 


n 


A" 


ddd 




A- 


uus 


S+ 


S*+ 


uds 


SO, AO 


S*+ 


dds 


E+ 


E*+ 


uss 




E*° 


dss 




S*~ 


sss 







Table 1: Baryons and their quark content within the quark model. Electric charges 
are indicated by superscripts. The names of strange baryons can be extended to those 
containing heavy quarks as discussed in sections 17.31 and 17.41 

a totally symmetric wave function. If there was only one species of u quark, this would 
violate the Pauli principle. 

High energy experiments also showed that nucleons have constituents. The existence 
of three colors of quarks showed up in measurements of the cross section for e^e~ decaying 
into hadrons. 

In studies of the structure of the nucleon in high-energy scattering experiments, con- 
stituents ('partons') were found which behaved as almost free particles in the high energy 
limit. Those can be identified with the quarks previously suggested by the quark model, 
and the gluons. 

Important for the understanding of QCD as the theory of the strong interactions was 
the discovery of mesons consisting of heavier quarks. In 1974, a bound state of a heavier 
quark flavor (charm) and its antiquark, the J/ip meson, was discovered. A fifth species 
of quark (6) was found in 1977 with the discovery of the T, abb bound state. It had 
been postulated earlier by Kobayashi and Maskawa to explain CP violation in the weak 
interactions, which had been discovered in 1964 by Christensen et al in decays of neutral 
K mesons into pions. Heavy-light mesons with b quarks are B~^(bu), B^{bd), B (db), 
and B~{ub). The sixth quark postulated by the Standard Model, the top (t) quark, was 
discovered in 1995. 

2.2 Gauge Theories: QED and QCD 

Interactions between elementary particles can be explained as resulting from local, i.e. 
space-time dependent, mathematical symmetries. 

For example, the Lagrangian of free Dirac particles is invariant under global unitary 
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t r ansf ormat ions 



q Vq, 

q ^ qV\ (1) 

where q are Dirac spinors and V G U{1) otV& SU{N). The symmetry group of the strong 
interactions is SU{3) corresponding to three colors of quarks, which sit in the fundamental 
representation of the gauge group. Postulating invariance of the Lagrangian under local 
U{1) or SU{3) transformations requires the introduction of gauge fields coupling to the 
quarks with the covariant derivative D^, and the Lagrangian takes the following form: 

C = q{t^^D^' -m)q + Cg, (2) 

where Cg is the gauge field Lagrangian. The different features of QED and QCD at tree 
level are compared in Table El 



Quantum Electrodynamics Quantum Chromodynamics 



[/(I): V = exp[m(x)] 



electric charge e 
Afj^ photon, neutral 
r — —Ip pp-'^ 

^9 4 P-^ 



5f/(3): V = exp[ia%x)f 
SU{3) generators, [t", 
D^ = d^- igAlt' 
g strong coupling g 
gluon, color charged 



d^A?, 



Table 2: Comparison of QED and QCD tree level Lagrangians. 



Proof of the mathematical consistency of quantum theories with non-abelian gauge 
symmetries [8] in the early 1970's prepared the stage for a description of the strong 
interaction as a gauge theory with a local SU{3) color symmetry, namely QCD. 



2.3 QCD coupling 

If a field theory of point particles is quantized divergent mathematical expressions can 
occur. If the divergences can be absorbed by a redefinition of the parameters of the La- 
grangian the theory is (perturbatively) renormalizable. For non-abelian gauge theories 
renormalizability has been proven at any order of perturbation theory [8]. In general, 
the quantum corrections to a given physical process depend on the momenta of the ex- 
ternal particles. A dependence of physical couplings is introduced on the momentum 
scales used to probe, i.e. look at, the physical system. For dimensional reasons the mo- 
mentum dependence introduces an additional parameter, namely a renormalization scale. 
Then one can derive renormalization group equations, which determine how couplings. 
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masses and Green functions evolve as functions of the renormalization scale parameters 
or, alternatively, on the momentum scales. 

Renormalization can actually be defined in more general terms than in perturbation 
theory only. The most interesting case occurs when limits exist where physical quantities 
such as dimensionless couplings and mass ratios become scale invariant. Renormalizability 
can then be viewed as the possibility to vary the bare parameters such that physics can 
be kept fixed while the distance scales at which the system is probed are changed. In 
perturbative QED, an electric charge appears weaker if it is probed at a larger distance 
because of screening of the charge due to generation of virtual electron-positron pairs. 
Since the QCD gauge group is non-abelian, the gauge fields themselves carry charges of 
the strong interactions and interact not only with the color-charged matter fields, but also 
among themselves. Perturbation theory predicts that the QCD coupling becomes weaker 
due to the self-interactions if it is probed at a smaller distance, i.e. at higher momentum. 
At infinite momentum the coupling vanishes. In 1973 it was shown that QCD has a 
free- field theory asymptotic behavior up to logarithmic corrections [9]. 

The scale dependence of the coupling as = /{Air) is described by the [3 function 

/x^ = 2/3(a,), (3) 
which in perturbation theory can be expanded in powers of a^: 

2/3(a.) = -6og-6,^ + ... (4) 
At one loop, the evolution of Us from a scale /io to a scale /i is given by 

with bo = 11 — ^Nf, where Nf is the number of flavors. If Eq. © is integrated, one can 
re-express the scale dependence (or momentum dependence) in terms of the integration 
constant Aqcd- 

At tree-level QCD without quark masses has no scale parameter. The renormalized di- 
mensionless gauge coupling can be re-expressed using the renormalization scale or the 
momentum and the dimensionful parameter Aqcd- This is called dimensional transmu- 
tation [10]. 

Eq. (jH)) indicates that at larger distance, the coupling becomes strong and diverges at 
/i = Aqcd- Of course, at strong couplings perturbation theory is not valid. For example, 
the physics giving rise to masses of light hadrons is related to energy scales 0{Aqcd) and 
has to be described using non-perturbative methods. 
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2.4 Chiral symmetry 

If tlie quark mass vanislies, i.e. m = 0, tlie Lagrangian (j21) is invariant under cliiral 
transformations 

q exp i'y^a'^t'^q, 

q qexpi-f^a^t", (7) 

where t°- are the generators of the symmetry group. The effect of chiral symmetry on 
the particle spectrum would be degeneracy of positive and negative parity states. This is 
however not found in nature. Instead, the spectrum contains pseudoscalar mesons which 
are extraordinarily light compared to the other hadrons. The lightest mesons are the 
pions TT"*", TT*^, 7r~ which correspond in this picture to the Goldstone bosons of an SU{2) 
chiral symmetry between the u and d quarks broken spontaneously and explicitly due to 
small mass terms. If the strange quark is also considered as almost massless, the chiral 
symmetry is SU (3) and the would-be Goldstone bosons are an octet of vr, K and rj. 

The approximate conservation of the axial vector current A^, the Noether current of 
the chiral symmetry, can be expressed using the PCAC relation. For SU{2), this is 

d.A^; = 2m,5'^^^ (8) 

where vr" is an interpolating field for the pion, ^^(x) = g(x)7^75t''g(x), being a gen- 
erator of SU{2), and rriq is the quark mass. In this section, Euclidean formulae are used 
for the discussion of the topology. PCAC denotes the partially conserved axial vector 
current. Eq. (jHJ can be considered as an operator identity in the sense that it can be 
extended to a relation for matrix elements of Af,. 

The spontaneous symmetry breaking leads to a vacuum expectation value of qq, the 
chiral condensate. This can be calculated e.g. in lattice simulations. In the sum rule 
method of Shifman, Vainshtein and Zakharov [11] the chiral condensate, the gluon con- 
densate {^G'Jj^j^G"''^'^), and condensates of higher dimensional operators are used to pa- 
rameterize nonperturbative contributions to the hadron spectrum and matrix elements 
and to make further predictions using these parameters. For a recent review see [12]. 

The tree-level QCD Lagrangian also has a global chiral U{1) symmetry which has to 
be broken since there are no degenerate states of opposite parity in the hadron spectrum. 
The Tj' meson has the correct quantum numbers for being the Goldstone boson of sponta- 
neous chiral U{1) symmetry breaking, however its mass is too large. Experimentally this 
additional breaking of chiral symmetry is visible e.g. in the decay rate of tt^ 27. The 
symmetry is broken due to an anomaly which can be explained from the noninvariance of 
the measure of the path integral under chiral transformations [13]. The anomaly implies 
that the divergence of the flavor singlet axial vector current = 97^755 is non-vanishing: 

d^A^{x) = Nfq{x), (9) 
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where g, the topological charge density, is related to the color-electromagnetic field 
strengths by 

^(^) = "i^^M.^;:.' (10) 

where G^^, = |e^,^a^G^^ is the dual field strength. With the anomaly, a close relation 
between the topology of the gauge field and the chiral symmetry of the lattice fermions 
is established. An index theorem relates the topological charge to the difference between 
the number of right and left handed zero modes of the Dirac operator D: 

n.-n+ = NfQ. (11) 

If the axial current vanishes at spatial infinity, one finds after integration of Eq. Q also 
that the topological charge is related to the difference of the axial charges at positive and 
negative temporal infinity by: 

j d^xAoit ^oo)- j £xAo{t -oo) = NjQ. (12) 

The topological charge density can also be expressed as the divergence of a current which 
depends entirely on the gauge fields: 



q{x) = d^K^{x 

1 1 
16^ 



where f^^"^ are the structure constants of SU{?)). If this again decays sufficiently fast 
at spatial infinity, the topological charge can also be expressed as the difference of the 
Chern-Simons numbers Nqs of the gauge field at t ^ oo and t —oo, where 

Ncs{t) = j d'xKoit). (14) 

The classical action of the Yang-Mills field is minimized if the gauge fields are self-dual. 
The lower limit to the action of the classical field is proportional to the difference of the 
Chern-Simons numbers. Instantons which are self-dual field configurations with finite 
action corresponding to a change of Chern-Simons numbers thus provide a possible classi- 
cal or semiclassical scenario for the occurrence of chiral symmetry breaking in QCD. For 
reviews about the physics of instantons in QCD see [14, 15]. 

There are indications that chiral symmetry in QCD is restored for a larger number of 
quark fiavors Nf. The restoration is predicted to occur for Nf > 7 by an instanton model 
calculation [16] and by lattice results [17]. 
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2.5 Confinement 



In nature there are no free quarks and gluons. Color charges seem to be confined. Aheady 
perturbation theory predicts a growing of the couphng constant at decreasing momentum 
scale. The one-loop relation is Eq. ® , which has a pole at the scale Aqqd- Perturbation 
theory is of course not valid at very low scales, but it already indicates that there is an 
energy scale where particles are very strongly bound and confinement occurs. Determining 
as at higher energy scales it is however possible to determine Aqcd using perturbation 
theory, using a higher-loop equivalent to Eq.®. 

The spectrum of hadrons actually fits quite well into the picture of an inter-quark 
potential which rises linearly at large distances, corresponding to a fiux string with a 
constant string tension. Using a space-time lattice it could be shown numerically that 
the potential between static color sources, whose world lines are just parallel to the time 
axis, indeed increases at larger separation of the color sources, which implies that QCD 
without dynamical quarks is confining at low momentum scales or large distances. The 
theory with fermions remains confining, although there is a probability for the string to 
break due to the generation of quark- ant iquark pairs. The authors of Ref. [17] find with 
lattice calculations that for Nf < 6 QCD remains confining at low energies, while for a 
larger number of fiavors confinement is lost and there are indications for a complex phase 
structure of QCD at large Nf. 

It seems that confinement can be related to a color-electric fiux tube between color 
charges. A physical picture which could explain this on a classical level is the dual 
superconductor model [18] where the fiux tube enclosed in the vacuum is characterized by 
the condensation of color-magnetic abelian monopoles. The condensation of monopoles 
is indeed found in lattice simulations. One model which could potentially describe the 
physics of the monopoles is the Ginzburg-Landau model. Lattice results indicate that 
the monopoles in lattice QCD may rather be arising from center vortex structures which 
imply a different monopole dynamics from the picture realized in the dual Ginzburg- 
Landau superconductor. A dependence of the abelian monopoles on the method used for 
the abelian projection is investigated in [19-21]. A recent review on lattice results on the 
physics of confinement is given in [22]. 

The previous discussion is related to zero temperature. If the temperature is in- 
creased, QCD becomes non-confining. In QCD without dynamical fermions the corre- 
sponding transition is a first order phase transition, while with more than two species of 
light fermions the transition is found to be a crossover. Recent results for the transition 
temperatures from lattice simulations are [23] 



where Nf is the number of quarks fiavors present in the lattice simulation. The range of 
values quoted corresponds to different simulations. The difference of the Nf = results 





(15) 
(16) 
(17) 



171(4) - 173(7) MeV, 
164(2) - 169(12)(4) MeV, 
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is due to different ways of determining the lattice spacing. For quenched gauge fields the 
transition temperature is thus approximately of the order of magnitude of the scale Aqcd- 

2.6 The Standard Model 

The Standard Model of the strong and electroweak interactions is based on an SU{3) x 
SU{2) X U{1) gauge symmetry with three generations of quarks and leptons as fermionic 
matter fields and a scalar field, the Higgs, which is responsible for the masses of the weak 
SU (2) gauge bosons and the fermions, assuming neutrinos to be massless. The SU{3) 
gauge group gives rise to the QCD interactions, with gluons as gauge bosons. SU{2) x f/(l) 
is the gauge group of the Glashow-Salam- Weinberg model of the electroweak interactions. 
The gauge bosons are photons related to QED and the massive W and Z bosons related 
to the weak interactions. The fermionic matter fields are three generations of quarks 
(m, (i, c, s, t, 6), and leptons (e, i^e, /x, z/^, r, Ur)- 

For a recent review about the status of the Standard Model and new physics see e.g. [24, 
25]. Although providing a remarkbly good description of elementary particle physics, 
the Standard Model has some unsatisfactory features. On the experimental side, recent 
discoveries such as neutrino mixings, the possibility of the existence of undiscovered matter 
in the cosmos and perhaps new results from accelerator experiments [26, 27] indicate 
possibilities where physics beyond the Standard Model may appear. The Higgs particle 
has not yet been found. Further, there are theoretical motivations to search for physics 
beyond the Standard Model which contains a relatively large set of coupling constants 
and masses as input parameters. The values of typical energy scales such as the masses 
of the weak gauge bosons are input parameters of the Standard Model and require new 
physics for explanations. 

A strategy in the search for new physics is to simultaneously measure as many physical 
quantities as possible, test the results for self-consistency using Standard Model formulae 
and search for indications for new theories which can imply the presently known phenom- 
ena. 

The Standard Model parameters relevant to the interactions of quarks are 

• Quark masses 

• Parameters of the weak currents (CKM matrix elements) 

For a review about recent results on quark masses from the lattice and sum rules see 
Ref. [28]. A recent comprehensive review on the Standard Model is given by Ref. [29]. 

2.7 CKM Matrix and Heavy Quark Decays 

Weak flavor changing decays of quarks are investigated to determine elements of the CKM 
matrix [30, 31] which parameterizes the mixings of quark generations in the Standard 
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Model. The corresponding Lagrangian density is given by: 



72 



92 




(18) 



where ul, d^, cl, sl, ti, bi are left-handed quark spinors, a charged weak gauge boson 
and g2 the weak gauge coupling. If there are three generations of quarks, the CKM matrix 
V should be unitary. Its entries are usually parameterized with: 



The CKM matrix elements most relevant to decays or mixings of the b quark are Vcb, Vub 
and Vtd- Since quarks are bound their decays take place on the level of hadron decays. The 
decays of interest for lattice calculations are those to a definite hadron state, the exclusive 
decays. In parallel one also uses inclusive decays, where one sums over all hadrons in the 
final state which can be the end products in a certain quark decay. The high precision 
of measurements of inclusive decay rates is in contrast to some model dependence in the 
theoretical interpretation of the results. is e.g. proportional to the rate for the 

leptonic meson decay B~ l^ui (see Fig. [TJ where / is a lepton (e, or r) and vi the 



corresponding neutrino. Since Vub parameterizes the quark decay b ^ u, its modulus 
can also measured by the semileptonic decays of the B meson B {it, p,uj)lVi. |V^6p 
can be determined e.g. from the semileptonic decay B ^ D. |VJrfp is proportional to 
the oscillation frequency between the mass eigenstates of the B^ — B^ mixing, which is 
described by the left and middle diagrams in Figure 01 in the electroweak theory. 

If there are only three flavors of quarks and no new physics the CKM matrix has to be 
unitary which leads to 9 relations of the matrix elements. Two of the relations relevant 
to b decays are given by: 




(19) 




Figure 1: Leptonic B~^ decay in the Standard Model. 



Vudv:, + v,v:, + VdV* 
Vubv:, + VusV,: + Vudv:, 



0. 







(20) 
(21) 
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To describe CP violation which has been discovered in — K oscillations in 1964 and 
has been established in meson decays in 2001 [32], within the Standard Model, it is 
necessary that some CKM matrix elements have an imaginary part. Eqs. (j^UI) and (PT|) 
describe triangles in the complex plane. Eq. ()2()|1 is in particular used in the so-called 
'standard analysis' of the unitarity triangle shown in Fig. El which focuses on the deter- 
mination of the complex phase of Vtd, related to CP violation within — B oscillations 
via t quark production. The procedure is e.g. reviewed in [33]. 

The CKM matrix parameterizes mixing between quark families. The relative smallness 
of the off-diagonal elements can be taken into account with the Wolfenstein parametrisa- 
tion: 

/ 1 - f A A\\p - tT]) \ 

V = i -X 1-f +0(A^) (22) 

V A\^{1 -p-iv) -A\^ 1 J 

with A = \Vus\- 

With more exact measurements of the matrix elements, a generalized Wolfenstein 
parametrization [34] has become useful: 

/ A + 0(A^) AX%p~^r^)\ 

V=l -X + A^^[l-2{p + ir])] 1-^-^(1 + 4^2) AX^ + 0{\^) \ +0{X'). 
\ AX^l-p-ir]) -AX^ + A^[l-2{p + i7])] \ - A% j 

(23) 

The parameters p and 77 differ from the not barred parameters by Oiy^y. 

'^=^(^"y)'^ = ''0"y)- ^^^^ 

iKifel) and |Vtrf| parameterize the sides of the triangle in the complex plane shown in 



A = (p,Tl) 




C=(0,0) B = (1,0) 

Figure 2: Unitarity triangle. 

Fig. El 
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Rt 



'l-p)^ + v'j\Vtd/V,t\ 



(25) 



One of its angles directly measures the phase of Vtd'- 

Vu = \Vtd\e-'^. (26) 
Lectures reviewing the CKM matrix and CP violation in particular in the context of B 



u,c,t 




Figure 3: Box diagrams describing — B^ mixing in the electroweak theory (left and 
middle) and the weak effective theory (right). 

physics and discussing future directions of research on CP violation in B decays can be 
found in Ref. [35]. A comprehensive review on the status of the elements of the CKM 
matrix is given in [36]. 

Processes at energy scales much less than the W boson mass can be calculated within 
the weak effective theory where interactions mediated by the W or Z particles can be 
described by point interactions. In other words, if the momentum of an exchanged W, Z 
boson is small with respect to its mass M, the propagator can be approximated by 



M2 



M2 



(27) 



The leptonic decay of B mesons (e.g. Fig. ^) can be defined from matrix elements of 
the heavy-light axial vector current going into the vacuum: 



{0\A^{x)\B{p)) = -zfBP,e-^^^, 
where Jb is the B decay constant. The branching ratio for the decay B^ 

GpUiBrnf 



is 



BR{B+ rui) 



Stt 



^1 (. r2\.r 



(28) 



(29) 



where Gp = g2/{^^w) Fermi constant and tb the B lifetime. If is known, 

\Vub\ can in principle be determined experimentally from this decay. There exist only 
experimental upper bounds on Jb and Jb^- Recent experimental measurements of fr, 
from IJ'^^ii exist from CLEO-c [37] with the result 



fo = 223(17)(3) MeV, 



(30) 
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and from BES [38]. Experimental results on fo^ are: 



/^^ = 280(17) (25) (35) MeV [39], and 
/d, = 285(19)(40) MeV [40]. 

In the weak effective theory exclusive semileptonic decays are described by matrix 
elements of currents between hadronic states. They can be parameterized in terms of 
form factors which depend on the momenta of the hadron states. The form factors receive 
contributions from long-distance QCD interactions, and therefore have to be evaluated 
nonperturbatively. This can be done from first principles using the lattice. — 
mixing is described by the matrix element of a current-current interaction, see the third 
diagram in Fig. 01 It is common to parameterize it as f^Bs, where the 'bag parameter' 
Bb quantifies to what extent the matrix element is described by i?-to-vacuum currents 
^ggg 6.g. [41]): 

^ 3 {BO\{hl,dL){hl,d,)\B') 

2{B^\lrf^^,d\0){0\lrf,^,d\B^)' 

In the vacuum saturation approximation, i?^ = 1. The oscillation frequency of the mass 
eigenstates is proportional to the mass difference and related to the current matrix element 
by 

AM, oc IV^IVuI'^bBb. (32) 



The status of the constraints on the unitarity triangle in the p — rj plane including 
results on |Kt6/Kfe|, and B^ — B and B^ — B^ mixings from lattice [42] at the CKM05 
Workshop is shown in Fig. |31 The parameter ex measures CP violation in the neutral K 
meson system. Again, the flavor eigenstates are not eigenstates to CP, rather they are 
superpositions of the CP eigenstates and of the mass eigenstates. Denoting the heavier 
(or longer-lived) one as Kl and the lighter (or shorter-lived) one as Ks, and normalizing 
+ = 1; they can be parameterized as: 

\Ks) = e-^-^\p\K') + qlK")) 

\K,) = e-^-\p\K')-q\K')), (33) 

where ■jl,s are complex, ex is the ratio of amplitudes of decays of and Ks respectively 
into two pions with isospin 0, and is related to the mixing parameters by 



b|2 - 2ReeK 



|p|2+|g|2 l + |g^l2- 



(34) 



The existence of an imaginary part in CKM matrix elements is now well established, 
and apparently the Standard Model can explain the presently known CP violation in 
hadron physics. For recent reviews of the impact of lattice results on the determination 
of parameters of the CKM matrix see [44]. 
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Figure 4: Status of the CKM matrix after the CKM05 Workshop [42]. The shaded regions 
denote the values which are allowed at 95% probability. The dotted line denotes the bound 
on the mass difference of the CP eigenstates. For recent updates see Refs. [42,43] 



The 4-dimensional Minkowski space-time can be substituted with a 4-dimensional Eu- 
clidean lattice by Wick-rotating real time into imaginary time. The Minkowski vacuum 
expectation values of products of fields can be reconstructed from the Euclidean Green 
functions [45]. Matter fields, e.g. quarks, sit on the lattice sites separated by a distance 
a, the lattice spacing. Gauge fields sit on links connecting the lattice sites. Here, a brief 
introduction to the features of the lattice formalism which appear relevant to present 
calculations is given. For detailed introductions see e.g. the textbooks [46] and the review 
articles [47-51]. 

A lattice can be considered as a gauge-invariant regulator. QCD has been shown to 
be renormalizable to all orders in lattice perturbation theory [52]. It is believed that 
continuum QCD can be constructed as a continuum limit of lattice QCD also beyond 
perturbation theory. 



We discuss an introduction into the discretization of a gauge theory, following the dis- 
cussion of e.g. [47]. To arrive at gauge invariant products of fields at different space-time 
points as they occur in actions containing derivatives, one connects them with path- 
ordered integrals of gauge fields 



3 QCD ON THE LATTICE 



3.1 Gauge Fields 




(35) 
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The L{x,y) are parallel transporters which transform under V{x) G SU{Nc) (where Nc 
is the number of colors) as 

L{y,x)-^V{y)L{y,x)V~\x). (36) 

Then q{y)L{y,x)q{x) is gauge invariant. Line integrals over closed loops, L{x,x), are 
called Wilson loops. In particular, TtL{x,x) is gauge invariant. 
The gauge fields on the lattice are represented by 'link' variables 

t/,,^ = Pexp[zga j dz^Al{z)t% (37) 

where the line integration starts at the lattice site which is shifted by one lattice spacing 
in the /i direction from the site x. In the continuum limit this is approximately equal to 

= exp[-igaA1^{x + anf,/2)t''], (38) 

where ^^(x) are the continuum gauge fields, is a unit vector in the /i direction. The 
U fields are parallel transporters between neighboring lattice sites transforming as 

U,,, ^ V,U,,,V^^^, V{x) e SU{N,). (39) 

We use the short form x + fi for x + an^. 

The most elementary product of links along a closed path is the (1x1) Wilson loop, 
or the plaquette, 

W^::' = Tr U.,,U,+,,mU.^,UI,, (40) 
sketched in Fig. El It is used to define lattice actions for the gauge fields. 



Figure 5: Plaquette W^^^^ in the (/i, i^) plane. The corners are separated by the lattice 
spacing a. 



By expanding the gauge links around a ~ 0, 

Ux,^ = 1 - igaA^{x + ^anf,) + 0(0^), (41) 

with A^{x) = f^A^^^x), one finds for the tree-level relation between plaquettes and the 
continuum field strengths: 



ReW'^::' = N,- ^a^ Tr G^^G^'' + 0{a'). (42) 
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Thus one obtains the Wilson (or 'plaquette') lattice action as a discretization of the 
continuum gauge field action of Eq. 



The constant term is sometimes dropped for simplicity. The action (j43|) has lattice spacing 
errors at O(a^). To reduce discretization effects further, actions can be improved. For 
gauge field actions, this consists in adding larger Wilson loops. Improvement can be done 
by removing discretization effects order by order in a (originally suggested by Symanzik 
for scalar field theory [53], and developed into an improvement program for on-shell 
quantities by Ref. [54]), or with a renormalization group method. The perfect lattice 
gauge action [55] is given as the fixed point of a renormalization group trajectory. For the 
renormalization group improved (RG) action [56], the coefficients are determined using 
block spin transformations along a renormalized trajectory. 

Some of the lattice calculations whose results are discussed in the following use im- 
proved actions with a 1 x 2 Wilson loop (W^j^'^) added to the plaquette term: 

The coefficients are given by ci = —0.331 for the renormalization group improved action, 
by —1.4088 for the DBW2 action [57], and by —1/12 for tree-level Symanzik improve- 
ment [58]. The other coefficient is fixed by the condition cq = 1 — 8ci. 

At typical values of /S in lattice simulations, there are large coefficients in perturba- 
tion theory due to tadpole diagrams on the lattice. The perturbative corrections can be 
reduced with 'mean-field' (or 'tadpole') improvement [4]: the gauge links U are divided 
by their expectation value 

U^{x) ^ U^{x)/uo, (45) 

where mq can be calculated in perturbation theory and determined in simulations. A 
gauge invariant definition is the fourth root of the plaquette 

^0 = {w'^::Y'- (46) 

The expectation value of the gauge link in e.g. in the Landau gauge has been used 
in many lattice calculations of NRQCD for the tadpole improvement of the heavy quark. 
NRQCD quarkonium data were found to have a better scaling behavior than using the 
plaquette [207]. 
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3.2 Lattice Fermions 



One can discretize the Euclidean continuum Dirac action by substituting the covariant 
derivatives by covariant lattice differences 

forward ^ ii.(lx = 

backward V*^^ = 

symmetric |(V/,g^ + V*g^) = 

Using the symmetric derivative, one can define the 'naive' lattice fermion action by 

x,y 

■^xy ^ ^ '^1^J,{Ux^fj,Sx+fj,^y Ux_^^^Sx — fj,^y) + UlSxy, (48) 

A* 

which is chirally symmetric if m — > and has O(a^) errors. However, the naive discretiza- 
tion leads to a flavor multiplication, the so-called 'doublers'. The fermion propagator is 
a periodic function of the momentum /c^ Considering for example the case m = 0, one 
finds that the free propagator has next to the pole at /c^ = which also exists in the 
continuum, poles at = n/a. The states corresponding to the poles are degenerate and 
have finite energies. So one has 16 species of fermions, which occur in pairs of opposite 
chirality. The occurrence of doublers is a problem with a quite general class of fermion 
actions. According to the Nielsen-Ninomiya theorem [59], for any lattice fermion action 
which is translation invariant, local and has a hermitian Hamilton operator, there is an 
equal number of right- and lefthanded fermions corresponding to each quantum number. 
Wilson's solution to the doubling problem is to add a term of the form 

Sw = 0''^'y^^QxDw,xy<ly, 

X 

Dw,xy = - —{Ux,tiSx+tj.,y + Ul_^^^5x-^i,y - 2Sxy) , (49) 

to the action. With the Wilson term the doublers obtain masses which remain finite in 
lattice units: m^a 7^ 0, i.e. ~ 1/a if a — > 0. The exact behavior of the doubler masses 
depends on the choice of r. The most common value in lattice simulations is r = 1. 

To normalize the term proportional to 6x^y in the action to one, it is common to 
parameterize the mass m in terms of the hopping parameter k. Setting r = 1, the action 
is: 

Y.'^xD'xr^y, 

X 

(50) 



^{Ux^fiQx+n Qx)i 

-Mx-ul 



2a 



-^l,^lQx-^J,) 

'X- 



(47) 



{Ux^fiQx+ii Ux — a,iiQ.x- 



Sf 



Y)hopp 
xy 
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with 



am 




(51) 



and 



q-x = a 



(52) 



The Wilson term introduces a shift in the quark mass parameter such that the massless 
free Lagrangian corresponds to k, = Hc = 0.125. 

In the interacting case, the critical hopping parameter is renormalized and the bare 
quark mass is given by 



also called the vector ward identity (VWI) quark mass. Now chiral symmetry is broken 
at a 7^ 0. Discretization errors of 0{a) occur and manifest themselves on dimensional 
grounds as 0{am), 0{ap) effects, where m is the quark mass and p the quark momentum 
which is in a hadron typically of 0{Aqcd)- The k value corresponding to zero quark mass, 
Avc, can be determined using the hadron spectrum. The pion mass is usually associated 
with explicit chiral symmetry breaking, therefore one defines Kc as the hopping parameter 
where the pseudoscalar meson mass vanishes. With a tadpole-improved Wilson action, 
— * U^/uq the hopping parameter is redefined to k = kuq. In the tree-level mean field 
improved theory the critical hopping parameter is 



such that K = 0.125. 

Another definition of quark mass used in lattice calculations [60,61] is to use the axial 
vector ward identity (AWI) Eq.(jHl): 



where Af^{x) is a lattice version of the axial vector current, which includes improvement 
terms where necessary. 0(a) errors can be removed from the quark action using the clover 
term 




(53) 



Kc = 1/(8mo), 



(54) 




(55) 




(56) 



where 0"^,^ = |[7^,7iy], and G^u denotes now a discretized version of the field strength 
tensor using four neighboring plaquettes [62]: 




(57) 
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where the sum goes over all four plaquettes in the (/i, u) directions which are connected 
to the lattice site x. The W denote the untraced plaquettes, and a, b are color indices. 
The coefficient of the clover term can be calculated in perturbation theory. In many 
of the calculations whose results are discussed later, csw is at tree-level with tadpole- 
improvement. Most recent calculations use a non-perturbative determination of the clover 
coefficient [63], i.e. the action is 0(a) improved to all orders in perturbation theory. 

For the calculation of matrix elements of currents, improvement of the field is also 
necessary. For 0(a) improvement this can be done by changing the quark fields by a term 
of the form 

q{x) — > Vl + bamqi^x), (58) 

where b can be determined nonperturbatively. The Fermilab normalization for massive 
quarks, Eq. ()2()4|1 of Section 16.3^ corresponds to 6 = 1 at tree-level. In our lattice cal- 
culations discussed in Section |H1 we used tadpole improvement, i.e. b = uq. The current 
operator itself has to be improved as well by adding a higher dimensional operator hav- 
ing the same symmetries with the appropriate coefficient. The improvement of currents 
due to higher dimensional operators in effective theories is very similar and discussed in 
Section KT^ 

A spontaneous breaking of flavor and parity symmetries in lattice calculations with 
dynamical Wilson fermions has been conjectured by Aoki [64]. With two flavors, the 
SU{2) symmetry breaks down to a residual U{1), and there are two massless Goldstone 
bosons and one Goldstone boson with mass proportional to the lattice spacing. Ref. [65] 
finds, using chiral perturbation theory, that for certain parameter values of the chiral 
Lagrangian an Aoki phase might exist close to the continuum limit, while the authors 
of Ref. [66] find indications with a dynamical lattice calculation that the Aoki phase is 
contained only within a region of relatively strong coupling. 

Another standard technique to reduce the number of doublers is to use staggered 
fermions, which is the method chosen in some of the lattice calculations discussed in the 
sections EJ [7| and IHl Staggered fermions are obtained from a spin-diagonalization of naive 
four-component fermions: 

Qx = IxXx , Qx = Xxll ' (59) 

with 7a, = 71^72^73^74'', where x = {xi,X2, X3, Xi). In the tree-level Lagrangian the spinor 
components of the x fields are disentangled, and chosing one component one obtains the 
staggered fermions. The staggered action is at finite lattice spacing invariant under a 
vector U{1) symmetry, and for zero mass there is an axial lattice U{1) symmetry which is 
in the continuum limit flavor non-singlet. For a discussion see e.g. the review article [47], 
and for a detailed explanation of staggered fermions and their symmetries [67]. In the 
continuum limit there is a SU{4:)l x SU{4:)ji x U{1)l x U{1)r flavor symmetry, where 
the axial part of the chiral U{1) has the usual anomaly. Discretization errors of staggered 
fermions are O(a^). Recent calculations use staggered quarks which are improved by 
adding higher dimensional operators, see e.g. [68]. 
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Continuum symmetry: 


Lattice realization 




Gauge symmetry 


exact 




Euclidean invariance: 


recovered for a 


0,V ^oo 




naive fermions: 


remaining, O(a^) errors, 15 






doublers 




Wilson: 


broken at 0(a), doubler 






masses lifted 




staggered: 


remaining, O(a^) errors, 3 


Fermion chiral symmetry: 




doublers 


clover: 


0(a) breaking terms removed, 






doubler masses lifted 




Ginsparg- Wilson: 


O(a^) errors, single chirality 






possible, invariant under lat- 






tice chiral transformation at 






a / 


U{1)l X U{1)r anomaly: 


exists at finite lattice spacing for Ginsparg- Wilson 




continuum anomaly reproduced for a ^ for all 




five formulations 





Table 3: Symmetries of lattice fermions. 



The Ginsparg- Wilson relation [69] expresses a modified chiral symmetry. If a lattice 
fermion Dirac operator D obeys the equation 



75L' + L'75 = aD-f^D, 
the corresponding action is invariant under a modified chiral transformation [70] 



(60) 



exp 
ip exp 



ze-T^^,(l - -D) 



(61) 



where D is the massless gauge covariant lattice Dirac operator, and one speaks of lattice 
chiral fermions or Ginsparg- Wilson fermions. For a — >• 0, it corresponds to the usual chiral 
symmetry. 

Under quantization, the measure in the path integral is not invariant under this trans- 
formation, and the lattice axial U(l) symmetry has an anomaly with an index theorem 
similar to the chiral U(l) anomaly known from continuum QCD. 

Chiral lattice fermions are e.g. domain wall fermions in the limit of infinite 5th di- 
mension [71], overlap [72] and fixed point fermions [73]. Their popularity in QCD lattice 
simulations is increasing. Renormalization effects are simplified due to the chiral symme- 
try. Ginsparg- Wilson fermions are on-shell 0(a) improved, and due to the difference in 
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the eigenvalue distribution of the Ginsparg- Wilson Dirac matrix compared to the clover 
Dirac matrix it is possible to simulate at small quark masses. 

A brief summary of lattice fermions and their symmetries is given in Table El An im- 
portant remark is that chiral fermions are automatically 0(a) improved, see the discussion 
in [74]. Chiral symmetry and chiral fermions on the lattice are reviewed in [75,76]. 

The Schrodinger functional [77] is the partition function of a Euclidean system with 
certain boundary conditions at the initial and final time. It provides a scheme to calculate 
renormalized quantities as a function of the lattice extent, which is a natural length scale in 
lattice calculations. In practical calculations one often uses periodic boundary conditions 
in space and Dirichlet boundary conditions in time [77]. 



3.3 Hadron Correlation Functions 

Green functions can be calculated by evaluating the path integral over the lattice degrees 
of freedom numerically. For example, a two-point function of a field 0{x), which can for 
example be composed of gauge and fermion fields, is given by 

{0\O\x)O{0)\0) = ^ J Pf/PgPgO^(a;)O(0)e-^«[^]-^/[^'^l (62) 

Ideally, these calculations are done at various values of the lattice spacing, and the con- 
tinuum estimate is obtained by extrapolating as a function of a to a —> 0. In practice, 
some lattice calculations are performed only at one or two values of a, in which case a 
continuum limit cannot be taken, and the discretization effects have to be included into 
the estimate of systematic errors. With NRQCD calculations, higher dimensional oper- 
ators are included as discussed in Section Therefore an a — * extrapolation cannot 
be done. Calculations at several values of a then serve to determine the discretization 
uncertainty. 

Correlation functions are expectation values of time ordered products of hadron oper- 
ators: 

C{t,x) = (0|T[O^(t,f)O(0)]|0). (63) 

One can now derive how masses can be extracted from lattice correlation functions (see 
e.g. [47]). We are interested in correlation functions of hadron operators for states with 
definite spatial momentum p, 

C{t,p} = (0|T[Ot(t,p)O(0)]|0). (64) 

The momentum projection is obtained by Fourier transforming, i.e. summing over all 
points of the spatial volume at the 'sink': 

C{t,p) = ^(0|T[O^(t,x)exp(-ip-f)O(0)|0) 

X 

^ {0\T[O\t,p}Om0) (65) 



25 



For the moment we consider a Euclidean 4-dimensional lattice with spatial volume V and 
infinite extent in time. The time evolution of the operator is in the Heisenberg picture 
described by exponentials of the Hamiltonian H, and the correlation function can be 
re-written as 

C{t,p) = (O|e^*0t(o,p)e-^*0(0)|0),t>0 

= (0|OV^*O|0). (66) 

Inserting a complete set of energy eigenstates \n) with energy En and using relativistic 
normalization, this gives 

p — Ent 

C{t)=J2mO\n)\'^^, (67) 

n " 

i.e. for sufficiently large times the correlation function decays exponentially with the 
ground state energy. For p = this is given by the mass. If the lattice is finite in 
the time direction, and depending on the fermion action, there may be a modification to 
the exponential behavior. 

To reduce discretization errors and to improve the resolution in the time direction, 
some calculations use anisotropic lattices with a finer temporal than spatial lattice spacing. 

Local interpolating operators for relativistic quarks to project onto pseudoscalar, 
scalar, vector and axial vector mesons are 

O^^x) = q''ix)Tq\x) (68) 

where F = 1, 75, 7/^, 757;^, and a and (3 are fiavor indices. A correspondence of lattice 
operators to continuum quantum numbers is given in Tabled 





0~ 


0+ 1" 1+ 


F 


75,7o75 


l,7o Ik 7fc75 



Table 4: Meson operators corresponding to a set of ground state quantum numbers. 

Orbitally excited mesons (e.g. = 0"*", 1"*" and 2"*"), can also be constructed using a 
combination of spin and derivative operators, see e.g. Table |H1 In Table |H1 also examples 
for baryon operators are given. Further examples are in Sec. El 

The correlators can contain contributions from all states with the same quantum 
numbers. To improve overlap with the ground state, one uses smearing, i.e. spatially 
extended meson operators. For example, working in a fixed gauge, one can use wave 
functions motivated from quark model considerations. For [79,80], hydrogen- like smearing 
functions were used. 

0(r) = exp(— r/ro) (69) 
'^W = 7^fl-7^V^P(-^/(2^o)) (70) 
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to optimize overlap onto the ground and excited state S wave state respectively. A gauge 
invariant technique is e.g. Jacobi smearing where the meson shape is approximated with 
the Green function of a Klein-Gordon equation for a scalar field coupled to gluons in a 
hopping parameter expansion. For an explanation of the method see e.g. [81]. 

A measure for the quality of the projection to a single state is the 'effective mass'. If 
the correlation function in Euclidean time is C(t), then the effective mass is given by 



In 



C{at + 071.4) 
C{at) 



(71) 



Examples of effective masses of B meson correlators are given in Fig. 



3.4 Dynamical Fermions and Quenching 

In full QCD, the partition function is an integral over gauge and fermionic fields 
J T>UT>'qDq. To decrease computational expenses, lattice calculations are often done in 
the quenched approximation, i.e. the path integral is calculated setting the determinant 
of the quark matrix + m to one. 

Simulations with two flavors of light fermions have been performed for a range of 
lattice spacings and various gluon and quark actions. Simulations in particular at small 
quark masses can be computationally very expensive. An estimate of how the computer 
time T in simulations including two flavors of clover fermions scales with quark mass, 
lattice size L and lattice spacing a based on CP-PACS calculations is given by [82]: 

where {nips/mpY is a rough measure for the quark mass (see the section on chiral per- 
turbation theory). 

It is computationally significantly less prohibitive to simulate dynamical staggered 
instead of Wilson fermions at relatively light quark masses. Because of the flavor doubling, 
dynamical simulations with staggered fermions use the square root (for Nf = 2) or the 
fourth root (for Nf = 1) of the fermion determinant. Although Ref. [83] finds that the 
square root of the free staggered fermion operator is non-local in the continuum limit, 
Ref. [84] finds a representation of the fourth root of the free staggered fermion matrix in 
terms of a local operator and discusses generalization to interacting fermions. 



3.5 Systematic Errors 

To use lattice results for phenomenology, it is necessary to estimate systematical errors 
as accurately as possible. The most important sources are: 

• Finite lattice spacing 
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• Finiteness of lattice volume 

• Quenching (unphysical number of dynamical quarks) 

• Extrapolation to physical quark masses 

As mentioned before, if the momentum scales are sufficiently smaller than the inverse lat- 
tice spacing, in QCD with relativistic actions one can apply improvement in an expansion 
in powers of a. The Schrodinger functional [77] method utilizes the finite lattice extent as 
a length scale of the theory, and has been proven to have a well-defined continuum limit. 
The results in the infinite volume are obtained using a step-scaling technique, i.e. rescaling 
the lattice size at fixed a in a first step and rescaling L/a, making the lattice larger and 
finer in such a way that the renormalized coupling is kept fixed. For an explanation see 
e.g. [78]. Quarks whose masses are larger than one in lattice units, and where discretiza- 
tion errors proportional to powers of the mass are expected to be large with conventional 
formalisms, can be simulated using lattice formulations of effective theories. If these are 
non-renormalizable, a continuum limit is not possible but instead one has to find a lattice 
spacing with minimal systematic error. The reliability is greater if there is an interval of 
lattice spacings where the results scale as a function of the lattice spacing. A theoretical 
handle to calculate extrapolations from the quark masses used in the lattice simulations 
to the M, d quark masses is provided by chiral perturbation theory, discussed in the next 
section. Chiral perturbation theory also provides a theoretical handle to calculate finite 
volume effects. Extrapolations of heavy masses are usually supported using predictions 
from HQET. Since the formulae partly depend on unknown parameters, the most reliable 
way to reduce the extrapolation uncertainty is to simulate at quark masses which are as 
close as possible to the physical quark masses, which requires, as mentioned above, effec- 
tive theories if the quark mass is very large. This will be discussed in Section |3J Before 
this, we look at various methods to set the scale, which is in particular closely associated 
with the discussion of quenching effects. 

3.6 The Static Quark Potential 

The potential between a static quark-antiquark pair (see also Section IU!T|) parameterizes 
confinement at large distances as well as the strong coupling as at short distances. Much 
interesting phyiscs is associated with it. For a review see e.g. [92]. Here we are just 
interested in the practical aspects. The potential contains length scales, and since it has 
small discretization errors (O(a^) or higher for quenched gauge fields), it is important 
for determinations of the lattice scale. Further it is related to the strong coupling 
which is a fundamental constant of nature and of interest here in particular for use in the 
perturbative matching calculations of matrix elements. Some general features and results 
are described in the following. 

Flux tube like field configurations in QCD give rise to a linearly growing part of the 
potential. The zero point energy of the string vibrations gives rise to a 1/r contribution 
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(see Refs. [14,15] in [93]). Using the string picture, the QCD potential is predicted to be 
of the form 

nr) = ^r + ^-^ + o(l), (73) 

where r is the distance between the quarks, a the string tension and \x a mass parameter. 
A usual parameterization of the potential is 

\/(r) = I/q - - + or. (74) 
r 

If the quarks mass is heavy but finite, the force between them also depends on the spins 
of the quarks and their angular momentum. Derivations of the spin-dependent correction 
terms to the static potential are given in Refs. [94, 95] . 

On the lattice, discretization errors of the static potential can be reduced using im- 
proved actions. Breaking of the rotational invariance of the potential is a signal of dis- 
cretization errors. Comparison of the rotational invariance with a renormalization group 
improved gauge action with respect to the Wilson gauge action is demonstrated on dy- 
namical lattices at an inverse lattice spacing of ~ 1 GeV in [96]. The result is also 
given in Fig. IHl With dynamical gauge fields one expects a flattening of the potential at 
larger distances since the string can break due to quark- ant iquark pair creation. An expla- 
nation why flattening is not visible in the CP-PACS result of Fig. IHl could be insufficient 
overlap of the operator projecting onto the static quark-antiquark pair to calculate the 
potential with the static-light meson-antimeson pair which is created as an end product 
of string breaking. 

The Coulomb part can be used to define the coupling in the V scheme [97,98] via the 
relation 

I/(g^) = -4vr^^^, (75) 

at large g^. Cp = 4/3 for SU{3). Their definition relates the strong coupling directly to 
the potential, while higher order corrections do not appear. The procedure is expected 
to minimize higher order corrections to other relevant physical quantities as well [99]. 
The scale of the running coupling is then also directly related to the momentum transfer 
between the static quarks. 

It can be argued that the plaquette is dominated by short distance effects and non- 
perturbative contributions are small. Therefore it can be expanded in a power series in 
the coupling. E.g. using the V scheme one has: 

-\nW'>^\Nf) =f2^^'^i^v'\llxi)l (76) 
1=1 

where gj^^i is the average gluon momentum in the perturbative calculation of the plaque- 
tte, explained below. Truncating after the second term one obtains the definition for the 
coupling in the plaquette or P scheme [100], such that the plaquette has contributions 
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Figure 6: Heavy quark potential for Nf = 2 with Wilson gauge and quark actions (upper 
left), RG gauge and Wilson quark actions (upper right), Wilson gauge and clover quark 
actions (lower left) and RG gauge and clover quark actions from [96]. The clover quarks 
are mean-field improved. Different symbols denote the axes in units of the lattice spacing 
along which the potential has been measured. 

only up to second order in perturbation theory. The P scheme has been used e.g. in the 
perturbative calculations for the decay constant calculation presented in [101]. 

To calculate perturbative corrections to a given physical process one can use the av- 
eraged exchanged gluon momentum, q*, as a scale for a.^ [97, 102] to minimize higher 
order perturbative corrections, q* depends upon the physical quantity / which is sup- 
posed to be calculated. 61, the one-loop correction to J, is an integral over a function 
^ of the momentum q. One can choose the scale for the running coupling such that the 
one-loop correction is equal to the value with ay for the dressed gluon integrated over 
the momenta [102]: 
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Using the one-loop formula for the running of ay from g* to g, this gives 

' - i'l^'i) ■ 

Other lattice schemes, e.g. the Schrodinger functional scheme are also used in recent 
determinations of the strong coupling constant or alternatively h^qcD from the lattice. 
For comparison to results from experiment for example one usually converts the lattice 
results to the MS scheme. The first important scheme for the strong coupling constant 
in lattice calculations is the bare lattice coupling ah = g^/{27i) with g given in Eq. 
Since perturbative calculation using this scheme usually suffer from large higher-order 
corrections from tadpole diagrams, it is more common to use 'boosted' perturbation theory 
with as = afj/uQ. In Refs. [89, 103] some renormalization constants calculated in the 
boosted scheme are employed. 

We now give some recent results for the values of Aqcd and as determined using the 
lattice. The ALPHA collaboration obtains [104] 

^QSD,Nf=o = 238(19) MeV; 

^Wd,n,=2 = 245(16)(16) MeV. (79) 
and the UKQCD-QCDSF collaboration [112] 

KSd,n,=o = 259(1)(20) MeV; 

^Wd,n,=2 = 261 (17) (26) MeV, (80) 

where the first error is statistical and the second systematical. The UKQCD-QCDSF 
result corresponds to a value of as at 5 flavors of 

«S=5(Mz)=0.112(l)(2). (81) 

That result can be compared with the world average including experimental measurements 
and lattice results as input [105] of 

cSf=5iMz) = 0.1182(27), (82) 

and a recent lattice calculation using 2 + 1 dynamical quarks [106] : 

af^^^,{Mz) = 0.1170(12). (83) 

Recent lattice determinations of as and Aqcd are reviewed in [107]. 

In many lattice calculations one determines the lattice spacing from length scales from 
the potential. For example one uses vq, related to the interquark force [108] with 



2 ^ 



= 1.65. (84) 

r=ro 
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It can be calculated on the lattice with high precision; determinations using quenched 
Wilson gauge fields are given in [109,110]. However for the string tension a and tq it is 
not quite simple to assign experimental values. Relying on potential models to extract 
the physical result from experimental quarkonium level splittings, one finds values of 
tq = 0.49 — 0.5 fm. One can also determine tq from the lattice. With a two-fiavor 
lattice calculation using the nucleon mass as input, Ref. [112] quotes tq ~ 0.467 fm. With 
Nf = 2+1 and using T 25*— IS" and IP— IS" mass splittings to set the scale, Ref. [113] finds 
values of ro in the range 0.462 — 0.467 fm. Unless noted otherwise, we use an experimental 
value of 0.5 fm when we use tq to convert lattice results to physical units. For the string 
tension a, usually experimental values of ^/a = 427 or 440 MeV are assumed. The relation 
between the potential parameters tq and ^/a is relatively independent of the gauge action 
and discretization errors. A dependence on the fiavor number can be observed e.g. in the 
results for r^^^fa (see TableE}. Ref. [114] finds indications for string breaking in two-fiavor 
staggered QCD. A recent study of string breaking in two-fiavor QCD with a sea quark mass 
around the strange mass and at zero temperature [115] finds clear indications for string 
breaking at a distance of r ~ 1.2 fm, which corresponds to roughly 6 GeV~^. One should 
note that in the CP-PACS result shown in Fig. IHl there is no indication of a fiattening 
at this distance. Indications for string breaking in two-fiavor QCD at finite temperature, 
however below the deconfinement phase transition have been found by Ref. [116]. 



Ref. 


/3 




Wilson, A^/ = 


[110] 


5.7 1.160(15) 


[110] 


6.0 1.175(5) 


[117,118] 


6.5 1.199(25) 


1-loop SI, iV/ = 


[119] 


1.160-164 


RG, A^/ = 


[111] 


2.187 1.158(18) 


[111] 


2.575 1.158(8) 


RG, A^/ = 2 


[111] 


1.80 1.143(24) 


[111] 


1.95 1.142(16) 


[111] 


2.10 1.129(5) 


[111] 


2.20 1.144(5) 


1-loop SI, A^/ = 2 + 1 


[119] 


1.114(4) 



Table 5: tq^Iqiy Wilson, renormalization group improved and 1-loop Symanzik improved 
(SI) gauge actions. The A^^ = 2 values for the RG action are from the data sets with the 
smallest sea quark masses. 
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3.7 Setting the Scale 



Masses and decay constants coming out of a simulation are at first dimensionless numbers 
in units of the lattice spacing. The value of the lattice spacing is determined by calculat- 
ing a suitable quantity aM on the lattice and adjusting the corresponding dimensionful 
quantity M to its physical value. If the calculation is free of systematic errors such as lat- 
tice spacing, finite volume and quenching effects, using any quantity should give the same 
result. In practical calculations all of those errors can occur. Then, 'suitable' means that 
systematic errors of the quantity used to set the scale and the quantity that is supposed 
to be calculated cancel as well as possible. 

Length scales from the static potential often used in lattice calculations are the Sommer 
parameter ro and the string tension, discussed in the next section. Hadronic quantities 
used to set the scale are nip, the nucleon mass, the decay constants /jr and fx, and 
charmonium and bottomonium level splittings. 
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Figure 7: Ratio of inverse lattice spacings from various physical quantities to the inverse 
lattice spacing from tq on quenched lattices. Results are from [85-89, 109, 111]. The 
overbar in the legend denotes an average over spin-orientations. 

As mentioned in the discussion of the light hadron spectrum, one expects a system- 
atic difference between the quenched and the real world. For example, in the quenched 
approximation the coupling runs differently, and the string between quark and antiquark 
does not break. 

This is in particular of interest for heavy-light mesons since their physics contains 
several rather different momentum scales. 

In Fig. 13 we show examples for the discrepancy between lattice spacings determined 
from different physical quantities in the quenched approximation. In the quenched ap- 
proximation the p does not decay, there is however a systematic uncertainty in the chiral 
extrapolation. The mp values in these plots result in some cases from simple linear fits 
in the quark mass, which might turn out as not an appropriate parameterization for very 
light quark masses (see Fig. ^2 in Section f5.2|l . To date many calculations with h and 
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Nf = 2iight + Istrangc flavors of dynamical quarks (right), a is determined using the T' — T 
mass splitting. 

light quarks do not reach light quark masses much below the strange, which gives rise to 
a systematic uncertainty which has to be taken into account. 

From HQET, typical momentum scales in B mesons are expected to be of similar 
magnitude as in light-light mesons, namely 0{Aqcd)- The use of rrip, /^r, and fx to set 
the scale in heavy-light meson calculations is motivated by the expectation that quenching 
errors would cancel optimally. In calculations of the heavy-light meson spectrum one 
indeed finds good agreement with experiment, if nip is used to set the scale. Using ro, one 
tends to find a smaller lattice spacing than e.g. using rrip. However, in recent quenched 
calculations with chiral and chirally improved fermions the ground state spectrum of light 
hadrons agrees well with experiment already in the quenched approximation, depending 
on the procedures how the lattice spacing is determined and the physical quark masses 
are fixed [120,121]. Their preferred quantity to set the scale is ro, but if their calculations 
are confirmed it would imply a rather good consistency of a values determined from ro, 
rrijv and rrip. 

With two flavors of sea quarks, usually a better agreement between scales is found 
than demonstrated in Fig. [3 Using Wilson gauge fields and two flavors of 0(a) improved 
clover sea quarks, Ref. [122] quotes an agreement of scales from rrip, fx and ro = 0.5 fm. 
With renormalization group- improved glue and tadpole- improved clover quarks, lattice 
spacings from rrip and ro agree in the chiral limit of sea quarks except for very coarse 
lattice spacings [111]. Agreement of a from rrip and tq is also quoted in [123]. However, 
in the two flavor calculations of [86, 124] (Wilson gauge fields, staggered sea and clover 
valence quarks) and of [87,111] (renormalization-group- improved gauge fields and tadpole- 
improved clover sea and valence quarks) at a ~ 0.5 GeV~^, there is a ~ 20% discrepancy 
between lattice spacings from rrip and Xb mass splittings calculated using NRQCD. 

With two flavors of light and one flavor of strange dynamical quarks, using a 1-loop 
Symanzik O(a^) improved gauge action and a tree- level tadpole O(a^) improved stag- 
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gered sea quark action, Ref. [90] finds agreement of a variety of physical quantities with 
experiment (see Fig. A recent calculation with 2 + 1 flavors using 0(a)-improved 
clover quarks finds that the continuum extrapolated strange meson masses agree with 
experiment [91]. 

4 EFFECTIVE FIELD THEORIES 

In many situations one can obtain quite a good description of physical phenomena if 
only interactions determined by energy scales which are low with respect to a cutoff scale 
are considered dynamical, and the interactions at higher energies are contained in the 
coupling constants of the low energy theory. 

If the low energy interactions are weak, they can be described using perturbative 
methods. An example for this is chiral perturbation theory which is an expansion in terms 
of masses and momenta of the lightest particles in QCD, namely the pions, which can be 
thought of as would-be Goldstone bosons of the spontaneously broken chiral symmetry 
of QCD. Certain constraints on the interactions can be derived from chiral symmetry. 
For example, the pion self-interactions vanish in the limit of vanishing pion mass and 
momenta. 

In hadrons with heavy quarks, mass splittings are often ©(fewAgc^i), i.e. much smaller 
than the mass of the heavy quark itself. Effective theories. Heavy Quark Effective Theory 
(HQET) and Nonrelativistic QCD (NRQCD), were developed to describe the interaction 
of soft gluons in heavy mesons with momentum scales much smaller than the meson 
mass. The heavy mass acts as a cutoff. The low-energy degrees of freedom are gluons 
and light quarks interacting via the QCD gauge coupling as in QCD, and modified heavy 
quark fields. The large momentum scale is given by the heavy quark mass, and low 
energy physics is again an expansion in terms of a small momentum scale. In NRQCD, 
this is a power of the heavy quark velocity in the meson, whereas in HQET it is the 
heavy quark recoil due to the interaction with soft gluons. Since the low-energy physics 
is non-perturbative, the expectation values of the corrections have to be calculated non- 
perturbatively, for example on the lattice. The effect of the high energy interactions 
modifies the coupling constants of these terms. Heavy meson or baryon chiral perturbation 
theory describes the dynamics of the low-energy interactions of a heavy hadron with a 
pion in an expansion around the limits of infinite hadron mass and small momenta, using 
similar modified fields for the heavy degrees of freedom as HQET. 

Constraints on the precision of a lattice calculation with heavy and light quarks result 
from the finite volume of the box and from the finite lattice spacing. On one hand is the 
particle mass the inverse of the correlation length ^. If discretization artifacts are supposed 
to be small, ^ should be large compared to the lattice spacing. As discussed in the previous 
section, discretization effects are in conventional lattice formalisms proportional to powers 
of the quark mass and momenta in lattice units. Using effective theories on the lattice 
it is possible to re-formulate physics in terms of the energy splittings of hadrons, where 
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energy scales and the associated discretization errors are smaller. On the other hand the 
volume has to be much larger that the wave functions of the particles to be described, or 
their Compton wavelength, such that finite volume effects are small. 

In the following table we show schematically the various relevant scales in a lattice 
simulation of hadrons with a b and a.u or d quark together with the boundaries achievable 
on present lattices with dynamical quarks: 



Nature 


= 140 MeV 
nib 5 GeV 
AE{B** - B) ^ 400 MeV 


Ideal 


„ . physical quantities , ^ 
nnite size , . ,^ cutott 
lT-i\ ^ (masses oc 1/4 <C , 

and energy splittings) 


Lattice 
Reality 


< 100 MeV 10 MeV < < 0.5 - 1 GeV 0-^=2-4 GeV 
(L = 2 - 3 fm) 



4.1 Chiral Perturbation Theory [xPT) 

The chiral symmetry of massless QCD is assumed to be spontaneously broken. As a 
consequence one would expect massless Goldstone bosons in the spectrum. One can 
identify the physical pions with the Goldstone bosons if small quark mass terms are 
included in the QCD Lagrangian resulting in a small explicit breaking of chiral sjTiimetry 
and a finite but small pion mass m^r. Typical momenta of low-energy interactions of pions 
will be 0(m7r), which is the momentum scale that dominates also the integrals over the 
pion loops. 

If hadronic interactions at low energies are to be understood within this scenario, it is 
of interest to be able to describe the interactions of the Goldstone bosons of spontaneously 
broken chiral symmetry, within a field theoretical framework. 

The Goldstone boson couples to the symmetry current with the pion-to- vacuum matrix 
element given by 

{Q\A,{x)\iT{p)) = -ip,Ue-'P-^, (85) 

where is the axial vector current. A relativistic normalization is assumed. The zeroth 
component of an axial vector current is an interpolating operator for pseudoscalar states 
in general. For Goldstone bosons the vacuum matrix element of d^A'^ vanishes, since 
the Goldstone boson mass is zero. The matrix element of d^A^ between two hadron 
states is proportional to the pion momentum squared, and therefore also has a vanishing 
zero-momentum limit if the symmetry is spontaneously broken. In the case of explicit 
symmetry breaking it is proportional to the mass-squared of the on-shell pion. 
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The idea is therefore to construct a Lagrangian with the required symmetries which 
can be used to describe the low-energy interactions in terms of hadron fields and hadron 
currents. One can formulate the low-energy interactions in terms of an expansion in 
powers of the pion mass and momenta p = 0{m.^), the so-called p expansion of chiral per- 
turbation theory [125-127]. Beyond tree level, quantum corrections generate new effective 
interaction terms with new coupling constants. To make a perturbative calculation in the 
framework of the p expansion meaningful, it should be possible to determine how each of 
the perturbative loop diagrams contributes to the various powers of the p expansion. In 
the systematic expansion of Ref. [129], n-loop corrections to amplitudes are 0{p^"') with 
respect to the tree-level amplitude. The lower order terms which are generated by each 
diagram can be mapped onto redefinitions of couplings of already existing terms in the 
Lagrangian, while higher order terms predict corrections to the momentum dependence 
of the amplitude under consideration. 

To give meaning to the new interactions and to be able to assign the correct order to 
them a careful choice of the renormalization procedure is necessary. 

The resulting formalism can be applied to low-energy hadronic interactions, for ex- 
ample pion-pion and pion-nucleon scattering and many others which are of interest for 
determinations of parameters of the Standard Model, e.g. weak interaction parameters. 
It also gives predictions for the expansion of hadron masses around the zero quark mass 
limit, which can be used in the analysis of lattice calculations to extrapolate hadron 
masses from simulations at larger quark mass to the physical light {u, d) quark mass. 

The quark condensate occurs within the chiral perturbation theory formalism as a 
parameter which incorporates the dynamics of chiral symmetry breaking and can be 
determined e.g. using xPT and comparison to experiment or lattice QCD. 

Beside the chiral SU (2) invariance, the quark Dirac Lagrangian is symmetric under 
transformations of global SU (2) vector transformations. The left-and right-handed com- 
ponents of the field q therefore have an invariance under the independent transformations 



and the corresponding transformations for and qj^. If chiral symmetry is explicitly 
broken, the Lagrangian is only invariant under simultaneous (ez, = eji) transformations of 
left-and right-handed quarks. 

We begin with considering an effective Lagrangian for the pions coupling to u and d 
quarks which should be chirally symmetric up to quark mass terms which act as external 
sources for the pions. The pions are massless except for non- vanishing quark masses which 
provide the couplings to the chiral condensate. It is common to use a nonlinear a model 
to represent the pions: 



where the vr* are the pion fields, and the r* are Pauli matrices, /tt is the pion decay 
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constant. U transforms according to the (2,2) representation of SU{2)l x SU{2)ii: 

U{x) LU{x)R\ (88) 

where R G SU{2)r and L G SU{2)l. The expansion of the pion Lagrangian begins at 
0{p^), i.e. quadratic in the momenta and hnear in the quark mass: 

= Tr [d,U^d>^U + x^U + xU^)] ■ (89) 

X = 2Bq}A, where M. is the quark mass matrix. If u and d quarks have masses m„ and 
irid-, the mass matrix is: 

M^iZ- (90) 





Bq is proportional to the chiral condensate: Bq = — {qq) / f^. Keeping only quadratic terms 
in the pion fields one obtains a free Klein Gordon equation for canonically normalized fields 
TT* with mass rupg = BqTt A4/2, i.e. 

J(W)I. 



mis = ^^ra„ (91) 



the Gell-Mann-Oakes-Renner (GMOR) relation. The relation has been verified by com- 
paring with experimental results to hold to a high accuracy [214]. The result is taken 
as confirmation that the chiral condensate is the leading order parameter of spontaneous 
chiral symmetry breaking. In lattice calculations, is used as a measure of the quark 
mass. For heavier quark masses, the GMOR relation is modified as can be seen from 
lattice results [215,216]. 

At the next order, the pion Lagrangian contains many additional terms [125] with 
couplings which have to be determined from comparison with experiment or lattice QCD 
predictions. If the couplings to external fields are not taken into account, the O(p^) terms 
Lagrangian are [213]: 

4) = J/i [Tild^U^df'U] Y + J/2 TT[d^U^d,U] TT[d''U^d''U] 

+ J(/3 + k) {Tt[x^U + Ux^]Y + \kTT[d,U^d,U] Tiix^U + xU^] 



{Tiix^U - Ux^Y + {Tt[x^U - Ux^Y -2Tt[x^Ux^U + U^xU^] (92) 



Extension of the formalism to quenched QCD, i.e. without dynamical sea quarks, is 
discussed e.g. in [135]) and references therein, and extension to partially quenched QCD 
with different valence and sea quark masses is discussed in [133, 134]. 

If the pions are light compared to the inverse lattice extent, but the lattice size is 
not too small, one can study fluctuations around the zero modes within the so-called e 
expansion of xPT [136]. 



38 



4.1.1 xPT for relativistic baryons 

The nucleoli (and other hadrons such as the p) remain massive in the chiral hmit, and 
their mass cannot be directly used as a small expansion parameter. However typical 
momentum scales for the interaction with pions can still be of the order of the pion mass. 
Here, we employ xPT to predict self-energy corrections to the nucleon mass at finite pion 
mass. 

The Lagrangian description of low-energy nucleon-pion interactions in chiral pertur- 
bation theory at lowest order (0(p^)) is 

£(1) = ^ (^^^DM _ rno) ^ + i(?^^7/.75W^^, (93) 

where ^ is the Dirac spinor of the nucleon. 

qa is the nucleon axial charge, related to the nucleon-pion coupling G-^jq and the pion 
decay constant through the Goldberger-Treiman relation 

= (94) 

J IT 

Although the nucleon is massive, the Lagrangian is invariant under the chiral transfor- 
mations 



'^{x) V{x)'^{x), 



(95) 



The transformation matrix V is a product of L, R and the pion fields U{x): 

V{x) = R [R^LU{x)\ ^'^ [U{xf'^] ^ = L [L^RU\x)\ ^'^ Uixf'^. (96) 
with = U, the covariant derivative 

D, = d^ + ^[uld^u], (97) 

and the pion current = iu^d^jUu^ . rriQ the nucleon mass in the chiral limit. 
The 0{p^) Lagrangian is given by 

£(2) = ciTr(x+)^^ + -%Ti(u.u^) C^D^D"'^ + h.c.) 

+ |Tr(nX)^^- j^7^7l^M'M^' (98) 
with x+ = u^x^'^ + ux'^u. 
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In xPT, the quark mass dependence of the mass and the finite size effect are de- 
termined by the mass renormahzation of the nucleon. The mass dependence in infinite 
volume has been used to calculate the chiral extrapolation formula, and the finite size 
effect is given by the difference of the self energy in finite and infinite volume. The renor- 
malized mass is determined by the pole of the nucleon propagator calculated to a certain 
order in the p expansion with the external momentum set on-shell, 

m = ^{(/=m), (99) 

where q is the momentum of the external nucleon line, and p = 0{(^ — m). In practical 
calculations, one will substitute for m its expansion in powers of p to the necessary order 
in Eq. (jHS))- Extrapolation of lattice nucleon masses to the chiral limit and the calculation 
of the effect of the lattice size on nucleon masses in QCD with Nf = 2 is discussed in 
Section 15.51 

In quenched QCD, chiral perturbation theory is modified because quark-antiquark 
pairs cannot be dynamically generated. Additional terms which are non-analytic (loga- 
rithmic) in the quark mass arise in quenched chiral perturbation theory. The modifica- 
tion for the nucleon mass is given after Eq. ()157j) . Modifications also occur in partially 
quenched QCD, i.e. if the sea and valence quark masses are not equal. Partially quenched 
chiral perturbation theory is reviewed in [222]. 

4.1.2 Heavy baryon xPT 

Since the nucleon mass is large, of the order of the cutoff scale of chiral perturbation theory, 
it is non-trivial to retain the power counting of the effective theory in calculations beyond 
tree-level within the relativistic formalism [130]. In Heavy Baryon Chiral Perturbation 
Theory [128] (HBxPT), an infinitely heavy field is used for the baryon, and finite mass 
corrections can be included within an expansion in inverse powers of Mb, where Mb is 
the baryon mass. 

If the baryon mass is considered to be very large, the interactions introduce only small 
fiuctuations k of the baryon momentum, and one can write: 

P = Mbv + k, (100) 

Corrections can be included in a 1/Mb expansion. The baryon Dirac fields are projected 
onto fields B^{x) which are projected to constant velocity: 

B^{x) = exp{iMBV ■ x)P+{v)B{x). (101) 

V is the velocity of the heavy nucleon. In the rest frame, the temporal component is one 
and the spatial components are zero: v = (1,0,0,0). P+{v) = |(1 -|- /;) projects onto the 
forward moving part of the nucleon field: 

P+B, = B,; (102) 
^B, = B,. (103) 



40 



Accordingly P-{v) = |(1 — /') projects onto the backwards moving part of the nucleon 
field. In the baryon rest frame, the operators P project onto the upper two or lower two 
components of the Dirac spinor. 

At 0{p^) the nucleon Lagrangian is given by: 

C^^l = N4tvD + gAS-u]N,. (104) 

with the definitions for the axial pion current and the covariant derivative = + 
as before. The axial vector coupling of the relativistic Lagrangian is substituted by the 
coupling to the Pauli-Lubanski spin vector: 

= ■^l5(^fiuv''. (105) 

The tree- level Lagrangian at O(p^), 

£gj = N {cTrix-,) + . . .} (106) 

includes the low-energy constant, ci, which has to be determined from experiment. The 
leading pion mass dependence is quadratic. In principle, Ci can also be fixed from the 
lattice, however we will see later that low order xPT does not fit lattice data at the 
present large quark masses well. 

A similar formalism exists for vector mesons [131]. 



4.1.3 Lattice heavy baryon 



There are various ways of regularizing xPT-, 6.g. with dimensional regularization or a 
cutoff scheme. One can also construct a lattice version of xPT for heavy baryons [132] (a 
modification for SU{2) and the representation (jS7j) for the pion fields is given in [152]). 
The inverse of the free pion propagator is then: 
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(107) 



The lattice baryon propagator in the HBxPT is 

'1 



1-e 



-iaq4 ^ 



:io8) 



Lattice regularization for chiral perturbation theory for pions is discussed in [220] . Use 
of a lattice formalism to calculate pion finite volume effects is described in [221]. 

At finite lattice spacing, the quark Lagrangian has a different symmetry structure 
from continuum QCD. Unless Wilson quarks are 0{a) improved, chiral symmetry is also 
broken by terms 0{a) at finite lattice spacing, and these terms have to be included in 
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the chiral Lagrangian. The leading terms have been derived for the unquenched theory 
in [65,223-225]. For a review on chiral perturbation theory for Wilson-like quarks on the 
lattice see e.g. Ref. [226]. The chiral symmetry group of staggered fermions is also reduced 
due to O(a^) corrections at finite lattice spacing. A chiral Lagrangian for staggered 
fermions was constructed in [227]. 

4.2 Heavy Quark Effective Theories 

In effective theories for heavy quarks, one makes use of the fact that although the b 
quark mass for example is around 5 GeV, the mass splittings of heavy hadrons are only 
of 0(few 100 MeV), and that the dynamics of a very heavy quark in a hadron can be 
described as a small correction to the infinite quark mass limit. Using an effective theory 
where momentum scales > M are integrated out, a 1/M expansion for various quantities 
can be derived. 

In Fig. El on the left we show mass splittings of heavy-light hadrons averaged over 
spin-orientations. H stands for a generic heavy-light meson where the heavy quark is 
one of s, c, b, and Hs for a meson with a. c or b quark and a strange quark. Overbars 
denote averages over spin orientations. -^3/2 denotes the P wave meson doublet with 
light quark angular momentum 3/2 averaged over the heavy quark spin orientations. For 
an explanation of the P wave mesons and of the baryons see Eqs. ()210j) and ()211|) . They 




l/(M[GeV]) l/(M[GeV]) 



Figure 9: Mass splittings of heavy-light hadrons. On the left, splittings of states averaged 
over heavy quark spin orientations. On the right, hyperfine splittings between states of 
different heavy quark spin orientation. 

are largely independent of the heavy quark mass and have a finite limit in the limit of 
infinite mass. Experimental results on spin splittings of heavy-light hadrons are shown 
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in Fig. ini on the right. Circles show ground state S wave hyperfine sphttings, squares 
ji = 3/2 P wave hyperfine sphttings and crosses ground state S wave ji = 1 baryon 
hyperfine sphttings, where ji is the hght quark angular momentum in the quark model. 
For an explanation of the P wave mesons and of the baryon states see Eqs. ()210|) and 

(EHD- 

One observes that the spin splittings behave very roughly proportional to the inverse 
heavy meson mass and appear to vanish in the limit of infinite mass. 

4.2.1 HQET 

A quark with infinite mass Mq does not experience recoil from interactions with gluons 
and moves at constant velocity with momentum Mqv, i.e. it is static in the rest frame of 
the hadron. The spinor components of the heavy quark decouple, which gives rise to a 
spin symmetry, and the quark can be regarded as a simple color source. Energy splittings 
are determined by the interactions of gluons with momenta 0{Aqcd)- Since in B and D 
mesons the energy splittings are much smaller than the heavy quark mass, one expects to 
be able to describe them to a certain degree of accuracy within an effective theory where 
energy scales above the heavy quark mass are integrated out. 

At finite heavy quark mass, exchange of gluons with momentum 0{Aqcd) exchange 
leads to a recoil momentum k of the heavy quark with k ~ 0{Aqcd) which is small 
compared to the momentum of an on-shell heavy quark which is non-relativistically given 
by Mqv. Recoil and spin fiip of the heavy quark give rise to corrections 0{Aqcd)/Mq to 
the expectation values of the physical quantities, which are described by the interactions 
of higher- dimensional operators consisting of quark and gluon fields. The coefficients of 
the operators have to be determined by matching to full QCD. Since energy scales larger 
than the heavy quark mass are O (several GeV), the physics they correspond to, which is 
integrated out in the effective theory, is to a large extent perturbative. 

For a quark of mass Mq moving at velocity v, one can introduce a velocity projected 
heavy quark field which is related to the heavy quark Dirac field Q through the unitary 
transformation 

We calculate the tree-level Lagrangians for at zeroth and first order in 1/Mq according 
to the description of [2]. Undoing the projection, the equation of motion for Q^, is 

[lP+Mq{^-1)]Q, = 0. (110) 

Using the projection operators introduced in the HBxPT section, one finds the forward 
(+) or backwards (— ) moving fields 

= (111) 

Multiplying Eq. (jllOj) with P+ and using = 2a ■ b — one finds 

iv-DQt = -P^ipQ-. (112) 



43 



An analogous equation can be obtained by multiplying Eq. (jllOj) with P : 

iv-DQ- + 2MqQ~ = P~i]pQ+. (113) 

If one assigns 0{Mq) to and 0{1/Mq) to Q~ , one finds for the equation of motion at 
zeroth order in 1/Mq: 

iv-DQy = 0. (114) 
For the zeroth order Lagrangian in the 1/Mq expansion one then obtains 

£(0) = Q^ivDQ^. (115) 
From this one can read off the free propagator for the Q^, fields, 

G^^\k) = / . , (116) 

and the tree-level vertex 

igt^v^". (117) 

In the following we derive the 1/Mq corrections to the Lagrangian at tree level. Re- 
ordering Eq. ()113p . one finds 

<5,7 = -^P^iJpQt - i^^^Q7, (118) 
2Mq 2Mq^" ^ ' 

The second term is 0(1/Mq) relative to the first term. Substituting this into Eq. pi()|l 
one obtains 

%v ■ DQt = -P+iIp^p-iIpQt. (119) 
This can be simplified further. Using 

-^P-^If^ = —L[lf^-2{D. vf] (120) 

and G^y = ^[-D^, D^], one finds 

]/^=^-ga^''G,. + D\ (121) 



where a^'^ = |[7^, 7^^]. The 1/M Lagrangian is therefore given by 

Qv. (122) 



2 

1 



D'-{v DY + -ga^'G,, 



In the B meson rest frame the spatial components of the velocity are zero, i.e. — {v 
DY = —D^. Then the two contributions to Eq. p22|l are recognizable as the 
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kinetic energy operator 



IMn 



• spin-magnetic energy operator 

Omag = -T7-a5^-5Q,, (124) 

Mq 

with the spin operator = |757°7'^ and the magnetic field B'^ = —^e'^^^F^''. In the 
Heavy Quark Effective Theory (HQET), the mass of a heavy-hght hadron H therefore 
consists of the contributions 



Mh = Mq + A ^ 



2M, 



Q 



{H\QD^Q\H) {H\Qa-BQ\H) 
2Mh ^ 2Mh 



+ 0(1/M^), (125) 



where A the binding energy of the meson for Mq — > oo, and the other two terms are 
the expectation values of the heavy quark kinetic energy and the spin-colormagnetic 
interaction energy respectively. 

The HQET predictions for the mass dependence of decay form factors in the heavy 
quark effective theory are useful for the analysis of experimental and lattice results of 
form factors, e.g. in fits of the mass dependence. We summarize the discussion of the 
infinite mass limit of Jb [2]. In QCD, the decay constant of a meson H with a heavy and 
a light quark can be defined through {0\A^{0)\H{p)) = —ip^fn (Eq- (EH))- We further 
assume that the meson states are normalized relativistically, according to 

{H{p')\H{p)) = 2/(27r)=^5(=^)(p - (126) 

The decay matrix element at infinite mass in the HQET should be similarly related to 
the heavy quark velocity 

{0\A^{0)\Hiv)) = -tfHV,, (127) 

where we have denoted the HQET quantities with ~ . Using a non-relativistic normaliza- 
tion for the HQET fields: 

{H{v')\H{v)) = 2v\27Tf6^''\p-^) (128) 
one obtains the heavy quark limit of fn'- 

fH = fHVM^, (129) 



i.e. fEvMn is Mg-independent up to l/Mg corrections, and the tree level mass depen- 
dence of /hV^h is expected to go as 

fn^/M^ = ^+^ + ^ + (0(Mg3)). (130) 



45 



In lattice calculations, when matrix elements are studied at a range of heavy quark masses, 
Eq. (jlHOj) or a renormalized version (see Eq. is used to fit the results as a function 

of heavy quark or meson masses. 

Beyond tree level, renormalization effects have to be taken into account. We attempt 
to describe some general features, following the discussion in [208]. Since we will discuss 
lattice results on Jb, we are interested in the matching of the axial vector current. In QCD 
with relativistic quarks, the axial current is protected by chiral symmetry and is therefore 
not renormalized. To calculate the QCD matrix element from the matrix elements in 
the effective theory, one requires that the matrix elements in full QCD and the effective 
theory match at an energy scale /x. This statement can be expressed by: 

A?^^ = (131) 

where C is a Wilson coefficient, at leading order in HQET. In general there can be several 
operators contributing to the right hand side which mix under renormalization, and we 
will discuss this in the context of matching the axial vector current in lattice NRQCD to 
continuum QCD. 

A typical momentum scale for the gluons exchanged in the interaction is yU. = Aqcd- 
The scale dependence of the coefficient C is determined by a renormalization group equa- 
tion: 

d 



^dil ~ ^^^^ ^ °' ^^^^^ 

where 7^ is the anomalous dimension of the axial vector current. The cutoff for the 
effective theory being the heavy quark mass implies the initial condition 

C{^l = Mq) = 1. (133) 

The Wilson coefficients C contain logarithms of fi/rriQ which are large for typical choices 
of /i. At leading log approximation the result for the scaling behavior of the coefficients 
is 

The large logarithms can be resummed if the coefficients are multiplied by the inverse of 
the factor p34j) . The axial current matrix elements /i-y/Mi and f2\/^2 of two heavy-light 
mesons with different heavy quarks of masses mi and m2 are then related in HQET by 



(135) 



/2V^ \asiM2) 

i.e. Eq. fll3()j) is modified by logarithmic corrections. 

Beyond tree-level, one expects a similar renormalization of the operators in the HQET 
Lagrangian. For the spin-magnetic operator, 

Omag -cb{ij)^QvS ■ BQy. (136) 

Mn 
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Scale invariance of physical mass splittings gives rise to a similar renormalization group 
equation as for the decay constant. Using 



CB{fi = MQ) = l, (137) 
the leading-log scaling behavior of the coefficient is [228] 



The kinetic energy operator is not renormalized due to invariance under reparameteriza- 
tions [229] 

e 

V V + 



Mr 



Q. - ^^^-Ml + ^l- (139) 



4.2.2 NRQCD 

Non-relativistic QCD (NRQCD) is an effective theory formulated for heavy quarks as- 
suming that their dynamics is non-relativistic, with correction terms which can be added 
within a systematic expansion. For quarkonia the higher order interactions are arranged 
in a f ^ expansion, where v is the heavy quark velocity (see e.g. [4]). NRQCD has been very 
successful in studying quarkonia. For a review discussing effective field theory methods 
in studies of quarkonia see [230]. One may say that lattice calculations of the T use to a 
large extent non-relativistic effective field theory, and it has been applied with very good 
success (see e.g. the results in [87,175,193,195]). We now discuss NRQCD for heavy-light 
mesons. 

In hadrons with heavy and light quarks it is an expansion in v or 1/Mq, where Mq 
is the heavy quark mass, similar to HQET. At infinite mass, the heavy quark is just a 
source of the color electromagnetic field, whereas at finite Mq, there is a recoil of the 
heavy quark due to the interaction with soft gluons with typical momenta of 0{Aqcd)- 
One can argue that the heavy quark momentum Pq and light quark momentum pq are 
equal due to momentum conservation in the rest frame of the meson: 

Mqv^Pq = Pgr^OiAQcn). (140) 

Therefore v = p/Mq ~ Aqcd/Mq ~ 0.1 in B mesons and should be a useful expansion 
parameter to specify corrections to the Mq oo (static) limit for quantities which involve 
a small heavy quark momentum, such as the spectrum or leptonic decay constants. 

Contributions at 0{1/Mq) are the kinetic and the spin-colormagnetic energy of the 
heavy quark; at 0{1/Mq) a heavy quark spin-orbit interaction and a Darwin term are 
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included. The tree-level terms in the 1 / Mq series can for example be obtained with Foldy- 
Wouthuysen transformations on the heavy Dirac spinor which disentangle the interactions 
between the upper and lower two components of the Dirac spinor (see e.g. the discussion 
in [231]). The Dirac Hamiltonian contains spatial 7 matrices which are off-diagonal. Using 
unitary transformations which act on the heavy quark Dirac fields Q as 

Chosing S to be proportional to 7 ■ D/Mq, one obtains decoupling of the Pauli equation 
for the heavy quark and a corresponding equation for the antiquark up to higher orders 
in I/Mq. Mq is the bare heavy quark mass and agrees with the previous Mq at tree level. 
The 0(1/M^) Lagrangian for the heavy quark is finally given by 

C = ^|J\Dt + H)^^J, (142) 

with the heavy quark Pauli spinor ijj and the Hamiltonian 

2Mo Mo 

, 'Igjn.E-E.n)-^^,.0.E-E.O)-^-^. (143) 

E is the color-electric field with E^'^ = —P°-^^'^^ where a is a color and i a space index. The 
last term is the first relativistic correction to the kinetic energy of the heavy quark, which 
is usually included in calculations at 0(1/Mq). The coefficients of the various terms can 
be found with matching calculations to full QCD in the continuum. 



5 THE LIGHT HADRON SPECTRUM 

Light hadrons, e.g. the vr, p, p and n and their excited states, are of course among the 
most interesting states to be studied on the lattice. To understand the dynamics of the 
forces between quarks and to study fundamental properties of QCD such as confinement 
nonperturbatively, one investigates the potential between heavy quarks. For a detailed 
discussion see e.g. [92]. Here the potential is discussed in the context of setting the lattice 
scale and to calculate the strong coupling constant used in perturbative calculations. Only 
a very brief introduction is given. 

In this section we are interested in hadrons consisting of u, d and s quarks in the quark 
model. Using quenched gauge fields, Wilson quarks and quark-model inspired lattice 
operators, it was found already some time ago that the features of the experimental light 
hadron spectrum are reproduced well by the lattice [138,139]. Progress in recent years 
was made using improved lattice actions and/or unquenched gauge fields. In this section 
we discuss in particular results on light meson masses, meson decay constants and the 
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light baryon spectrum from [111] using improved quark and gluon actions. We give a brief 
outline of the method used in [111] as an example for a standard lattice calculation, and 
then discuss some aspects which could not have been discussed in detail in the paper [111], 
re-analyse the comparison of quenched and dynamical results, and compare to other new 
lattice calculations. The chiral extrapolation and the finite size effect of the nucleon is 
discussed using numerical results from QCDSF, UKQCD, CP-PACS and JLQCD. 

The lattice simulations discussed here have been performed around the strange mass. 
To avoid extrapolation uncertainties we compare not extrapolated, i.e. strange lattice 
hadron masses, with experiment, and keep in mind that with using strange sea quark 
masses a systematic uncertainty may be included. 



5.1 Technical Details 

We now discuss some technical details of the calculation described in Ref. [Ill], as an 
example for a standard lattice simulation with two flavors. The explanation is kept close 
to the paper [111]. 

The operators employed to project onto the mesons and bayons are constructed of 
quark fields such that they project onto the quantum numbers of the state under consid- 
eration. It is usual to construct them in the same way as quark model wave functions to 
ensure a good overlap with the hadron state. 

The baryons under consideration in this study are the ground state spin 1/2 and spin 
3/2 baryons. The operators used for the spin 1/2 octet baryons are 

0„ = e'^'^ [qfC^ml] qla, (144) 

where C = 7472 is the charge conjugation matrix. The subscripts 1, 2, 3 denote the 
different quark flavors. The index a denotes the spinor component of the quarks which 
are not summed over and spinor components of the baryon fields, a, b, c are color indices. 
The quark model wave function is antisymmetric under interchange of color indices. The 
spin and flavor wave functions are symmetric under interchange of two quarks, but not 
totally symmetric. The operators for the octet baryons except for the A are given by 

-l=e^^^{[qtq'M+[q^,ql]ql), (145) 

with [^1^2] = qiq2 — qiqi, where the operator C75 is not written out explicitly. For the A, 
the operator 

is used. The spinor components of the baryon are given by the spinor components of the 
third quark. For the spin 3/2 decuplet baryons, operators projecting onto specific values 
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of the three-component of the baryon spin depend on the spinor component of the diquark 
and of the third quark: 

O3/2 = e'^'^{qfC^+4)ql, 

0-1/2 = e^'''l[{q''^c^3q2)ql2-[{q^^c^-q2)qli] 

0-3/2 = e^'^{q^^C^.4)ql„ (147) 

where the first subscript denote the flavor and the second the Dirac spinor components. 
7 = (71 T ^72)/2. The gauge fields used in this simulation are renormahzation group- 
improved [56]. Both sea and valence quarks are clover quarks with the clover coefficient 
csw being tree-level mean-field improved using the plaquette calculated at one loop per- 
turbation theory [140]. This choice agrees rather well with the one-loop result, as shown 
in Table El 





(tad) 




1.8 


1.60 


1.69 


1.95 


1.53 


1.60 


2.1 


1.47 


1.53 


2.2 


1.44 


1.49 



Table 6: Tadpole-improved and one-loop clover coefficient csw for renormahzation group 
improved glue with Nf = 2. 

The hadron masses are extracted using fits of the hadron correlation functions pro- 
jected onto zero spatial momenta, as a function of Euclidean time (the idea is outlined in 
Sec. 13. 3|) . The masses are obtained for the values of the simulated bare gauge couplings 
and quark masses and have to be interpolated or extrapolated to the physical quark 
masses, defined by the condition that certain hadron masses or mass ratios take physical 
values. 

As mentioned before the pion mass squared is a rough measure for the quark mass, 
and the meson masses from [111] can indeed be described as functions of the pion mass 
squared. It is found that vector meson masses can be fitted for a given sea and valence 
quark mass combination (^sca, i^vi, i^v2) as functions of the degenerate valence pion mass 

/iv = ^ (ah + f^h) (148) 
and the pion mass for degenerate valence and sea quark masses 

/is = fTT'PS (^sea; '^sca; ^^sea) ; (149) 
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including up to quadratic terms in the quark masses: 

m = A + Bsfis + -Bv/iv + Cs/ig + Cv/^v + C'sv/isA'-v (150) 



|3=2.1,K=0.1382 p=2.1, K =0.1357 




1/(2k) 1/(2k^) 

Figure 10: Squared pion mass as a function of the averaged valence hopping parameter k^, 
for the hghtest = 0.1382) and the heaviest = 0.1357) sea quark at /3 = 2.1. Squares 
denote the actual pion for a given k combination, and diamonds the average squared mass 
from degenerate valence quarks. Dashed lines show linear fits to the averaged, solid lines 
to the non-degenerate data. 

The decuplet baryons are also fitted to the functional form given in Eq. p50j) . where 
/iv corresponds to the valence quark mass averaged over all three quarks. 

The octet baryon masses are described as functions of fi^i and yUv2 (corresponding 
to the masses of both valence quark flavors) and the sea quark mass, including up to 
quadratic terms in /ij. One should note that the mass of the actual pion with valence 
quark parameters k^i and agrees with the averaged mass of two degenerate pions 
with the respective mass parameters k^i and only if the quark masses are similar. 
In Fig. ^1 we show the averaged and real pion masses in lattice units as functions of the 
VWl quark mass shifted by the critical mass. 

The lattice spacing and the quark masses can now be fixed to the physical point. For 
the calculation of quark masses and spectrum in [111], the lattice spacing is determined 
from rrip, and the determination of the lattice spacing and the determination of the k values 
corresponding to the physical u, d and s quark masses are related. The procedure used in 
Ref. [Ill] is described in the following. First the pion and rho masses in lattice units are 
extrapolated to the physical ratio m^/mp using a diagonal extrapolation (setting fiyai = 
fisea in Eq. ()150j) . This gives the pion mass in lattice units, which is in turn substituted in 
Eq. p5()j) to determine arup, from which one obtains the lattice spacing a. ni, the k, value 
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corresponding to the light {u, d) quark mass is calculated using a polynomial fit of m^ps 
in terms of the VWI quark masses Eq. and setting m1>g{ni; k;, ni) = m^. The k value 
corresponding to the strange quark mass (/t^) is determined using strange meson masses. 
One method is to use the K meson and to fix the ratio 7ii?pg{Ki \ ki, Ks)/^^ to the physical 
value m\;/m\. Then one can define a fictitious strange-strange pseudoscalar meson, the 
rjssi from the relation m^^^ = m?pg{ni\ Hs, i^s)- '^'^iVss) can then be used to parameterize 
the strange quark physical point in the fits of the other hadrons in terms of the /I's. The 
strange quark mass is alternatively determined using the 0, an ss vector meson. 

Parameterizing the using Eq. ()150|) with ^g^a set to and fixing the ratio m^/mp 
to the physical value one obtains a result for n^ai = Vss which is used to fix the strange 
point for baryon spectrum. Using the (f) to determine the strange quark mass, the K 
meson mass is a prediction, determined by the equation m\ = rn?pg{Ki; ki, k^). 

The continuum results are obtained using linear extrapolations in a. 

5.2 Light Meson Masses 
5.2.1 Quenched results 

The quenched approximation is only an approximation, and some physics is expected to 
be different in the quenched world from real QCD. 

Since the p mass is an important result of QCD and CP-PACS's favorite quantity to 
set the lattice scale, we first compare the quenched vector meson mass to experiment, 
using To to set the scale. Plotting the chiral extrapolation of the p mass from the coarsest 
and finest lattice of [111] in Fig. ^2 one sees that the result does not agree with ex- 
periment. To investigate whether this may be due to discretization errors and whether 
the behavior changes qualitatively if yet smaller quark masses are used, we compare to 
a set of recent quenched results using improved gauge field and fermion actions, which 
are expected to scale well, at small quark masses. Refs. [121,141,142] use chiral lattice 
fermions. We compare the results using chirally improved fermions of [120], and results 
using clover fermions with a fattened-link Wilson term which also enables simulations at 
small quark masses [143]. We find an agreement of the CP-PACS lines with most of the 
new lattice results at larger quark masses. For smaller quark masses masses both overlap 
results seem to extrapolate to higher values of rrip. A possible reason for the difference 
of the results from overlap fermions is that [142] uses unimproved Wilson glue on coarse 
lattices and [141] the Liischer-Weisz gauge action [54] on fine lattices. The plot however 
indicates that the overlap quark results might extrapolate to a relatively larger value. The 
fixed point and chirally improved fermions [120, 121] appear to extrapolate towards the 
experimental result. Using a chiral extrapolation guided by the functional form predicted 
by quenched chiral perturbation theory they find quenched p and masses in agreement 
with experiment. 
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Figure 11: Quenched vector meson masses as a function of ?Tip^ in units of tq. The 
hues denote hnear fits to results from [111] at the coarsest (dashed) and finest (sohd) 
lattice spacings. Triangles up use overlap quarks from [141], diamonds fixed point quarks 
from [121], circles overlap quarks with coarse Wilson gauge fields from [142], squares use 
chirally improved quarks from [120] and triangles down use FLIC quarks from [143]. The 
experimental value is denoted by a star, with the error bar depicting the variation if 
tq = 0.47 fm instead of 0.5 fm is used. 

5.2.2 Unquenched results 

In FigurefT^lthe vector meson mass is shown as a function of the pion mass at finite lattice 
spacing, for equal sea and valence quark masses. The agreement of data from various jS 
values is rather good, and a continuum extrapolation seems to be not necessary. The 
result is rather linear in the quark mass. A potential systematic error in this relation is 
due to the fact that the unquenched p can undergo hadronic decays. For the given set 
of masses, the decay thresholds of 2^vn?pg + {2t\ jVf' are about 400 MeV larger than the 
vector mass. 2mps is also higher but for small quark mass close to the vector masses. 

We have seen that the p mass agrees with the experimental value if the lattice scale 
is set using tq in unquenched QCD. It is of interest to study also extrapolations in the 
valence quark mass at finite fixed sea quark masses. In Fig. we compare the chiral 
extrapolation in the valence quark on a relatively fine lattice for two different values of 
the sea quark mass to the quenched extrapolation. Again, the scale is set using tq. As 
in the previous plots, we can use the tq values determined at the respective sea quark 
masses, and find agreement with experiment already at finite sea quark mass. In a mass- 
independent renormalization scheme, the lattice spacing depends only on the coupling, 
not on the quark mass [144] up to discretization errors, and we compare the meson masses 
also with chirally extrapolated tq values, Those were determined in [111] assuming 
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Figure 12: The p meson mass from [111]. Filled symbols denote the simulation results. 
Errors in the y direction are roughly of the size of the symbols. Open points indicate 
the corresponding thresholds given by 2ro Y/?7Tj^^r+~(27r7Zp. Diamonds stand for /3 = 
1.95, circles for (3 = 2.1 and squares for (3 = 2.2. The dashed curve indicates romy = 
2ro'mps- The solid line indicates a linear fit to the data at /5 = 2.20. The star denotes 
the experimental value. The error bar gives the variation if tq = 0.47 fm instead of 0.5 
fm is used. 

a linear quark mass dependence: 

— (/^sea) = i + C'/is. (151) 

We find that partially quenched results using the chirally extrapolated rg into the com- 
parison in Fig. ^1 also move towards the experimental value, while the agreement is less 
obvious at finite sea quark mass. 

5.3 Light Meson Decay Constants 

Lattice results for light decay constants are important tests of QCD as a theory of light 
hadron physics and also a calibration of the lattice method for calculations of hadron 
matrix elements. Since they are likely to be less sensitive to systematic errors, it is of 
interest to calculate ratios of light and heavy-light decay constants. We discuss the results 
of [111]. 

The decay constant of a pseudoscalar meson P of mass M in Euclidean space-time 
can be defined as: 

{0\A,{0)\P{p = 0)) = fpM. (152) 
For the chiral extrapolation of the matrix elements, the same fitting form as for the 
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Figure 13: Partially quenched vector meson mass as a function of {mpsY at /5 = 2.1. 
Diamonds denote K^ca. = 0.1367, squares Kgea = 0.1382. Filled symbols correspond to the 
ro values with equal sea and valence quark masses. Open symbols correspond to chirally 
extrapolated values of tq. Error bars are smaller than the symbols. The star denotes 
the experimental value. Fat dot-dashed and dashed lines denote the quenched fits on the 
finest and coarsest lattice respectively. Thin dashed lines denote quadratic fits to the 
various data sets. 

vector meson masses is used, Eq. ()15()|1 . The quenched decay constants can be described 
by a linear function in fi^ only. 

In [111] the decay constant was calculated with one-loop tadpole-improved renormal- 
ization constants. The axial vector current could then be expressed by the relation: 



with one-loop values for Z^, and ca- The effect of the improvement term proportional 
to Cyi on the decay constant is a few MeV. 

5.3.1 Quenched results 

The quenched results for and fx appear to have reasonably small discretization errors 
at lattice spacings up to around a ~ 0.15 fm. The plateau occurs however at a roughly 
20% higher value than the quenched result of [139] using Wilson quarks. 

Using the result from the finest lattice, one can quote for the quenched light meson 
decay constants 




(x) = 2hiUoZA (1 + hAm/uo) {A^{x) - CAd^P{x)) , 



(153) 



fK 



145(5) MeV, 

162(4) MeV, K input, 

166(4) MeV, (p input. 



(154) 
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i.e. 



fx 



1.12(5) , K input 



1.15(4) , input 



(155) 



Errors are statistical and extrapolation errors only, fx is close to experiment, while is 
~ 10% higher than experiment. 

5.3.2 Unquenched results 

Scaling violations in the unquenched decay constants are large, i.e. 50% between 0.2 > 
a > 0.1 fm. To avoid systematic errors from the continuum extrapolation, we quote 
results from the second finest unquenched lattice at /5 = 2.1, since at /? = 2.2 the error 
bars are very large. 



where only statistical errors are quoted. The results are to be compared to the experi- 
mental values of 132 MeV for and 160 MeV for fx- 

In the ratio of vector and pseudoscalar meson decay constants, fp/ f-K-, discretization 
errors largely cancel, and the continuum extrapolation of the ratio is in agreement with 
experiment, however within large error bars. 

5.4 Light Baryons 
5.4.1 Quenched results 

First we would like to compare the quenched baryon masses from [111] to other recent 
quenched calculations and experiment. The values were extracted using the method 
described in Sec. 15.11 In the quenched calculation, the decuplet baryon masses were 
extrapolated using just a linear fit in the pion mass, while for the octet states a quadratic 
term was included. The results were continuum extrapolated using a linear fit in a. The 
quenched spectrum from this procedure is found to agree with the quenched continuum 
extrapolated results using Wilson gauge and fermion actions, see Fig. El The lattice 
spacing has been set with in both cases. The features of the experimental spectrum 
are reproduced, while strange baryons with K input and with input may differ among 
each other and with experiment by several a. Before trying to make definite statements 
of whether this is a quenching effect, we should consider possible other systematic errors. 
It is likely that discretization artifacts will affect also any unquenched calculation with 
the same action at similar quark masses and lattice spacings, and we therefore try to give 
some more thoughts to the question. 
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Figure 14: Quenched light baryon spectrum from CP-PACS from lattices with a ~ 0.1 
fm [111]. 



In the present case, the quark action is tree-level tadpole improved. The dominant 
discretization effects should be O(a^) and 0{asa), i.e. in principle there are linear and 
quadratic contributions in a. In Fig. UHl we show the nucleon mass as a function of the 
lattice spacing. 

We find that the data from the finer lattices scale well, i.e. they are independent of the 
lattice spacing, and a linear continuum extrapolation including all data points does not 
appear more reliable than using the result from the finest lattice as continuum estimate. 
We plot the quenched light baryon spectrum using the results from the finest quenched 
lattice at /3 = 2.575 in Fig. El The result changes by several a, and in most cases, 
agreement with experiment is better for the results which are not continuum extrapolated. 

Chiral extrapolation uncertainty could be another source of error. The extrapolated 
nucleon masses from CP-PACS using naive polynomial extrapolations agree quite well 
with experiment. The results at finite quark masses are not available, therefore we show 
only the chirally extrapolated values from the finest and the coarsest quenched lattice in 
Fig. El The fixed point [121] The overlap results from QCDSF [141] and the chirally 
improved quarks [120] seem to extrapolate to the experimental value. For the fixed point 
results plotted here [121] it is not so clear whether the experimental value is reached, 
however their results from coarser lattices and smaller quark masses do extrapolate to the 
experimental value [120]. For the overlap results from [142] and the FLIC fermions we 
have used here [143], the extrapolation is again not obvious from the data on this plot. 
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Figure 15: Quenched nucleon masses from CP-PACS [111] as function of the lattice 
spacing (squares) with the continuum extrapolated value (diamond) on the left. The line 
denotes the experimental value for the nucleon mass. 

5.4.2 Unquenched results 

Ref. [85] works at a fixed lattice spacing and thus cannot perform a continuum extrapo- 
lation. For the chiral extrapolation, a similar formula as in [111] is used. In the baryon 
mass calculation of Ref. [113] with two light and one strange sea quarks, the continuum 
extrapolation is performed first, and then a chiral extrapolation is done. 

In Fig. El "we compare the the baryon spectrum from the recent unquenched simula- 
tions of [85,111,113]. Ref. [85] assigns additional systematic errors of up to 25 MeV from 
the chiral extrapolation uncertainty and the determination of r^. We do not add them to 
Table [7| and just note that their additional error bars would point towards the opposite 
direction to the experimental value. We find discrepancies with experiment of up to ~ 2a 
for the nucleon and A, where improvement with respect to the quenched spectrum is not 
visible. Chiral extrapolation uncertainty might be relevant. An attempt to use chiral 
perturbation theory to extrapolate the nucleon mass is discussed in the next section. In 
Table [7| we give results for light baryon mass splittings corresponding to the results shown 
in Fig. El For the splittings, the agreement with experiment is at the 1 — 2a level. A 
recent calculation [113] with Nf = 2iight + Istrange finds a A — splitting which agrees 
well with experiment. The mass ratio tjin /mp has also been calculated using domain-wall 
fermions with Nf = and Nf = 2 [123]. Their quenched results seem to extrapolate to a 
value in agreement with experiment. Their unquenched value is in quite good agreement 
with experiment. 

These lattice studies do not take into account that in the unquenched theory the A can 
decay strongly; it is an A^vr resonance with a Breit-Wigner width of 115 — 125 MeV [146] 
and one would expect a shift in the mass. 
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2 4 6 

Figure 16: Quenched nucleon mass as a function of in units of tq. The chirally ex- 
trapolated results from [111] from the coarsest and finest quenched lattices are denoted by 
crosses (slightly shifted to the right for clarity), triangles use overlap fermions from [141], 
diamonds fixed point fermions from [121], circles overlap fermions with coarse Wilson 
gauge fields from [142], squares chirally improved fermions from [120] and triangles down 
FLIC fermions from [143]. The star denotes the experimental value, with the error bar 
giving the variation if tq = 0.47 fm instead of 0.5 fm is used. 

To avoid the uncertainties due to extrapolation in the comparison beween theory and 
experiment, we now consider relations between strange hadron masses in two-flavor QCD. 

In Figure ^1 we plot the decuplet baryon mass from three different /5 values versus 
the vector meson mass in units of tq. Sea and valence quark masses are degenerate. The 
physical point corresponds to the position of the meson and the Q~ baryon mass and 
agrees with the lattice result. Mixing between the and the uj has been ignored at this 
point. 

In Fig.Hniwe show the decuplet /octet mass ratio as a function of the pion/octet baryon 
mass ratio for three (3 values from CP-PACS. The experimental points correspond to quark 
masses averaged over all three valence quarks. Sea and valence quark masses of the lattice 
data in the plot are taken to be equal. The 'experimental' ssu pseudoscalar meson mass is 
determined by interpolation using the physical values m'j^ and m^. Although the plot has 
to be interpreted with some caution since the sea quark mass is not extrapolated to the 
physical point, the lattice octet baryon masses cannot be well described as a function of 
the average quark mass, and the baryons can be affected by strong decays, we nevertheless 
find the agreement with experiment remarkable. 
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Figure 17: Unquenched light baryon spectrum from recent unquenched lattice calcula- 
tions. 

5.5 Nucleon Mass: Finite Size Effects and Chiral Extrapolation 

In this section we discuss chiral extrapolation and calculation of the finite size effect of 
the nucleon mass using chiral perturbation theory. 

We begin with the Heavy Baryon xPT. The self energy of the nucleon propagator 
is considered on-shell, i.e. with external momentum zero. The 0{p^) correction to the 
quadratic dependence on mps from HBxPT is represented by the loop diagram in Fig.l^ 
on the left. Evaluating the diagram one finds for the mass at 0{p^): 



In quenched chiral perturbation theory, the m^pg contributions are modified, and new 
terms proportional to mps and In mp^ are introduced [150]. 

One can also attempt to calculate the finite size effect using lattice xPT . The lattice 
provides a way to regularize integrals, if one use substitutes the continuum HBxPT 
actions for the nucleons and pions with lattice ones. In Fig. 15.51 we study the dependence 
of the finite volume effect on the lattice cutoff for various pion masses. For rups around 
500 MeV, one finds a 10% difference between a lattice cutoff of = 2 GeV and the 
continuum limit a ^ 0, whereas the difference is almost 50% if a lattice with a^^ = 1 
GeV is assumed. For rrips = 200 MeV, the differences are slightly less than 10% and 20% 
respectively. 

HBxPT at 0{p^) yields a pion mass dependence rather different from the nucleon 
mass lattice data [151,152]. Moreover, finite size effects predicted by HBxPT are found 
to be clearly too small with respect to the finite size effects of the lattice data. The 
contribution of high pion momenta in the internal loops is of particular interest since 
on a lattice of 1.6 fm extent the lowest mode satisfying periodic boundary conditions 
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Table 7: Light baryon mass splittings. The first error is statistical, the second is the 
difference from fixing the strange quark mass using the K or (p meson where applicable. 
The quantity used to fix the lattice scale is shown on the left. 



has a momentum of ~ 0.8 GeV, which implies that there is only one mode within the 
domain where x^T is theoretically expected to be valid. Using lattice regularized chiral 
perturbation theory one sees an indication that the contributions from pion modes at 
high momenta can indeed be sizable, in particular for larger pion masses. An example is 
plotted in Fig. 15.51 Including a A-nucleon coupling, HBxPT predicts roughly 50 — 60% 
of the actual finite size effect for pion masses around 500 MeV [152]. 

Also, since for the unquenched clover lattice results, quark masses are rather large 
such that mps/niN ~ 0.5, the non-relativistic approximation does not appear particularly 
justified. The calculation was therefore repeated with a relativistic formalism. 

In the relativistic formalism at 0{p^), the nucleon mass can be written as [153]: 

m = nriQ — Acini\,g + + einipg. (158) 

nia is the one-loop contribution to the mass renormalization. The corresponding Feynman 
diagram is shown in Fig. |^on the left. The term proportional to m'j^g comes from the 
tree-level Lagrangian C^'^\ and the term proportional to ei is necessary as a counterterm 
to regularize the loop integral, although it contributes to the Lagrangian formally only at 
0{p^). Using the infrared regularization scheme for baryons explained in [154], 

.3g\momls f d^P r 2 22 2 n n , -1-2 

= ^2 J [P ~ ~ - ^) + ^ej (159) 

in Minkowski space. To summarize the idea of infrared regularization one might say 
that the integral is decomposed into an infrared finite and an infrared singular part; the 
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Figure 18: Unquenched decuplet baryon mass as a function of the vector meson mass 
corresponding to the same gauge couphng and quark mass. Different symbols correspond 
to different values of the inverse bare gauge coupling. Valence and sea quark masses 
are equal. The experimental point corresponds to the meson and the Q baryon, with 
ro = 0.5 fm. 



infrared singular part is retained in the calculation which is relevant to the calculation 
of the low energy structure of the theory. The Feynman parameter then runs from to 
oo instead from to 1. Calculation of the integral yields the following expression for the 
nucleon mass at 0{p^): 



rriN 



mo — icinipg + 



ei(A) + 



3^1 



m 



PS 



4mQ 



IT 

— h arctan ■ 



2 In 



mps \ 
A J 



m 



PS 



m 



PS 



2 ^Aim?pgml - m%g 



(160) 



At 0{p^) there is an additional tadpole diagram with an insertion from the Ci term in 
C^'^^ (see Eq. EH and in Fig. EH middle), and the nucleon mass then reads: 



rriN 



mo — Acimj:,g + + rrih + eirripg. 



(161) 



The tadpole contribution, m^, is given by the integral 



rrib 



3 f d'^k 2cimls - (C2 + C3)(A;°)2 + c^k'' 



fl J (2vr) 



mpg — k"^ 



le 



(162) 



For this diagram, infrared and dimensional regularization give the same result. The 
diagram on the left of Fig. |^ also receives a correction due to a third diagram with a ci 
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Figure 19: Unquenched mass ratio of decuplet and octet baryons. Different symbols 
correspond to different values of the inverse bare gauge coupling. Sea and valence quark 
masses are equal. The experimental points respectively correspond to the physical light 
quark mass and quark mass averaged over the ssl quark triplet, with ro = 0.5 fm. 



insertion to the nucleon propagator (the right diagram in Fig. 1211) . If the nucleon mass is 
redefined using rriQ — ^ rriQ — Acim'j:,^, the contribution of the additional graph is absent. 
The arctan can be expanded in powers of mps, with the first terms being: 



rriN 
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+ 0{m 



6 1 

PS) 



(163) 



If chiral perturbation theory is a valid low-energy description of QCD, this formula should 
be predict the relation between the pion and nucleon mass if the quark mass is small. 

A collection of lattice nucleon masses for Nf = 2 from clover-improved Wilson quarks is 
plotted in Fig. The QCDSF and UKQCD collaboration [148] uses non-perturbatively 
0(a) improved clover quarks and Wilson gluons. The UKQCD analysis is discussed 
in [149]. The JLQCD collaboration uses the same quark and gluon actions [122], and 
the CP-PACS collaboration uses tadpole- improved clover quarks and RG-improved glu- 
ons [111]. 

The lattice spacing is set using tq = 0.5 fm at each finite value of the sea quark mass. 
The data are not on a universal curve. If we compare results at various values of /5, it 
appears that discretization effects in the relation between pion and nucleon masses are 
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Figure 20: Finite size effects on the nucleon mass from lattice heavy baryon chiral per- 
turbation theory at 0{p^) as a function of the inverse lattice spacing. 




Figure 21: Nucleon self-energy diagrams in xPT at O(p^). 

small. In Fig. 1221 we plot results from QCDSF and JLQCD as a function of the spatial 
lattice size at fixed values of simulation parameters /? and rrisea- For smaller pion masses, 
one finds an increase of the mass on smaller lattices. If only data on larger lattices are 
considered (here mpsL > 5), they do lie on a universal curve and we use the results with 
pion masses < 800 MeV to fit to the chiral perturbation theory formula. 

Setting C2 = 3.2 GeV~^ which is in reasonable agreement with the experimental values 
and C3 = —3.4 GeV~^, the other parameters can be determined by a fit to the lattice results 
from large lattices and are 

mo = 0.89(6) GeV, (164) 
ci = -0.93(5) GeV-\ (165) 
ei = 2.8(4) GeV-^ (166) 

(167) 

with gA = 1.267, A = 1 GeV. 

As seen in Fig |^ the formula with these parameter choices is in good agreement 
with the data, and the chiral extrapolation is in agreement with experiment. Fitting 
the lattice results to the 0{p^) formula one finds relatively similar values for the fitting 
parameters [151], although the coefficient of the additional term is numerically not very 
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Figure 22: Nucleoli mass from QCDSF and UKQCD [148], CP-PACS [111] and 
JLQCD [147]. 

small: (8ci + Ca + 4c3) -18-25 and 3/(327r2/^) ^ 1- 

It should be noted that with a simple three parameter fit allowing for constant and 
quadratic and cubic or quartic contributions from the pion mass, one can also obtain a 
reasonable description of the lattice data with x^/dof ~ 2. However then the nucleon 
mass extrapolates to a chiral limit of ~ 1.2 GeV. 

In Ref. [155], a good description of nucleon data up to pion masses ~ 600 MeV could 
be achieved using a non-relativistic formalism at O(p^). Octet baryon masses in QCD 
with 2iight + Istrange flavors of sea quarks have been calculated in Ref. [156] at O^p"^) using 
relativistic heavy baryons with infrared regularization. The octet baryon mass has been 
calculated for 3 flavors at O(p^) in Ref. [158], where partially quenched heavy baryon xPT 
has been first formulated. Partially quenched results have been presented for the decuplet 
baryon mass at 0{p^) for 3 flavors in Ref. [157], for 2 and 3 flavors in Ref. [159] and for 3 
flavors at 0{p^) in Ref. [160]. In 2-flavor partially quenched xPT the nucleon mass was 
calculated by Ref. [161], and the nucleon and decuplet masses to O(p^) by Ref. [162]. 

The corrections to the p meson mass within chiral perturbation theory using infrared 
regularization are described in [163]. 

In the following we attempt to describe the finite volume effect using x^T in the p 
expansion. The temporal extent of the lattice is assumed to be infinite. 

On a finite lattice, the integral in the self energy correction is substituted by a sum over 
the allowed modes. The difference between sum and integral is not ultraviolet divergent 
and can be brought into a closed form according to [164] 
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Figure 23: Finite size effect of the nucleoli mass for JLQCD [147] (open circles) and 
QCDSF/UKQCD data (filled circles) [148]. Lines connect sets belonging to the same /5, k 
values. 
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where Ky{x) is a modified Bessel function. 

The contribution of diagram a to the finite size effect is 



5ma{L) = ^^l"^^^}' j^dx 1^ K, [L\n\^mlx^ + mls{l-x)^ . (168) 



The contribution of diagram h is 



4 



Ki{L\n\mps) , K2{L\n\mps) 



(2Ci - C3) — — — h C2 



L\n\mps {L\fi\mps) 



(169) 



The xPT result at O(p^) indeed finds a good agreement with the finite size effects of the 
lattice results. An example for a pion mass around 550 MeV is given in Fig.|2Sl Agreement 
is also found for larger pion masses, rrips ^ 700 MeV. is considered to be an effective 
field theory valid for momentum scales < 1 GeV. Using a cutoff regularization, Ref. [155] 
calculates a validity bound of x^T for pion masses less than 600 MeV. Therefore further 
verification of the nucleon mass formula in comparison to lattice results would be very 
important, e.g. by comparison with simulation data at smaller quark masses once they 
will be available. It should be mentioned that the result is not very sensitive to the exact 
experimental value of tq. 



66 



> 

Q) 1.3 
O 



Figure 24: Chiral extrapolation of nucleon masses on large lattices in two-flavor QCD. 
The crosses denote the lattice results. The dashed line denotes relativistic xPT at 0{p^). 
The plus denotes the experimental value. 

Good agreement of lattice nucleon data at relatively high masses with chiral perturba- 
tion theory is also achieved if a flnite nucleon radius is taken into account using a dipole 
form factor [165]. 

Finite volume effects on J'k and Bk are discussed in [166] in dynamical, partially 
quenched and quenched QCD. Finite volume effects on heavy-light meson matrix elements 
are discussed in [167]. In the next section we discuss formalisms which we used to simulate 
b quarks on the lattice. 

6 HEAVY QUARKS ON THE LATTICE 

HQET predicts 1/M expansions for various mass splittings, decay constants and heavy 
meson decay form factors, where M denotes the heavy quark mass. The expansions are 
useful for guiding extrapolations or interpolations of experimental or lattice results as a 
function of the heavy quark or meson mass. 

Moreover, in simulations of massive quarks using a relativistic action, discretization 
effects of the form 0((aM)"') where n is a positive number, can occur. This effect can 
be avoided by simulating the heavy quarks on the lattice using Heavy Quark Effective 
Theory or Non-relativistic QCD. 
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Figure 25: Volume dependence of the nucleon mass relativistic X-P^ compared with Nf = 
2 lattice results, from Ref. [148]. Solid line: 0{p^) result, dashed line: only 0{p^) terms. 

6.1 Lattice HQET 

The HQET Lagrangian for the static quark can be put on the lattice directly [168]. In 
the rest frame of the heavy meson, the static tree-level lattice Lagrangian is given by: 

^ty*^ = 0. (170) 

Beyond tree-level there is an additive mass renormalization. Corrections to the infinite 
mass limit break the spin symmetry by terms 0{1/Mq). Relativistic corrections are 
incorporated in a systematic expansion in terms 0{Aqcd/Mq). 

The formalism enabled many studies of B meson matrix elements. For one of the 
first very extensive ones see [169] . Recent work using non-perturbative renormalization is 
discussed in [170,171]. 

6.2 Lattice NRQCD 

In NRQCD on the lattice one simulates a discretized version of the action corresponding 
to Eq. ()142|) which includes the terms to the desired order in 1/M. One solves a discretized 
Pauli equation with a covariant coupling to the color gauge field and relativistic correc- 
tions. In this section we summarize some general features of the calculation of masses 
and decay constants in NRQCD on the lattice, in particular based on the discussions in 
Refs. [4,172,173]. 

6.2.1 NRQCD lattice action 

Analogous to the Dirac fermion action, the Hamiltonian ()143|l can be put on the lattice 
by substituting the derivatives by lattice derivatives. The dominant discretization error 
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in the spatial derivative of the heavy quark action is 0{a?pf) ~ Agj^-^^a^, and 0{aMv^ /2) 
and the dominant discretization error in the temporal derivative should be proportional 
to the kinetic energy: -EkinO = p^/(2M)a = AQ^j-,a/{2M). Since in typical NRQCD 
calculations a ~ M~^, the discretization corrections are introduced for consistency in 
calculations which are supposed to be correct at 0(1/M^). Discretization errors in the 
hyperfine interaction are of higher order in the expansion in 1/M and a. 
The lattice Lagrangian density used for Refs. [80, 101] is 

L = 

where the Hamiltonian is H = Hq -\- 6H. The contributions of the 1/M^ terms to the 
low-lying spectrum and Jb is small, therefore also some recent calculations include only 
the 1/M terms, i.e. the kinetic energy and the spin-magnetic operator, (e.g. [174]), and 
using Mq for the bare lattice mass, the kinetic energy operator becomes 

A(2) 

Ho = , (172) 

° 2Mo' ^ ^ 

with the second lattice derivative given by 

A(2) = J2 V,V* = J2 ^i'^- (l'^3) 

i i 

Mq is the bare lattice heavy quark mass. 
SH = ^-a-B 



2M, 







+— ^f-^ (W-E-E-W)- ^T^^o ■ {VE - EV) 



3(Mo)2 V ; 8(Mo)2' 



8(Mo)3 ' 24Mo 16n(Mo)2 ' ^^^^^ 

The parameter n is introduced to avoid instabilities due to large momenta contributing to 
the kinetic energy operator Hq which cause an increase of the correlation function in time. 
For Pi ~ 1/a, the instability can be avoided if n > 3/(2Ma). In all previous calculations, 
the coefficients are set to their tree-level value (one) with mean field improved gauge links. 
The last two terms in Eq. ()174|1 are corrections for the leading discretization errors of the 
kinetic energy operator and the temporal derivative respectively. The fourth derivative 
in the kinetic energy correction is defined as 

A(^) = ^(V.V*)'. (175) 

i 

One can use this action to simulate h quarks directly on the lattice, since 0((aM)") 
discretization errors do not occur. 
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6.2.2 Meson mass in NRQCD 



The exponential decay rate of tlie correlation function in NRQCD or HQET, denoted 
by -Esim; is the unrenormalized meson binding energy whose contributions should after 
renormalization be related to the energies appearing in the 1/M expansion of the HQET: 
the binding energy of the light quark due to the static color source, the heavy quark 
kinetic energy, the spin-colormagnetic energy and so on. 

Also in NRQCD, the inverse propagator at one loop, S{p), is related to the tree-level 
propagator S^^\p) by the formula [172] 

S~\p) = S^'^~\p)-^{p), (176) 

where S(p) is the heavy quark self-energy. By matching the heavy quark dispersion 
relation to the continuum [172], one can find the one-loop on-shell mass and wave function 
renormalization constants. CmMQ is identified with the pole mass of the heavy quark, Mq. 
At one loop. Cm = l + 0{as). The quark mass receives an additive renormalization which 
starts at 0{g'^) and has the form 

Eo = -log(l + «,A). (177) 

This cancels the leading 1/a correction in the exponential decay rate of the heavy quark 
propagator. The wave function renormalization also starts at 0(1), Cq = 1 + 0{as)- 

The mass Mhl of a heavy-light hadron with an NRQCD heavy quark is related to the 
exponential decay rate of the hadron propagator by 

Mhl = Mpert = ^sim + A, (178) 

where A contains the renormalized heavy quark mass and zero point energy 

A = CmMo + Eo, (179) 

with CmMo being the above-mentioned heavy quark pole mass. 

In hadrons with several heavy quarks or antiquarks of bare mass Mq, the mass renor- 
malization should contribute in the same way for each heavy quark, i.e. the quarkonium 
mass Mhh can be calculated using: 

Mhh = E^^,HH + 2A. (180) 

The meson mass can be determined in simulations by considering the momentum 
dependence of the meson binding energy. Formulae p78|) and ()180|) should also hold for 
the mass of excited states and of meson states at finite momentum. Then, the energy 
difference between zero and finite momentum should be described by a dispersion relation: 

E.Up') - ^sim(O) = ^/m^~+^ - Mkin. (181) 
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Mkin is the kinetic mass of the heavy meson. 

Mass sphttings are simply given by the difference of the binding energies of the states. 

If states are close together, statistical fluctuations of the individual correlators may be 
larger than the energy difference. If the overlap with excited states is neglected, one can 
however extract the energy difference from a fit of the ratio of correlators to an exponential 
decay [175] 

^ocexp(-AEt), (182) 

where AE is the energy splitting of the ground states belonging to the correlation functions 
Ci and C2. 

6.2.3 Currents in NRQCD 

If matrix elements of operators are to be calculated in this formalism, the 1/M corrections 
to the operators have to be taken into account as well. To determine for example 
from the lattice, it is necessary to calculate the renormalization factors to match the 
unrenormalized lattice matrix element of the axial vector current to the corresponding 
matrix element in continuum QCD. The matching is discussed e.g. in [173,177], and we 
briefly repeat some features of interest here. The bare current operators at 0(1) in the 
1/M expansion are the same as in full QCD. For the time component of the axial vector 
current this means 

= q{xh,i,Q{x), (183) 
where the heavy quark spinor Q is projected onto the upper two components: 

= ( '^^0^ ) . (184) 

The tree-level corrections occuring at higher order can be derived applying a Foldy- 
Wouthuysen transformation on the Dirac heavy quark spinor. At 0(1/M) this is 

4'^ = -^9(a:)7M757-Vg(a;). (185) 

Beyond tree-level, all operators with the appropriate symmetries can contribute. Some 
of them are related to each other by equations of motion. For practical calculations it is 
sufficient to consider a linearly independent set of operators. At 0(1/M), there is only 
one additional operator and one can choose: 

A? = -^g(x)f V7o7/.75Q(a;). (186) 
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The corrections at 0(1/M^) are given by 



4'^ = ^Qixh.l^V'Qix), 

A? = -^g(x)«-EQ(x), (187) 

with a = 7o7, E = diag(a, a). They have also been calculated on the lattice [80,178,179], 
but without inclusion of the matching renormalization constants. 

In [173], the continuum matrix element is first matched to the matrix element in the 
continuum effective theory which is an expansion in terms of the inverse heavy quark pole 
mass. Each of the continuum matrix elements then equals the renormalized lattice matrix 
element at the same order: 

(4^^) = E^^(<en) = (188) 

i i j 

where the lattice renormalization factors have the structure 

Zij = Sij + as (^2^Cq + Cg)6ij + CmSijSij + Q-^ . (189) 

Because of the mixing of lattice operators, the renormalization constants depend on two 
indices Expanding the renormalization factors in terms of their lattice and continuum 
contributions, one finds [79] 

(A,qcd) = [l + a.Po](4°^) + [l + asPi](4'^) + a.P2(4'^ (190) 
with the one-loop coefficients 

\{CQ + Cq). (191) 

Bi denote the continuum matching factors for HQET to QCD, the lattice vertex renor- 
malizations (which contribute to mixing of matrix elements for i 7^ j). Cm the heavy quark 
mass renormalization and Cq and Cq the heavy and light quark wave function renormaliza- 
tions. For tadpole-improved calculations, additional contributions to the renormalization 
factors have to be included [173]. 



Po 


= Bo 


— Coo 


-ClO 


pi 


= B, 


-Coi 


-Cii 


P2 


= B2 


— C02 


-C12 
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The continuum matrix element is independent of the lattice spacing. Since NRQCD 
results are always at finite lattice spacing, this is only true up to lattice spacing corrections. 
However one can derive an approximate renormalization group equation [173]: 

«T-Ec"(«s,aMo)(Af ) =0 + O(a2), (192) 
da ^-^ 

k 

with C'^(as,aMo) = 1 + UsPk- One can derive the scaling behavior of the coefficient 
functions from this: 

a^5'^ - 7^^'^ a = 0. (193) 

7 is the anomalous dimension matrix of the current operators. Since the coefficients are 
functions of the lattice scale and as, and as is also a function of the the lattice scale, the 
renormalization group equation is given by 



d d 
a— jj- 



da das 

where 7 is the anomalous dimension matrix of the axial vector current in NRQCD, and 
the (3 function determines the running of as with the matching scale, which is inversely 
proportional to a: 

P = -a-^as(a-^). (195) 
da 

In the leading-log approximation [173], 7 is found to be: 

(7^^) = diag(-^,-^,-^) (196) 
\ vr vr vr / 

Therefore there is a simple scaling relationship of the Wilson coefficients as a function of 
the lattice spacing or the inverse heavy quark mass 

^'(a^Mg) _ f as{l/a,) Y''' ..q7^ 
C^(aiMo^) \asil/a,)J ' ^ ^ 

i = 1,2, where the mass on the left-hand side is the bare b quark mass. 



6.2.4 Systematic errors 

Inclusion of higher dimensional operators in the action and currents leads to non- 
renormalizability, and unrenormalized lattice results may contain divergences going as 
powers in 1/a. A continuum extrapolation is therefore not possible. Since results are 
obtained at finite a, they also contain 0(a") discretization errors which have to be ac- 
counted for in the error estimate. To be reliable it is necessary to compare simulations at 
various values of the lattice spacing. One can attempt to make a theoretical estimate of 
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the error size at each lattice spacing, and continuum estimates of physical quantities can 
be taken from simulations at parameters which minimize discretization effects as well as 
1/a effects. The lattice with the minimal error can be used for the continuum estimate. 
The renormalization group arguments of the previous paragraph apply only if a scaling 
region of lattice spacings can be found where the result remains rather constant as a 
function of a. 

To compare lattice results at a fixed heavy quark mass from a range of lattice spacings, 
the leading lattice spacing dependence of the coefficients is given by Eq. ()197j) in a scaling 
window of the resummed, renormalized currents as a function of a. In the calculations 
of [174, 178, 179] it is found that the scaling of the currents resummed using Eq. ()197|) is 
indeed flatter. 

6.3 Relativistic Heavy Quark Actions 

Another method to avoid large discretization errors used in many lattice calculations is to 
simulate heavy quarks at smaller masses (e.g. around the charm mass) and extrapolate 
to the h. If the simulated masses are much smaller than the h mass, this might give 
rise to substantial uncertainties which might however be resolved if the extrapolation is 
substituted by an interpolation to an 0{a) improved static value. Static current matrix 
elements may contain sizable discretization errors as has been observed for the axial 
current in Refs. [79,173]. Renormalization constants for (partly) non-perturbative 0(a) 
improvement have been calculated recently by Ref. [180]. If the simulated heavy mass 
is not very small in lattice units, one might worry about discretization errors also if a 
mass-independent renormalization scheme is implemented non-perturbatively. 

Ref. [181] formulates an improvement program for on-shell quantities in heavy quark 
physics, here briefly denoted as Fermilab formalism or FNAL. An important result of their 
work is the realization that a non- relativistic re- interpretation of the Wilson or clover ac- 
tion leads to smaller discretization effects in simulations of very heavy quarks. Amplitudes 
are polynomials in ap, where p is the three-momentum, while they are bounded functions 
of the heavy quark mass in lattice units. Similar to NRQCD, the rest mass of the heavy 
quark does not affect physical quantities. Relevant are the mass calculated from the 
dispersion relation and mass splittings. An important difference to NRQCD is that the 
continuum limit is well-deflned. We briefly present some of the important formulae for 
massive Wilson and clover quarks from their work, see e.g. [181,182] and [259]. 

Expanding the energy-momentum dependence of a particle around spatial momentum 
zero one flnds 
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where equal indices are summed over. One can define the kinetic mass, M2 by the relation 



(199) 

p=0 



M2 = Ml for a relativistic particle, but at finite lattice spacing the relativistic dispersion 
relation is modified. 

The tree-level relation between the quark mass parameter in the Lagrangian and Mi 
can be derived from the energy- momentum relation of a free Wilson fermion. If the bare 
mass is denoted as m and the Wilson parameter is set to r = 1, this is given by 

1 (am + ia2p2)2 + ^25^2 

cosh aE = l + ? 200 

2 l + am+|a2p2 

with Si = ^ sin api and Pj = | sin(apj/2). The energy at zero spatial momentum equals 

Ml = -ln(l + am). (201) 
a 

Ml, the pole mass, determines the exponential decay of the quark propagator. 
M2 is given by: 



1 2 1 



aM2 am{2 + am) 1 + am' 

sinh(aMi 



or 



aM, = axp(aMO °'""\7';, . (202) 

1 -|- smh(aMi) 

with Ml = M2 only for small am. 

The Wilson and clover Lagrangians already contain the interaction terms of the heavy 
quark Lagrangian at 0(1/M). The temporal derivatives of the Wilson term, combined 
with the -D474 part of the Dirac kinetic term, project onto the upper or lower two com- 
ponents of the Dirac field, see e.g. [182]. The temporal term is therefore the temporal 
derivative of the heavy quark spinor in the HQET formalism. The spatial derivatives of 
the Wilson term contribute to the kinetic term of the heavy quark Lagrangian. The clover 
term can be split up into a color-electric and a color-magnetic contribution. 



Ssw = ^cswi^ I ^ ip^eijkaijBkipx + ^ip^(TQiEiilj^ j , (203) 

\x,ijk x,i 



and therefore contains the spin-colormagnetic interaction term of the HQET at 0(1/M). 
The heavy quark limit of the clover action indeed contains the kinetic and the colormag- 
netic terms of the HQET Lagrangian proportional to I/M2 with coefficients of 0(1) [181]. 
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The D ■ 7 term of the naive Lagrangian also contributes to the kinetic and the spin- 
magnetic operators, however with a factor proportional to Agc^/M where M is the 
heavy quark mass relative to the contribution from the Wilson term [181]. 
The residue of the pole of the free quark propagator is mass dependent: 

Z2 = e"'^^^^ = (1 + am)"\ (204) 

The canonically normalized field at non- vanishing quark mass is therefore given by 

Q. = e"*'^/2g., (205) 

In the hopping parameter representation the canonical normalization for the massive free 
quark is given by 

= v/r^6^a-=^/Vx, (206) 

which approaches the usual norm for Wilson fermions Eq. ()52|) at vanishing quark mass. 
To remove the tree-level p/M mass dependence of the heavy quark field Q, an additional 
application of the spatial derivative is necessary. The complete tree-level field redefinition 
at 0{l/M) is then 

Q; = e'^*^^/2(^l+adi7-^)gx, (207) 

with 

^1 = ^ - Ar^ (208) 

am{2 + am) 2aM2 

i.e. for am = 1, di = 1/12, while for am —^{],di^ 0. 

The Hamiltonian of heavy Wilson fermions can then be written as an expansion in 
terms of the inverse quark mass: 

where the pole mass Mi just acts as an additive zero point of the heavy quark mass. 
The different masses parameterize the departure of the lattice action from the Foldy- 
Wouthuysen transformed Dirac action. As mentioned before, Mb = M2 for clover fermions 
at tree level. In general, the coefficient of the spin-magnetic term is related to the pole 
mass and kinetic mass through = 1 + sinhaMi. 

Ref. [147] discusses a program to improve heavy quarks with am ~ 1 on-shell to 0{a). 
It is assumed that the discretization errors are of the form fn{am){aKQc d)"' with fn{am) 
analytic functions of am. The /o and /i terms are removed in their procedure. In [183], 
results are given for tree-level and one- loop 0{a) improvement. 

Anisotropic lattices with a finer lattice spacing in the time direction are also used in 
simulations of heavy quarks, to reduce discretization errors from the quark mass at least 
at tree level and to allow for more points in the time direction to be used in the fits. 
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Different flopping parameters and cfover coefficients corresponding to time- and spaceiike 
iinlcs are needed, wfiicfi ieads to tfie necessity of tuning many parameters if one wisfies 
to improve the actions. Programs for simuiations of fieavy quarks can e.g. be found in 
Refs. [184,185]. Discussions of earlier work are given in Refs. [186] and [187,188]. 
For reviews about lieavy quark lattice formalisms see e.g. [189,190]. 

7 THE b HADRON SPECTRUM 

The success of NRQCD in calculations of the quarkonium spectrum (for the T, see 
e.g. [100,175,192] and the subsequent studies [87,193-195]; for charmonium see e.g. [196]), 
motivated us to calculate the spectrum of hadrons with b and light quarks using NRQCD 
at 0{1/M) with a ~ 0.1 fm. First calculations of the ground state hyperfine structure, 
a P wave excitation and the Af, baryon obtained a good signal for the lattice data and 
agreement with the experimental spectrum [176, 197]. This was followed by a calculation 
at 0(1/M^) at a ~ 0.1 fm (pubhshed in Ref. [80] in 2000), which obtained good sig- 
nals for many spectral states and resolved several mass splittings for the first time. Also 
preliminary results on the Be spectrum were obtained [199]. The results were supported 
by a study of the B and D meson spectrum at larger and smaller lattice spacings [88]. 
From a further calculation of the Nj = 2 and quenched B spectrum on a range of lat- 
tices using a renormalization group improved gauge action only preliminary results were 
obtained [198]. 

In this section we discuss our results on the 6-light spectrum, in particular Refs. [80, 
198]. 

To demonstrate that a lattice calculation using NRQCD can reproduce the 6-light 
hadron spectrum, we first show the results of [80] compared with the available exper- 
imental data in Fig. [2^1 At the same time, predictions for a number of ground state 
baryons can be made. In the next subsections we discuss the results in more detail in 
comparison with other lattice and with potential model calculations. 

In this section we discuss results on the spectrum of hadrons with b quarks from 
[80,198] and compare with [88,124]. 

The spin operators used to project onto different hadron states in [80,198] are listed in 
Table IHl Operators from two representations of the cubic group have been used, denoted 
as T and E, which correspond to the continuum quantum numbers in the ^Lj notation 
used for heavy-heavy mesons. The baryon operators are constructed using only the upper 
two components of the quark spinors coupled with Clebsch-Gordon coefficients. The 
operators projecting onto the maximal z component of the spin are shown in the Table, 
the other components are defined in an analogous way. For naming conventions of the P 
wave heavy-light mesons see Eq. ()21()j) and of the heavy-light-light baryons Eq. ()211|) . The 
same operators are used for the heavy-heavy-light baryons. For their naming conventions 
see Eq. (ITT^ . 

The lattice scale has been set using rrip. In [80], the b quark mass is fixed using the 
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Figure 26: The B meson spectrum from quenched lattice NRQCD [80]. Long dashed 
hues denote the experimental error bounds from the Particle Data Book [146], dotted 
lines experimental error bounds from DELPHI [200]. We use the collective notation B** 
for the P states. 



spin-averaged B{1S), in [198], both the B and Bg masses were used for consistency. More 
simulation details can be found in the respective papers. 



7.1 S Wave States 

The splitting between the B^ and B^ mesons is expected to be dominated by the strange 
and light quark mass difference. The spin-averaged splitting depends in addition only on 
the difference of heavy quark kinetic energies in the cloud of gluons and relativistic light 
or strange quarks. In Fig. |2Z|on the left, the spin-averaged energy difference from [80] is 
plotted for various heavy meson masses. The dependence is relatively fiat, which is also 
found in experiment. The experimental Dg — D splitting is roughly 10% higher than the 
experimental Bg — B splitting. 

In HQET and NRQCD, the hyperfine splitting is generated by the spin-magnetic 
interactions. The leading mass dependence should therefore be ~ 1/ Mq or ~ 1/M where 
M denotes the spin-averaged meson mass. In Fig. 123 on the right we show that results 
indeed are consistent with this expectation for strange as well as for chirally extrapolated 
light quark masses. 

Results are shown in Table ^1 We find the results to be 30 — 50% lower compared to 
the experimental values of 45.8(4) MeV and 47.0(2.6) MeV [146] for the B and Bg splitting 
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operator 


-1 _L /I T— » \ 

1+ (ipl) 

2 ('P2)(^) 


u Q 
u a Q 
u AQ 
u a ■ A Q 
u a A Q 

u {aiAj + a,jAi) Q 
u {cTiAi - (JjAj) Q 


Aq (s, = +1/2) 

Sq (s, = +1/2) 
(s, = +3/2) 
{s, = +1/2) 


u'^d Q^ 

{u'^azd — v^Wa+d Qi)/V^ 
u'^a+d 

{^/2u^a,d Q^ + u^a+d Qi)/V3 



Table 8: Ground state S and P wave heavy-light meson and S wave baryon operators 
used in [80]. 

respectively. An error source could e.g. be radiative corrections to the hyperfine coefficient, 
since the spin-magnetic coefficient in the calculations discussed here is only set to its tree- 
level tadpole- improved value. Results from other quenched lattice NRQCD calculations 
of the B*-B and B* - P, splittings fixing the scale with (e.g. [88, 178, 198,201]), and 
ro = 0.5 fm [202] are around 25 - 35 MeV. Refs. [80,88], and [178,203] do not resolve 
a lattice spacing dependence of the hyperfine splitting. Using relativistic 0(a) improved 
heavy quarks, Ref. [89] obtains a splitting around 10 — 20 MeV with an error of 10 MeV. 
The quenched calculation of [241] gives AE{B* — Bs) ~ 40(10) MeV, but using a different 
quantity to set the scale, i.e. the charmonium IP — IS splitting. 

In the CP-PACS preliminary result of [198] we do not find a significant unquenching 
effect. The unquenched result is with 33 MeV still much lower than experiment. The sea 
quark mass is around 3x the strange quark and the lattice is quite coarse, and perhaps 
the observed absence of unquenching is due to a large systematic error. An NRQCD 
calculation with three dynamical flavors [204] and staggered light quarks does find the P* — 
Bs splitting to increase from 25(5) to 42.5(3.7) MeV, in agreement with the experimental 
value. An appropriate choice of gauge action might improve quantities which could be 
sensitive to discretization and perturbative effects. E.g. for quarkonia, there is a ~ 50% 
variation of the quenched T — rjh splitting between simulations using only different gauge 
actions, as shown in Fig. EHl Using relativistic quarks on quenched anisotropic lattices, 
Ref. [187] finds a ^Si - ^Sq splitting around 50 - 60 MeV which is in agreement with the 
potential model result of [205] of 60 MeV and close to the NRQCD result given in [195]. 
There is one experimental candidate for the rjb, which is roughly 160 MeV lighter than the 
T. It is quoted to have around 20 MeV statistical and 20 MeV systematic errors [146]. 
For a detailed review of theoretical and experimental information on quarkonium physics 
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Figure 27: S wave meson splittings from [80]. On the left, the Hg — H, on the right 
the H* — H hyperfine splitting. Circles denote results for mesons with a light quark 
[H = Hu,d), diamonds for mesons with a strange quark [Hs). The solid line denotes a 
linear fit to the [H) mesons with the errors denoted by dashed lines. The experimental 
value for the B (and on the right also for the D) meson is given by a burst. 

see [206]. 

7.2 Orbitally Excited B Mesons 

The spectrum of P wave B mesons is also e.g. of interest for the experimental study of 
CP violation in B[B) experiments. Resonance structures will be noticeable in the decay 

B* B^TT. 

The fact that in the limit M — > cx3 the heavy quark is a static source of a color field in 
which the light quark moves relativistically, motivates to treat heavy-light mesons similar 
to the hydrogen atom. The energy levels can be characterized by the light quark angular 
and radial quantum numbers. The light quark angular momentum, denoted as ji, is then 
coupled with the heavy quark spin to the angular momentum of the meson, J. This can 
be sketched as follows: 
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(210) 



Experimental knowledge of the P state level structure is still sparse. The Particle Data 
Group [146] lists two candidates. One is a relatively broad resonance, the 5}(5732), which 
is believed to come from several narrow and broad P wave states, and the relatively nar- 
row B*j{5850) which can be interpreted as stemming from excited Bg states. Recent 
experimental results from DO [238], who measures the Bi mass and for the first time finds 
separate signals for the Bi and the -Bg, and from DELPHI [200] are given in Table 01 
Heavy-light mesons can be characterized by the quantum numbers J^. The lattice oper- 
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Figure 28: Quenched T — rjb hyperfine splitting from NRQCD at 0{Mv'^) and 0{Mv^). 
From left to right, the points are from Refs. [194] at a ~ 0.08 and 0.06 fm, and [87,193,195] 
(from left to right). Different symbols denote different gauge actions. The scale is set 
using the T' — T splitting except for [87, 193] which prefers to set the scale using the 
Xb — T splitting to set the scale but find the ratio Am(T' — T) / Am{xb ~ ^) within errors 
in agreement with the experimental value already in the quenched approximation. 



Ref. 


state 


M or AE[MeV] 


DELPHI [200] 


Bi 


5732(20) 




b; 


5738(14) 




b:^ 


5852(5) 


DO [238] 


Bi 


5724(4) (7) 




B* - B, 


23.6(7.7)(3.9) 


Particle Data Book [146] 


5}(5732) 


5698(8) 




i?:,(5850) 


5853(15) 



Table 9: Experimental measurements for P wave B meson masses. 



ators corresponding to the ^Pi and ^Pi states in the LS coupling scheme both onto the 
= 1+ states. Since in the correlators of ^Pi ^Pi errors are very large, we cannot 
disentangle the splitting between the ^P states. If one assumes that the energy difference 
of the ^P states is very small (potential model calculations give e.g. ~ 40 MeV [191]), we 
can get a reasonable estimate of the spin-averaged J = 1 energy [239] with the present 
lattice data. Since the spin-average of the measured ^P energies is close to the energies 
in the ^Pi channel, the ^Pi might be used as an estimate for the spin-average of the P 
states. 

The spin-averaged P — S splittings from our calculations are given in Table nUl 
Refs. [197] and [124] give results for ^Pi — S splittings which are larger than 500 MeV. 
The former use a very small set of configurations. The latter might have increased finite 
volume effects because they are unquenched. Our lattice simulations (see Table HOI) are in 
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Ref. 


NRQCD 


Bs — Bd 


B* - B 


P-S 


b;-b 


75* 75* 
B2 - Bq 


iV; = 


light 


[80] 
[88] 
[88] 
[198] 


1/M2 

1/M2 
1/M2 
1/M 


90(9)(gi) 
84.1(20)(n 

109(26)aD 


24(5)(i) 
29.5(15)(fD 

32(3) 


457(31)(i) 


526(45) (ID 
426(17) 


155(32)(?1) 


strange 


[80] 
[88] 
[88] 


1/M2 

1/M2 
1/M2 




27(3)(i) 
28.3(10)(?^) 
27.3(20)(i) 


428(27) (ID 
385(70) (iD 

395(45) (ID 


493(30) (i) 


136(23)(1°) 

41(94)(1D 

179(65)(D 


Nf = 2 


light 


[124] 
[198] 


1/M 
1/M 


103(8)(?D 


33(2) 




461(14) 





Table 10: Summary of own results on the b meson spectrum. The light quark content and 
the flavor number are denoted in the headers. The errors are taken over from the papers, 
adding the systematic errors in quadrature where applicable. 

general in agreement with the preliminary experimental results, although statistical and 
systematical uncertainties are both large. We attempt to give an estimate of discretization 
and perturbative errors in Section 17.61 

An important result of [80] is to resolve the B2 — Bq energy splitting on the lattice 
for the flrst time. In the scaling study of [80,88] (see Table ITUj) . there are indications of 
discretization errors in the P states. 

A source of systematic error in our calculation [80] are flnite volume effects, although 
the lattices are quenched, since the lattice extent is only around 1.6 fm. Ref. [240] flnds 
that the Bethe-Salpeter wave function of a static-light P wave meson decays to around 
1/4 of its nodal peak at a distance of 1.6 x tq, which might affect the energy levels. 
CP-PACS [198] flnd a lower quenched result for the B2 mass. An advantage of their 
calculation is that the volume is larger, on the other hand the discretization effects are 
expected to be larger too, and the NRQCD action is only 1/M improved. 

A comparison of results from the lattice, model calculations and experiment for P — S 
splittings of the Bg meson is given in Fig. 1221 

There are only few dynamical simulations so far. In unquenched calculations, the 
ji = 1/2 lattice states can decay into or mix with S wave Btt or B*7i systems if the mass 
of the P wave states is above or near the threshold. This would result in a shift of the 
meson mass. Decays of the other two P states occur at higher angular momentum and 
are therefore suppressed. 
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Figure 29: Comparison of splittings between P wave and the ground state of Bg mesons 
fromRefs. [80,88,202,203,234,235] (lattice) and Refs. [191,242,248] (models). The lines 
denote the experimental value of the narrow 5*j(5850) resonance believed to come from 
orbitally excited Bg mesons [146]. 

Among potential model calculations (e.g. [191,242,243,248]) the sign of the B2 — B^ 
mass difference is disputed, see Table IT^ and Fig. 1211 Some lattice calculations [203,234] 
and [88] at a ~ 0.2 fm find a small positive splitting which is within errors compatible 
with the potential models predicting small negative splittings. In Ref. [234] the spin- 
chromomagnetic interaction term of the NRQCD action is not taken into account in the 
calculation. From potential models one would expect the P wave hyperfine splitting to be 
close to ~ 10 — 20 MeV [191]. Potential model calculations would predict [191,245] that 
including the hyperfine splitting rather increases the B2 — Bq mass difference instead of 
decreasing it. In Refs. [80,235], and on the finer lattice of [88], a positive splitting with 
2a significance or more is found. Refs. [80,88,235] are in close agreement with each other, 
and the static value of Ref. [235] is in agreement with the M — s> 00 extrapolation of [80] . 
The present lattice average of the B2 — Bq splitting is positive with a central value around 
100 MeV. The chiral extrapolation of the splittings has still large errors, but since the 
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Figure 30: Comparison of P wave fine splittings of Bg and B mesons from Refs. [80,88, 
203,234,235] (lattice) and Refs. [191,242,243,248] (models). The dashed line is to guide 
the eye. 

potential model predictions for the (j/ = 3/2) — (j/ = 1/2) splitting depend on the ratio 
of heavy quark masses, an inverted level ordering would already have been predicted in 
the static-strange system. 

7.3 b Baryons 

According to potential models, baryons with one b and two light quarks can be described 
as a light diquark of spin zero or spin one which then couples to the heavy quark [246]. 
The resulting states are: 

|: A, 

S*, El ni (211) 

St,, .TLfo, fifo, 

where the quark content is {lib), {Isb), [ssb) from left to right. Own results for the baryon- 
meson splittings are summarized in Table El A summary of lattice and model results is 
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Figure 31: The Aq — H splitting from quenched lattice NRQCD [80]. The experimental 
values for the Ab — B and A^ — D sphttings are denoted by bursts. 

given in Table El 

The spin-averaged splitting A;, — B is plotted in Fig. |^ Since in the quark model 
the light quarks in the Af, are in a ground state with total spin zero, the expectation is 
that the splitting is dominated by the constituent mass of the additional light quark. The 
small slope of the splitting suggests that the heavy quark kinetic energy in the baryon 
is close to that in the meson, for experiment as well as for the lattice data. In Fig. |^ 
we compare quenched and Nj = 2 lattice results and results of model calculations on the 
Afy — B splitting. On quenched lattices, the result is ~ 2cr higher than the experimental 
value. Preliminary results with two flavors of dynamical quarks around the strange quark 
mass [124, 198] are even higher. On coarse lattices, we find an increase by ~ 15% if Nf is 
changed from zero to two dynamical quarks of around 3x the strange quark mass [198]. 
Chiral extrapolation might be a source of uncertainty. Ref. [247] estimates the mass 
splittings of unquenched heavy baryons as a function of the light quark mass using x^T 
and predicts deviations from linear behavior. 

500 
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I 300 

< 

IW° 200 
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Figure 32: Mass dependence of the Sq — Aq splitting using quenched lattice NRQCD 
from [80] and [236]. The spin-averaged meson mass is denoted by M. 
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The splitting between the spin-averaged ji = 1 and the ji = state is also expected 
to be relatively heavy quark mass independent. This is indeed found in the lattice data 
shown in Fig. 1221 The light quark mass dependence of this splitting is very small. We 




l/M|GeV"'] 



Figure 33: The Sg — Sq splitting from quenched lattice NRQCD [80]. The solid line 
denotes a fit to the three lightest meson masses. Dashed lines show the fitting errors. 

now try to compare it with the results from other lattice calculations and results using a 
constituent quark model [145], a Skyrme model [244] where the baryon is described as a 
bound state between a soliton and a heavy quark, and a quark model with a relativistic 
description of the light quarks [246] . Ref . [244] predicts the relation between the heavy- 
light and light baryon splittings to be AM (Sq - Ag) / AM {A - N) = 2/3, for Q = c,b. 
For Q = c, the equality is experimentally well satisfied. The lattice results for the ratio 
with Q = b still have large errors and vary between 0.5 and 1. 

In Fig. EHl above we summarize lattice results for the spin-averaged — splitting. 
Where results from one collaboration at several lattices are plotted, the leftmost point is 
from the coarsest lattice. 

The expectation from HQET for the Eg hyperfine splitting (for baryons containing 
one heavy quark Q) is that it is generated by the spin-chromomagnetic interaction (see 
Eq. ()125j) ). It should therefore in first approximation be linear in the inverse heavy quark 
or meson mass, up to logarithmic corrections. The data is compatible with this within 
errors, as shown in Fig. EHl On larger lattices and with 0(a)-improved light quarks, 
Ref. [236] finds a similar behavior. Their results for the b meson and baryon splittings are 
in general agreement with ours [80]. Some trends in the data are slightly different. For 
example, their ratio {B* — Bs)/{flb — ^b) is closer to 3.5 instead of our value of ~ 1.5, 
and their Sf, — A;, splitting is smaller; however errors are still large. 

Assuming that the — splitting is proportional to 1/Mq up to logarithmic cor- 
rections, one can rescale the experimental value of AM(Sc — Sc) — 67 MeV by the ratio 
of b and c quark masses and finds AM(Sb — T,b) ~ 20 MeV, a value compatible with the 
lattice results. 
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Figure 34: Aj, — B splitting from quenched (left) and Nj = 2 (middle) lattices (Refs. [80, 
124,198,201,233,236]). The lattice average shown on the right is quenched. CP-Q and 
CP-2 denote quenched and unquenched results from [198], respectively. Model results are 
from Refs. [145,244]. 



7.4 bb Baryons 

It is expected that the level structure of baryons containing two b and one light quark can 
be understood as arising from a b diquark color antitriplet with size much smaller than 
A-QCD which interacts with the light quark in a similar way as the heavy antiquark in a 
heavy-light meson. 

7,-0 7-1.^/0 r>/- 

„■ _ 1 / ^ = 2 ■ ^bb^^bb ' ^bf) ' (212) 

- ^ \ 7 _ 1 . ^0 ^- O- 

~ 2 ■ ^bb^ "fob' 

For two identical heavy quarks the ji = state does not occur because of the Pauli 
exclusion principle. The three states for each J value with ji = 1 correspond to the quark 
content {bbu),{bbd) and (bbs) respectively. The mass difference between the J = 1/2 
and J = 3/2 states is a hyperfine splitting depending on the heavy and light quark spin 
orientation. Its leading behavior is expected to be linear in the inverse heavy quark or 
meson mass up to logarithmic corrections, which is also found in the lattice data shown 
in Fig. ESI 

The results are: 



S*,-S,,: 20(6)(tl)MeV 
20(4)(2)MeV 
Ebb : 10314(46) (;i?) MeV 
Qbb : 10365(40) (;i^) MeV 



The lattice calculation of Ref. [201] quotes very similar values for the hyperfine split- 



87 



300 - 

> - 
S 200 - " 



100 



30 



20 ■ 



10 - 




GLOK Mathur JLQCD UKQCD 



GLOK Mathui- JLQCD 



30 ■ 

> 

S 20- 

s ■ 

10 - 



E. -A. , models 

b b 




Jenkins Capstick Ebert Lat avg 



E. -Z. , models 

b b 



Capstick Jenkins Ebert Lat avg 



Figure 35: Above: S^ — A;, splitting from the lattice (left) and model calculations (right). 
The dashed line shows the experimental value for the Sc — Ac splitting. Below: Sb* — J^b 
splitting. Lattice data from Refs. [80,201,233,236], model results from Refs. [145,244]. 



tings of doubly heavy baryons, namely 22(6) (3) MeV and 20(6) (4) MeV with spatial lattice 
spacings of a = l.lGeV"^ a = 0.9GeV~^ for Am(S^^ - Ebb)- For Am(i7^j, - Qbb) they find 
20(4) (3) MeV and 19(4) (3) MeV respectively. 

The result for the Qbb mass is in agreement with the potential model calculation 
of [282]. The hyperfine splittings are slightly smaller than the potential model results of 
~ 35 MeV which might be due to a quenching effect in the lattice results. The potential 
model predicts the Ebb mass to be around 100 MeV (2cr) lower than the lattice result. 

Ref. [283] finds a value for the El^^ — Ebb of 34(4) MeV within the framework of potential 
NRQCD, in agreement with the potential model [282]. The most plausible explanation 
for the difference is again that the lattice result for the splitting is slightly underestimated 
in the quenched approximation. 

7.5 The Be Spectrum 

If the mass of the u quark in the meson is increased towards the charm quark mass 
one obtains the meson. In Fig. EH we summarize results on the B^ spectrum with 
clover charm and NRQCD b quarks [199]. The lattice scale is set using the charmonium 
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Figure 36: Eqq — Eqq (circles) and Eqq — E'qq (diamonds) splitting from quenched lattice 
NRQCD [80]. The solid lines denote linear fits, 
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Figure 37: Be level splittings [199]. The error bars are purely statistical. 



P — S splitting from clover quarks [252] . The B^. mass can be calculated using the non- 
relativistic interpretation for the clover charm as well as for the NRQCD b quarks and 
relation (HHOl), 

clover, J?) 

(213) 

with the result 6.3(1) GeV, in agreement with the experimental value of 6.4(4) (1) GeV 
from [146], and a recent measurement from CDF [253] of 6287(5) MeV. 

The quenched calculation of [255]) uses the same lattice actions as our calculation, 
however at a finer lattice with a ~ 0.06 fm if the Xb ~ ^ splitting is used to set the 
scale. The bare quark mass of their heavier quark is fixed to the h while the charm 
quark mass is reached using an interpolation. They find a Be ground state mass of 
6.386(9) (99) GeV, where the first error is statistical and the second is an estimate of 
uncertainties due to quenching and estimation of M^^, added in quadrature. Determining 
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the charmed n value by fixing the ratio rriDj {^m{T' — T)) to the experimental value 
and interpolating the ground state energy to the charmed k value, we obtain an estimate 
of the energy splittings to the ground state from their data. Their P — S mass splittings 
are in general larger than the lattice ones by ~ 100 MeV (or 20 — 25%). In particular, 
they find Am(i?*Q — i?c) ~ 500 MeV compared to our value of ~ 360 MeV. The difference 
can to a large extent be explained by the fact that in the quenched approximation the 
lattice spacings from the T are larger than the ones from charmonium. At /? = 6.0, the 
difference between the value used in our calculation and the Xft ~ ^ splitting is ~ 20%. 
The B*2 — B*Q splitting is found to be around 80 MeV in both calculations. An earlier 
calculation of [254] uses NRQCD for h and c quarks on a coarser lattice. The bare h 
and c quark masses are fixed using hh and cc quarkonium splittings respectively, using in 
each determination the different values for the lattice scale from mass splittings of the 
respective quarkonium. They fix the lattice spacing by assuming that the experimental 
value for the B^. P — S splitting agrees with the one for the T and obtain a value similar to 
the lattice spacing from the T' — T splitting. Their results for the P — S mass splittings are 
between ours and the results from [255], while their P state splittings tend to be smaller, 
e.g. Am(i?*2 — -B*o) ~ 50 MeV. A recent simulation using three Nj = 2 + 1 lattices ranging 
from 0.1 < a < 0.2 fm with NRQCD h and Fermilab c valence quarks finds a Be mass in 
agreement with the previous results: m^^ = 6.304(12) (q*) GeV, where the first error is 
a quadratic sum of statistical and systematical errors, and the second is an estimate of 
the heavy quark discretization uncertainty. The results for the ground state masses are 
therefore not very sensitive to quenching effects and support the idea that heavy quark 
energies are physically relevant. Once the quark mass has been fixed, e.g. by matching 
the dispersion relation of one meson to the physical value, the mass of another meson 
containing the same heavy quark can be predicted reasonably well. 

7.6 Comparison to other work 

In this Section we give technical details of the other lattice calculations occuring in the 
comparisons in the previous sections and describe how the lattice averages presented in 
the previous sections were determined. 

Results on P wave states and baryons from other authors are summarized in the 
Tables fTT| IT^ and IT^ Masses averaged over spin-orientations (spin-averaged) are denoted 
by an overbar. The first error on the individual lattice results includes statistical errors 
and uncertainties fixing the masses to the physical values, the second, where applicable, is 
a chiral extrapolation uncertainty. Comparison plots of the results, showing the statistical 
errors only, were already shown in Sec. 

For comparison to experiment and to model calculations one would like to have an 
estimate for the 'lattice average'. In [249] it is attempted to calculate the error using a 
similar method as for the averages of experimental data [146]. The assumption is that 
sets of lattice calculations using similar methods have the same systematic error, or that 
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the systematic error of sets of lattice results can be divided into a common part and a 
rest which is treated as independent. The common part is taken to be of the order of the 
error of the calculation with the smallest uncertainty. The error on the average is rescaled 
by r = a/x^/ ^ 1) if r > 1, where is the number of results. 

Among the states we could resolve, we found in the previous subsection that the 
spectrum changes only little if the NRQCD action is changed from 1/M to There 
have been calculations of the effect of higher NRQCD corrections on ground state B 
mesons. Ref. [203] quotes a 6% increase of the B* — B splitting when the action is 
changed from 1/M to 1/M^, and, depending on the exact discretization of the NRQCD 
action, a — 4% decrease when the action is changed from 1/M^ to 1/M^. For the Bq —B 
splitting, Ref. [203] finds no change between 1/M and 1/M^ and a — 6% increase from 
1/M2 to l/M^. 

Most of the calculations use NRQCD, except for [233] which uses heavy clover quarks 
and [235] which simulates Bg mesons in the static approximation and interpolates between 
the static and experimental Ds mesons. Discretization errors with non-perturbatively 
0(a) improved clover light quarks (finer lattice of [88] and [236]) are 0{a?K^(jQ), whereas 
the tadpole-improved light clover action has 0{a?K^(^jj) and O^asaKqcD) errors (coarser 
lattice of [88] and [80,198]). Refs. [201,203] use 0{a?) tree-level tadpole-improved clover 
light actions respectively. Ref. [202] uses staggered light quarks. 

The NRQCD action has errors 0{as^QCD/ from corrections to the spin- magnetic 
coefficient. Errors on spin splittings are treated as being dominated by an error on the 
spin-magnetic coefficient of 0{as) ~ 20 — 30%. 

It is obvious that this way of estimating the systematic error is rather schematic. Using 
this method one would expect that the discretization errors of the light hadron spectrum 
with staggered quarks (0(a^)) are rather smaller than those of tadpole-improved Wilson 
quarks 0{a?) and O^agO) at the same lattice spacing. In calculations of the light spectrum, 
discretization errors of staggered quarks appear not so small. However this is the usual 
method to estimate the systematic error in lattice calculations of weak matrix elements 
using a non-relativistic effective theory. 

The scale is determined using except in Ref. [198] which sets the scale with ^/a = 
427 MeV. Refs. [235] and [234] use tq with physical values of 0.525 and 0.5 fm respectively, 
and [202] uses tq = 0.5 fm and quarkonia at and around the charm [237]. All other 
calculations from the set discussed here use rrip. 

The variation due to the 10% ambiguity between using a from rrip and a from ro = 
0.5 fm in the quenched case, and asymmetric chiral extrapolations where applicable, are 
included in the second error on the averages. The Xft ~ ^ mass difference is not included in 
the estimate of the scale variation since it gives values for spin-independent mass splittings 
which are much higher than experiment. For example, Ref. [124] quotes an increase of 
the Bs — B and the Ab — B splittings from 98 and 560 MeV by ~ 20% if the scale is set 
with the Xb — ^ sphtting instead with nip. The experimental values are 90 and 345 MeV 
respectively [146]. 
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For the error estimates we use nominal values of A.qcd = 400 MeV, M = 5 GeV 
and as = ay(l/a), where av{q*) is defined in the potential scheme [98,102] at the scale 
q* of Eq. (f75j) . Since the lattice results have rather varying central values we do not 
calculate the error in percent of the individual lattice splittings but of the experimental 
splittings or nominal estimates thereof. We use guessed values where no experimental 
value is available: 400 MeV for B* - B and B^q - 5„ 500 MeV for B* - B and B^^ - B^ 
and 100 MeV for B* - B*. 

The lattice results for Bq — B splittings are summarized in Table ^2 for B2 — B 
splittings in Table and for B2 — Bq splittings in Table ^| 

In Table IT^ we summarize results for the spin-independent splittings Ai, — B and S;,— 
with M{%) = [2M(Sfe) +4M(Eb*)]/6 and M(B) = [3M{B*) + M{B)]/4, and the S^* - 
hyperfine splitting. 

8 /b with non-relativistic heavy 

QUARKS 

Since B meson decay constants are experimentally undetermined, it is important to cal- 
culate them on the lattice. Their lattice calculation also serves as the calibration of the 
lattice methods for the calculation of more complicated matrix elements. 

To implement the heavy quark directly at the b quark mass on the lattice, we used 
NRQCD. The first result was an unquenched {Nf = 2) study of the tree-level contribu- 
tions to the pseudoscalar and vector decay constants at 0(1/M) [217] (1997). Then (1997) 
we published a quenched calculation at a lattice spacing a ~ 0.1 fm including most of the 
renormalization factors at 0{as/M) [79], accompanied by a calculation using two fiavors 
which was published in 1999 [124]. To estimate the effect of the truncation of the 1/M 
expansion, and including a very broad range of states in the spectrum, we performed a 
quenched calculation which included 1/M^ terms at tree-level [101] (1998). Renormaliza- 
tion constants correct to 0{as/M) were included. Since our previous unquenched lattice 
calculation [124] used only one lattice spacing, a small physical volume and a different 
sea and valence quark action, the major remaining systematic uncertainty appeared to 
be quenching. We then performed a calculation comparing two and zero flavors, with 
equal sea and valence quark actions, at a range of lattice spacings. NRQCD and Fermilab 
heavy quarks were simulated in parallel to minimize remaining systematic uncertainties 
from the heavy quark action (see Refs. [137,174] of 2001). 

In this section, we discuss the calculations of Refs. [79,101,137,174], including sys- 
tematic effects, and compare to other lattice results and results using different meth- 
ods. A typical analysis procedure is discussed in the context of the unquenched simula- 
tions [137,174]. 
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8.1 Quenched calculations at an intermediate lattice spacing 

8.1.1 NRQCD at 0(1/M) 

Matching of the current matrix element was performed at one loop perturbation theory, 
using the strong coupling constant «<j in the potential scheme [98]. It is plausible that 
the matching scale q* is in the region of the lattice cut-off 1/ a < q* < tx / a. A value of 
aq* = 2.18 for the matching of static-light axial currents has been calculated by Ref. [218], 
while the newer calculation of Ref. [219], which retains parts of the momentum integral 
which were set to zero in the calculation of finds aq* = 1.43. Renormalization constants 
for Refs. [79] and [101] were calculated for g* = 1/a and q* = n/a. 

Numerical results presented here are from the calculation with tadpole-improved clover 
light quarks (Run B) of [79]. Since the mass renormalization constant Cm was not included 
in the renormalization of the matrix element J^^^ in [79], we list f\/M as a function of 
the heavy quark masses in Table El including the renormalization. Along with this we 
summarize the bare heavy quark masses in lattice units, the mass shift A^l and the 
strange meson masses in lattice and physical units, nip was used to set the scale. The 
values for the shifts are calculated from the quarkonia shifts Ahh with a 1/M^ NRQCD 
action. Values for aMo = 1.71,2.0 are taken from [175] and were newly calculated for 
aMo = 4.0. For oMq = 8.0, perturbation theory was used. To determine the errors on 
the meson masses in Table El the individual errors were added in quadrature. 

aMo = 2.0 can be used as a reasonable estimate of the bare b quark mass. A linear 
fit of the meson mass to the form ^^j^ = A + ^ for the three smaller mass values of 

Mo Mo 

Table Uni using uib^ = 5.37 GeV as physical input, gives cMq = 2.12. 

Using the matrix elements at aMg = 2.0 and including the mass renormalization 
constant from [173] which was not used in the matching calculation for [79], we find for 
the decay constants: 

q* = 1/a : fe = 182(28) (27) (16) MeV, 

fn^ = 207(8) (30) (17) MeV, 
q* = Tr/a: Jb = 190(30)(29)(16) MeV, ^ ^ 

fB, = 216(8) (32) (17) MeV, 

where q* denotes the scale of used in the matching calculation. We quote the same 
error bars as the aMo = 2.0 result of Ref. [79]. The first error denotes the statistical and 
fitting error, the second the uncertainty in the determination of a and the third, the error 
due to perturbation theory and higher orders in the 1/M expansion. 

8.1.2 NRQCD at 0{1/M'^) 

For the results from [101], the current operators were included up to 0(1/M^) while 
the matching renormalization was calculated only to 0{as/M). The decay constants are 
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found to be: 

fB = 147(11)(1?2)(9)(6) MeV, 

fs. = 175(8)(lIo)(ll)(r)(^o') MeV. ^'^"^ 
The first error is statistical, the second reflects the variation if the lattice spacing is 
changed between = 1.8 and 2.0 GeV. The third error gives the estimate of higher 
order perturbative corrections (estimated from the variation of the one-loop result using 
q* = 1/a and ir/a), and uncancelled corrections (estimated as ^^/(aMo)^), and the 

fourth error bar is an estimate of the 0(0^ Aq^^) = 4% discretization errors. 



8.2 Investigation of Unquenching and Scaling 
8.2.1 Technical details 

First we discuss a calculation of fs and fs^ using NRQCD at 0{1/M) at various lattice 
spacings, with quenched and dynamical {Nf = 2) renormalization group-improved [56] 
gauge fields at high statistics. The current and mass renormalization constants were 
calculated by Ishikawa et al. [251]. B decay constants have been calculated with dynamical 
configurations at two (3 values corresponding to lattice spacings of respectively = 

1.3 and 1.8 GeV in the chiral limit of the sea quark mass and at four values of the 
sea quark mass each to study an extrapolation to physical u, d sea quark masses. In 
parallel, quenched configurations at 4 (3 values are studied at a comparable range of 
lattice spacings. using five values of heavy and five values of light valence quark masses 
for reliable interpolations and extrapolations to the physical quark masses. Correlations 
functions are smeared at the source. The contributions to the bare decay constants are 
related to the matrix elements of the local lattice operators A^^^ which contribute to the 
time component of the axial vector current by 

a'/'fl^^/M^=^{0\Af{0)\B{p = 0)), (216) 
V Mb 

where Mb is the B meson mass. Assuming a relativistic normalization of the meson states 
\B{p)), one can define a local amplitude by 

Z« = -^{0\AfmBip = 0)) ^ a'/'fl^V^=V2Zi\ (217) 

The normalization of the pseudoscalar operator used in the simulation for the B meson, 
which is in general smeared and denoted as P, is taken into account by measuring the 
pseudoscalar density 

Zs^^^{0\PmBip = 0)). (218) 

The following correlation function of the current ^4 ^ is calculated on the lattice: 

{0\4\t,p)\P{p = 0)) = Z«Z5e-^-*, (219) 
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Figure 38: 'Effective mass' of heavy-light pseudoscalar-pseudoscalar correlation functions 
on a dynamical lattice at the smallest value of the gauge coupling, the lightest sea quark 
mass, the valence heavy quark slightly heavier than the b, and the lighter valence quark 
mass close to the strange quark mass [174]. On the left, the source and the sink are 
smeared, on the right, a smeared source and a local sink are used. 

inserting a complete set of mass eigenstates and assuming that the decay is dominated 
by the smallest mass at large times. The exponent is the same for the various A^^\ 
The pseudoscalar density Zs can be extracted from the correlator of the smeared meson 
operator at source and sink: 

{0\P{t,p)\P{p = 0)) = Z|e-^— *. (220) 

Masses and amplitudes for the current matrix elements are extracted from single ex- 
ponential fits to smeared-smeared and smeared-local correlation functions. Note that for 
NRQCD heavy quarks the correlator {Af^\P) is identical to {P\P) since spinors where 
only the upper components are nonzero do not distinguish between 70 and the unit matrix. 
Plots of typical heavy-light correlation functions are shown in Fig. EH 

The bare heavy quark mass corresponding to the b meson mass is determined from the 
difference of energies at finite and zero spatial momentum, Eq. ()18H) . with |p| = 27r/L, 
where L is the spatial lattice extent. Since the energy difference is small, the states are 
expected to be strongly correlated, and the splitting was extracted using a one exponential 
fit to the ratio of finite and zero momentum correlators (Eq. (jl82j) ). That the result is 
rather independent of the momentum can be seen in Fig.HHlfor the finest quenched lattice. 

The meson masses are fitted as a linear function in the VWI light quark mass to find 
the physical points corresponding to the light and strange quark. The result is fitted as 
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Figure 39: Dynamical meson masses from [174] at the smaller value of the gauge coupling 
{(3 = 2.1) and a sea quark mass around the strange. On the left, fit of the meson mass as 
a function of the VWI light valence quark mass. The heavy quark mass is close to the b 
mass. On the right, fit of the heavy-strange meson mass as a function of the heavy quark 
mass to determine the b quark mass. 

a linear function in Mps/Mq. The heavy quark mass parameter corresponding to the b 
quark is determined by setting the heavy-strange meson mass to the Bg mass. Setting 
the heavy-light meson mass to the physical B meson mass gives the same result, however 
with larger error. Fits of the meson mass as function of the light and heavy quark mass 
are shown in Fig. EHl 

In Fig.l^we show comparisons of the meson mass using the perturbative shift A with 
the kinetic mass. Usually the agreement is good; the worst disagreements are still less 
than 2 combined a and occur for the quenched configurations at /5 = 2.575 (shown in the 
Figure) and for Nf = 2 at /3 = 2.1, Kgea = 0.1357. The dependence of the matrix elements 
on the heavy and light quark masses is analyzed in a partially quenched setup, i.e. setting 
the valence quark masses to their physical values (light, s and b) at fixed value of the 
sea quark mass, and then extrapolating to zero sea quark mass. In Fig. |321 on the left 
we show the light quark mass dependence of the zeroth order current at a typical set of 
parameter values. The light quark mass is parameterized as l/{2Kyai) — ^/i^Kcvai), where 
i^cvai denotes the point of zero valence quark mass at a given fixed sea quark mass. The 
data can be described well by a linear behavior. 

To correct the normalization of the light quark in the lattice operators for A^^\ we mul- 
tiplied Eq. ()219|1 with R = a/1 — Sk/ (4kc), the tadpole improved equivalent of Eq. ()2()()|1 

— (j) 

using uq = 1/{8kc). After interpolating the bare matrix elements fyM ,z = 1,2 to the 
strange or extrapolating to the light u, d quark mass, the renormalized matrix element 
is interpolated as a quadratic function in l/(aMo) to the b quark mass. An example is 
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Figure 40: Kinetic meson masses determined from ()18H1 as a function of momentum, 
from [174]. 

shown on the right in Fig. 

The dependence of the decay constant on the sea quark mass is mild. At f3 = 1.95, no 
real dependence can be resolved, whereas at /? = 2.1 there appears to be a slight increase 
if the sea quark mass is taken to zero. A linear fit is shown in Fig. 0^1 

In addition, fs and fn were calculated using clover heavy quarks with the Fermilab 
non-relativistic reinterpretation [181]. Current operators at zeroth and first order in the 
1/M expansion were used: 

= qixH^M^), (221) 

= g(a:)75747-^gh(x). (222) 

Local and smeared hadron operators were used at the source; operators at the sink were 
always local. Amplitudes and masses were extracted from two-exponential fits to axial 
vector-pseudoscalar and pseudoscalar-pseudoscalar correlators. Since the meson energy 
at finite momentum was affected by rather large statistical errors, the mass of the meson 
was calculated by adding the tree-level difference of the kinetic and pole heavy quark 
masses and to the meson pole mass, instead of using the kinetic mass of the 
heavy meson: 

aM2 = aMi + {aM^ - aM^), (223) 

where the masses with superscript Q are quark masses and those without superscript 
meson masses. 

Resumming the logarithmic mass dependence of the matrix elements (see Eq. ()135|l ) 
one can define the matrix element $ which is, according to HQET, in the heavy quark 
limit constant up to power corrections in 1/Mq. It is normalized such that it agrees with 
fsV^B if the heavy quark mass is the b mass. For a general heavy-light meson mass Mp 
it has the form: 

$(Mp) = {as{Mp)/as{MB)f''f^/M^, (224) 
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Figure 41: Dynamical meson masses determined using the dispersion relation II8II (Mi.4„) 
compared to the masses using one-loop values for the mass shift A (Mpgrt) at the smaller 
gauge coupling, the lightest sea quark mass, and a light valence quark mass around the 
strange. 



with bo = 11 — ^Nf. The heavy and light quark mass dependence of the resummed matrix 
element $(Mp) was fitted as a polynomial in the light quark mass ruq and the inverse 
heavy meson mass: 



i^/^<l>{Mp) = Ao + Aiarug + A2{amgf + — — [Bq + Biam, 

alVIf 



{aMp) 



;Co. (225) 



8.2.2 Results for fp and fp. 

The results for decay constants are listed in Table El The first error is statistical, the 
second a and perturbation theory errors, and the third error comes from setting the 
strange quark mass using K or (p input. The agreement of the results from NRQCD and 
Fermilab heavy quarks is an important result since one might expect that discretization 
and perturbative errors are different between the two formalisms, although naively one 
would expect them to be of a similar magnitude. 

In ratios of unquenched and quenched matrix elements calculated at the same value of 
the lattice spacing one may expect that most of the discretization and perturbative errors 
cancel. From the NRQCD calculation we find 

%po = 1.07(5); ^ = 1.10(5). (226) 

f B f Bs 

where only statistical errors are considered. If the scale is set with mp, we find a ~ 10% 
unquenching effect on the decay constants from the NRQCD calculation. 
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Figure 42: Dynamical b decay matrix elements from [174] at the smaller gauge coupling 

{/3 = 2.1) and a sea quark mass around the strange. On the left, fit of fv^M at a 
heavy quark mass close to the 6 as a function of the VWI quark mass (left) and the 
renormalized matrix element fy/M with the light valence quark mass interpolated to the 
strange (right). The b quark mass corresponds to l/(aMo) = 0.36. 

In the calculation with Fermilab heavy quarks we find 

1.11(6); %^ = 1.14(5); 

%=2 (227) 
1.03(6); %^ = 1.07(5), 

Again, the scale is set with rup. 

With the Fermilab heavy quarks, the effect of setting the scale with /^r has been 
studied. This leads to a better scaling of the decay constants as function of a [137]. If 
one takes the central value at the smallest lattice spacing for the result, the unquenched 
value remains the same as with the p scale, whereas the quenched result for fs with a 
from changes from 188 to 175 MeV. 

8.2.3 The flavor breaking ratio 

The extrapolation to the chiral limit could introduce a considerable uncertainty. E.g. 
Ref. [122] find a decrease of Jb by almost 20 MeV if they use a chiral extrapolation 
supported by instead of a simple polynomial fit. We use linear extrapolations in the 
partially quenched NRQCD analysis and a linear and quadratic fit for the simultaneous 
extrapolation in sea and light valence quark masses in the analysis with Fermilab heavy 
quarks. The result is 

fBjfB = 1.179(18)(23), Nj = 2, 



Nf=2 
Nf^O 
Nf=2 
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Figure 43: Sea quark mass dependence of at /3 = 2.1, from [174]. 

fsJfB = 1.150(09)(20), iV^ = (228) 

from NRQCD, and 

fsJfB = 1.203{29){28){tl')Nf = 2 
fsJfB = l.U8mmi-f)Nf = 
fojfn = 1.182(39)(25)OiVy = 2 

toJlD = 1.138(05)(18)(g9)iV^ = (229) 

from heavy clover quarks. The first error is statistical, the second systematical (using 
only discretization effects of the clover quarks since O(a^) and 1/M errors should largely 
cancel in the ratios), and the third from the difference in setting the strange quark mass 
using the K or (j) meson. 

The flavor breaking ratio fss/fB has with both heavy quark formalisms no significant 
unquenching effects, setting the scale with nip. 

8.3 Discussion 

For the decay constants and some of the mass results presented here, matching has been 
performed at 0{as). In calculations of 1/M corrections in HQET and of higher twists 
in operator product expansions relatively large higher-loop corrections may be expected. 
On the lattice, higher dimensional operators generate power divergences under renormal- 
ization which are difficult to cancel in perturbation theory [284]. One might worry how 
precisely 1/M corrections can be quantified using effective theories on the lattice. The 
size of the higher order corrections in perturbation theory depends on the choice of renor- 
malization scheme and scale used at one loop, ay at the scale q*, described in Sec. 13. 6[ 
should be a rather optimized choice. 
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The 1/a terms affect the simulation energies as well as the matrix elements. For 

example we fit the matrix elements f\^M^^^ from Ref. [174]^, with the light quark mass 
interpolated to the strange with equal sea and valence quark masses, according to 

;(:.Vm; = a + A + ^. (230) 

where Mq is the bare heavy quark mass, which yields the fit parameters given in Table IT71 
The terms are unphysical since the action is only 1/M improved, they are however 

also on finer lattices not very large. The lattice spacing dependence gives an idea of 
to what extent the 1/a terms contribute to the matrix element. Whether the power 
contributions make out a substantial part of the decay constant can only be seen if the 
renormalization constants are included. 

In the calculation [101] using 1/M^ NRQCD and clover quark actions and Wilson 
gluons at a ~ 0.1 fm, we find for the ratio of the unrenormalized Bs matrix elements 
/^^^)yM^//5?y^ ~ 0-12 and for ~ 0-02 or less for z = 3,4,5. 

With renormalization group improved gluons [174], using 1/M NRQCD and clover quark 
actions, the ratio of matrix elements has a similar value, f^^ Mb^ / f^] \/ Mb, ~ 0.1. 
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Figure 44: Ratio of decay constants with rotated and unrotated currents: 
yrotated^ yunrotated _ X as a fuuctiou of the iuversc heavy meson mass. 

To estimate the accuracy of the perturbative matching calculation it is of interest to 
how well the cancellation of the as/ (aM) contribution works at one loop. As discussed 
for NRQCD in more detail in [179], for ~ 0.2 — 0.25 and the perturbative coefficient to 

^The unquenched matrix elements in Tables XIV- XIX of Ref. [174] arc mistakenly arranged in the 
order opposite from the valence quark masses 
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the radiative correction of A\ proportional to l/(aMo) around —(0.3 — 0.4), the physical 
^qcd/M contribution of the current correction is less than 50% of the bare current 
correction. However, there is no indication of very large higher order corrections. 

The change from the bare to the one-loop decay constant is around 6 — 10% in [79] 
using an 0(1/M) action and 9 — 14% in [101] with an 0{1/M^) action and Wilson gluons 
at = 2 GeV, however less than 5% with a 1/M NRQCD action and renormalization 
group-improved glue [174]. 

In contrast for relativistic heavy quarks [137], the 1/M current operator Eq. ()222|1 
contributes for B as well as for D mesons around 4 — 6% to the decay matrix element, 
see Fig. EH 
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Figure 45: Fit results for the heavy-strange matrix element $ (Eq. ()225p ). with the 
strange quark mass determined using the K meson from [137], as a function of the inverse 
heavy meson mass. Lines with different symbols describe fit functions from different 
lattices. The symbols aare just to giude the eye and do not correspond to actual data 
points. The two vertical lines indicate the Bg and Dg meson masses. 

To see the extrapolations to the static values and to be able to estimate the 1/M 
corrections at a range of values of the heavy quark mass, we show in Fig. |^ the fits of 
the quenched FNAL results as a function of the heavy meson mass for all the quenched 
lattices. This will be compared to other quenched lattice results later. 

Concerning the theoretical estimates of the discretization and perturbative errors used 
in the NRQCD calculations one might object that it is not obvious which numerical value 
of Aqcd should be used for the error estimate (e.g. quenched perturbative values from 
the determination of Og are around 200 — 300 MeV, while generic momentum scales of 
QCD are given by hadron mass splittings, typically ~ 500 MeV). Therefore we compare 
the for Aqcd = 600 MeV (used for the error estimate in the paper) and 400 MeV. The 
unquenched results in both formalisms are shown in Fig. H^l Finally we summarize our 




102 



450 

400 
350 

> 

S 300 - 
250 - 
200 ■ 
150[ 



■ NRQCD 
□ FNAL 



0.2 



0.4 



0.6 
a [GeV"' 



0.8 



450 
400 

350 - 



> 

at 

2 300 

^" 

250 
200 
150 



■ NRQCD 
□ FNAL 



0.2 0.4 0.6 

a [GeV"' 



0.8 



Figure 46: Jb as a function of a. Solid error bars: statistical errors, dashed error bars: 
systematical errors using Aqcd = 400 MeV (left) and Aqcd = 600 MeV (right). 



results on Jb and fs^ in Fig. |37| A very important finding of the studies [137,174] is that 
the B and Bg decay constants using NRQCD and Fermilab heavy quarks agree on the 
same set of lattices on a range of lattice spacings. 
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Figure 47: Summary of own results on fs (left) and fs^ (right), from Refs. [79, 124] 
(SCO) [101] (GLOK), [137] (CP-R) and [174] (CP-NR). Diamonds denote quenched, 
squares Nj = 2 results. 



8.4 Comparison to other work 

We now compare our results to the ones from other lattice simulations. Lattice values 
for Jb and Jb^ since 1998 are listed in Table E| The first error shown in the Table is 
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statistical, the second is the systematical uncertainty, adding all the contributions quoted 
by the authors in quadrature. 

In the quenched approximation, a wealth of lattice results on Jb is available, and we 
compare our results from NRQCD with two calculations using non-perturbative renormal- 
ization using the Wilson gauge action, from the ALPHA [210-212] and the UKQCD col- 
laboration [89]. The work of ALPHA [210-212] uses the Schrodinger functional approach 
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Figure 48: Resummed matrix elements from non-perturbatively 0(a) improved clover 
quarks. For the details of the resummation see the text. The result at the Bs mass 
is obtained by interpolation between charm and static heavy quarks for the ALPHA 
data [211] or extrapolation from the charm for the UKQCD data [89]. The data in the 
Figure use a set with tq. 

to simulate non-perturbatively improved clover charm and static quarks and obtains fs^ 
from an interpolation between the values at infinite and at the charm mass. The results 
are continuum extrapolated. The scale is set using r^. 

Ref. [89] uses two lattice spacings. The result from the finer lattice is their contin- 
uum estimate. An interpolation to a static value is not performed, since non-perturbative 
renormalization was not available for the static action at the time of their calculation. In 
Fig. I3H1 we compare the slopes of both nonperturbative calculations. Both use a mass- 
independent renormalization scheme. The matrix elements themselves agree around the 
charm mass, their slope is slightly different. Ref. [89] uses a one- loop resummation accord- 
ing to Eq. ()224|1 . and [211] uses the three-loop anomalous dimension of the static current 
in terms of the renormalization group invariant quark masses. To calculate the anomalous 
dimensions for the plots, the value of Aqcd = 238 MeV is used. Using the three-loop 
anomalous dimension results in a slightly steeper running than using the 1-loop formula. 
Adding results from finer lattices and larger masses to the calculation of [89] would be 
useful to clarify the situation. 
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In Fig. |3ni we compare quenched heavy-strange decay matrix elements from 
NRQCD [101] with the non-perturbatively renormahzed relativistic results. The results 
have been rescaled to take the quenched scale ambiguity into account. Since the lattice 
spacing of the NRQCD calculations is determined using the p mass, the ALPHA value in 
the Figure is for comparison decreased by 12%, which is the change they quote for fs^ if 
a is changed by 10% [211], and their static result by 15%. A static value obtained from 
extrapolation of the 1/M NRQCD results from [79] to infinite mass aX j3 = 6.0, and the 
non-perturbatively improved static value from ALPHA [210] aX j3 = 6.0. The difference 
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Figure 49: Comparison of NRQCD and perturbatively improved static decay matrix 
elements from [79,101] with non-perturbatively improved results from [89,211]. The non- 
perturbatively improved static value from j3 = 6.0 of [210] is denoted by the star. The two 
data sets from NRQCD at 0(1/M^) correspond to 1-loop renormalization with q* = 1/a 
and q* = 7i/a. 

of the NRQCD resuhs at 0(1/M) and 0(1/^^) is slightly less than 2 combined a, which 
we do not consider as significant since one of the ensembles is very small. Ref. [178] com- 
pares NRQCD results 0(1/ M) and 0(1/M^) on identical gauge field ensembles, and finds 
very good agreement. The effects seem to be small. With a quadratic fit including 

results in a region around the charm mass, [212] finds fs^ = 198(8) MeV which is close 
to the published value of 205 MeV from [210]. 

NRQCD results are obtained at a finite lattice spacing, and using different gauge 
actions can lead to different discretization and perturbative errors. In Fig. we see that 
our results from renormalization group improved gluons and FNAL heavy quarks [137] 
extrapolate to a relatively high quenched value also on finer lattices. As mentioned before, 
the NRQCD results agree at the b mass with the FNAL results. 

The extrapolation to the static limit is higher than the nonperturbatively renormahzed 
static result from the ALPHA collaboration using Wilson glue which can be taken as an 
indication of remnant discretization and perturbative errors. If the scale is set using f^^ 
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instead of nip, the value decreases by 7%. 

Results and global averages are listed in Table ^1 For the averaging, we employed 
the method used in the Particle Data Book for averaging the experimental results (see 
Ref. [146], Introduction). The statistical and systematical errors were added in quadra- 
ture. Taking into account that sets of results have the same systematic error, the weights 
were computed and the averages obtained. The estimation of the systematic errors of all 
lattice results is described in [249]. To briefly introduce the other calculations, Ref. [178] 
uses NRQCD at 0{1/M) and 0(1/M^) with one-loop tadpole-improved clover light quarks 
and quenched Wilson glue, and [202] uses quenched NRQCD with an 0{a?) tadpole im- 
proved gauge action and staggered light valence quarks. Ref. [179] uses quenched NRQCD 
at 0(1/M^) and tadpole-improved and non-perturbatively 0(a) improved light clover 
quarks. The quenched calculations of Refs. [103,257,258] use Wilson gauge and clover 
quarks at a = 0.35 — 0.37 GeV~^ {(3 = 6.2) simulated in the charm mass region, simi- 
lar to [89]. Ref. [257] uses a non-perturbatively 0{a) improved clover quark action and 
a partly non-perturbative current renormalization. Ref. [258] uses tree-level tadpole- 
improved clover quarks without including 0{as x a) current correction terms. Ref. [103] 
uses nonperturbative 0{a) improvement except for a perturbative value for the 0{asamq) 
quark mass correction to the renormalization constant. Refs. [259-261] use heavy quarks 
in the Fermilab formalism and extrapolate their results to a — 0. Refs. [262,263] uses 
a step scaling function with the Schrodinger functional scheme and clover heavy quarks. 
A part of the renormalization factors is calculated nonperturbatively. Their results are 
continuum extrapolated. Among the two-flavors results, Refs. [122,267] use NRQCD and 
Fermilab heavy quarks respectively, at a fixed value of the gauge coupling. 

Refs. [259-261] use heavy quarks in the Fermilab formalism and extrapolate their 
results to a ^ 0. The quenched calculations of Refs. [89,103,257,258] use Wilson gauge 
and clover quarks at a = 0.35 — 0.37 GeV~^ (/? = 6.2) simulated at the charm quark 
mass. Ref. [258] uses tree- level tadpole-improved clover quarks without including 0{asa) 
terms in the renormalization. Ref. [257] uses a non-perturbatively 0(a) improved clover 
quark action and a partly non-perturbative current renormalization. Refs. [89, 103] use 
nonperturbative 0(a) improvement except for a perturbative value for the 0{asamq) 
quark mass correction to the renormalization constant. Although different degrees of 
improvement are used, and the scaling behaviour is found to be different, the results for 
agree at /? = 6.2. A calculation from CP-PACS using the heavy quark formalism 
of [147] gives a preliminary result of around 210 MeV for Jb^ [268]. 

The unquenching effects on fs depend on the details of the calculation, e.g. how 
the scale is set. In Table ^1 ratios of decay constants from quenched and two-flavor 
simulations with the same gauge field and valence quark actions are listed. One finds an 
increase of ~ 10% with unquenching if a is set with 10 — 20% if a is set with nip. On 
the other hand, Ref. [269] quotes no significant difference if the scale is set with tq for 
fjj^, and Ref. [174] finds no significant difference between quenched and two-flavor results 
if the scale is set with the P — S splitting of the T, however on coarse lattices. A new 
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result with two flavors of domain- wall quarks on the ratio fss/fB of [270] is somewhat 
higher than the CP-PACS results of [137,174], but not statistically significant. 

To investigate the effect of an additional dynamical strange quark, we compare our 
result to the published NRQCD lattice results with 2 + 1 flavors of staggered sea quarks 
[264,265]: 

fB, = 260(7) (28) MeV, 

/b = 216(9) (20) MeV, (231) 

where the flrst error is statistical and the second is from the systematic uncertainties 
added in quadrature. Calculating at small quark masses, Ref. [265] succeeds in achieving 
a precise result for the flavor breaking ratio of 1.20(3). Calculations of the decay constants 
with the staggered Nf = 2 + 1 conflgurations using FNAL [271] heavy quarks are in further 
progress. 

So far there is no indication for large flnite volume effects in the decay constants. 
Chiral perturbation theory predicts a flnite size effect of < 1% for lattice sizes ~ 1.6 
fm and pion masses around 500 MeV for quenched QCD and much less than 1% for 
Nf = 2 [167]. Ref. [261] estimates the flnite volume uncertainty to be of the order of 
several percent; the estimate is however from a comparison with very small volumes. 

We also use the experimental value for f^,^ from unquenched lattice to fs^ and the 
lattice ratio JbJ Jds from two-flavor calculations which work directly at the b and c quark 
masses without using extrapolations to calculate Jb^- Taking the experimental value 
= 283(45) MeV using Eq. (p?T|l and /bs/ fo^ = 0.93(|) from the unquenched results 
from Table IT^ one flnds Jb^ = 263(gf) MeV, which is in agreement with the lattice result 
of [264]. 

Other recent review articles [272-276] quote lattice estimates for Jb and fs^ which are 
within errors in agreement with the averages quoted in Table IT9l 

Within errors the lattice results are also in agreement with recent sum rule [277-280] 
and potential model calculations [281], see Tab. QUI 

9 SUMMARY AND CONCLUSIONS 

This work presents a study of decay constants and the hadron spectrum in lattice QCD. 
The physics of hadrons with b and light {u, d, s) quarks is of great relevance for the deter- 
mination of parameters of the Standard Model, but both quarks are difficult to simulate 
because the light quark physics is sensitive to flnite volume effects and the heavy quark 
physics is sensitive to lattice spacing effects. Heavy quarks can be simulated effectively 
using NRQCD and HQET. The remaining difficulty is then to relate the light quarks 
simulated on the lattice to the physical light u and d quark, since extrapolations over 
a more or less large range of masses are required. Chiral perturbation theory provides 
a useful description of hadron masses directly as a function of the masses of the pseu- 
doscalar mesons built from the quarks under consideration. The simulation parameters 
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corresponding to the physical quark masses can be determined using the hght hadron 
spectrum. Of course the hght hadron spectrum itself, in particular the nucleon mass, is 
one of the most basic quantities of which the lattice has to be able to show that they can 
be explained by QCD. 

Light hadrons were studied in the quenched and two-flavor approximations using im- 
proved gluon and light quark actions. There are no indications of disagreement of the 
lattice results with experiment within the estimated systematic errors. Our results on 
the low-lying hadron spectrum with experiment agree surprisingly well with experiment 
already in the quenched approximation. Now further clarification of the question whether 
there is really a discrepancy between various improved fermion actions, would be of in- 
terest. Further one should note that quark model interpolating operators for the hadrons 
work very well. 

Chiral extrapolations and finite volume effects of the nucleon mass find agreement 
with O(p^) chiral perturbation theory predictions at rather large quark masses, while 
chiral perturbation theory is assumed to be valid for pion masses up to ~ 600 MeV or 
less. In particular in theories with dynamical quarks, lattice results at smaller quark 
masses would be of great interest to be able to make solid chiral extrapolations of hadron 
masses and matrix elements. 

The agreement of the spectrum with experiment is an important result corroborating 
the reliability of the lattice calculations of the corresponding weak matrix elements. Most 
of the lattice results on the heavy-light spectrum have been obtained using non-relativistic 
effective theories, in particular NRQCD. The agreement with experiment is a considerable 
success of the physical pictures used in these effective theories. For example it implies 
that the strange-light quark mass difference determined from light physics is the same as 
in heavy-light physics. 

For decay constants, comparison with results using different lattice discretizations and 
0, 2 and 2+1 flavors of light sea quarks were made. Effective theories, either those following 
the ideas of NRQCD and the Fermilab method, or HQET together with interpolations to 
smaller masses are presently still indispensable in lattice calculations of b hadrons. Using 
both methods next to each other it should be possible to obtain quite reliable estimates 
of hadronic matrix elements. 

For heavy-strange decay constants an error not much larger than 10% appears real- 
istic. While the magnitude of the error has not been reduced much compared to lattice 
errors quoted around 10 years ago, the error has become significantly more reliable. In 
unquenched calculations with heavy hadrons, effective theories will continue to play an 
important role also in the coming years. While using NRQCD a continuum limit is not 
possible and renormalization constants are only calculated perturbatively, the scaling 
studies of previous calculations have shown that the error estimates are reliable. It would 
be important to support the simulations at the b quark using non-relativistic effective 
theories also for dynamical quarks with simulations of static and charm quarks using 
non-perturbative renormalization, e.g. using the Schrodinger functional approach. Since 
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continuum extrapolations also involve uncertainties and since dynamical simulations on 
fine lattices are presently still expensive, nonrelativistic methods such as NRQCD and 
the Fermilab method are still important in precision calculations in B physics. 

Despite many lattice studies of a very broad range of observables there is no indication 
that QCD is not a consistent theory of the strong interactions. Comparison of of chiral 
perturbation theory in comparison with lattice data and experiment supports the present 
picture of chiral symmetry breaking in QCD. Increase of precision on a number of masses 
and matrix elements calculated on the lattice will be necessary to resolve the signatures 
of new physics. Important directions to investigate are unquenching and the chiral limit. 
Further, the lattice has achieved a significant potential for calculations of many quantities 
in hadron physics such as masses of excited hadrons, exotic states, a variety of form factors 
and structure functions, and the investigation of confinement and chiral symmetry. 
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Ref. 


AAf(fi)[MeV] 


AAf(5,)[MeV] 


B*-B 


Lattice 


[88], a ~ l.lGeV-i 

[88], a ~ 2.6GeV-i 

[80] 

[203] 

202 11 

[202], 11 

[202], 11 

[202], 11 

[202], Asq 

[202], Asq 

[202], Asq 

[235] 

[234] 

average 


374(37) 
475(1^) 

408(67) 
402(41) (f) 


400(30)(J9) 
295(45)(Ji) 
357(27) (0) 

451(li) 

442 f 56) 

471(25) 

315(105) 

425(60) 

285(78) 

523(94) 

403(56) 

386(31) 

419(37) 

382(23)(?D 


Model calculations 


[242] 


450 


440 


[191] 
[248] 


453 
427 


466 
431 


b;-b 


Lattice 


[88], a ~ l.lGeV-i 

[88], a ~ 2.6GeV-i 

[80] 

[203] 

[198] 

[235] 

[234] 

average 


526(45) 

493(1) 

426(17) 

440(77) 
498(49) (f) 


402(78)(g0) 
474(62) (J6) 
493(26) 
478(i) 

534(52) 
455(41) 
487(31)(r) 


Model calculations 


[242] 


450 


490 


[191] 
[248] 


453 
435 


469 

447 



Table 11: B2 — B splittings. Only statistical errors and, where quoted by the authors, 
errors due to chiral extrapolation and fixing the b quark mass are shown. 
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Ref. 


AM(5)[MeV] 


AM(S,)[MeV] 


Lattice 


[88], a ~ l.lGeV^i 

[88], a ~ 2.6GeV"i 

[80] 

[203] 

[235] 

[234] 

average 


155(32) 
18(i) 

32(87) 
98(47)(J0) 


41(94)(J4) 

179(65)(^) 

136(23) 

27(i) 

148(61) 

36(55) 

101(25)(ii) 


Model calculations 


[242] 
[243] 
[191] 
[248] 


40 

-155 

-5 

7 


50 
3 

16 



Table 12: P state fine structure of B mesons. Errors include statistical uncertainties and 
uncertainties due to chiral extrapolation and in determination of the h quark mass. 
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Nf = 


lib 


Ref. 




Ab-B 


^b-Ab 


n - 


SGO 

GLOK 

CP-PACS 


[197] 
[80] 
[80] 


410(10) 
370(67(1^) 
361(22) 


221(71)(1^) 


19(7)(^) 


Isb 


Ref. 










GLOK 


[80] 


392(50) (?5)(i^) 


186(51)(l}?)\l?o) 


19(5) (1) 


ssb 


Ref. 








— ^b 


GLOK 


[80] 




18(4)(1) 




Nf = 2 


lib 


Ref. 




Ab-B 


Sfe-A, 




SGO 

CP-PACS 


[124] 
[80] 


545(40) (22) 
417(19) 







Table 13: Own results on b baryon masses and mass splittings for various light quark 
flavor combinations. The first error is statistical. Where applicable, the second error 
comes from the variation of a"^ between 1.8 and 2 GeV, and the third from fixing the 
strange quark mass using the K* or $ instead of K. 
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Ref. Ab-B[MeV] S;, - A^MeV] 


Lattice 


[233] 
[80] 

[201], as = l.lGeV-i 
[201], as = 0.9GeV-i 
[236] 

[198], quenched 
[124],iV/ = 2 
[198], Nf = 2 
quenched average 


338(1) 
370(67) 

389(44) 
361(22) 
545(40) (22) 
417(19) 
372(33) (f) 


221(71) 
156(1) 
191(3?) 
122(65) 

174(28)(i5) 


Models 


[244] 
[145] 

[191,246] 308 
[145] c quark 


312 


196 
217 
202 
212 


Experiment 


[146] 


311(10) 




Ref. S*-S5[MeV] - fi5[MeV] 


Lattice 


[80] 

[201], as = l.lGeV-i 
[201], as = 0.9GeV-i 
[236] 
average 


19(7) 
22(12) 

mi) 

10(12) 
18(8)(^) 


18(4) 

18(^) 
20(9) 
7(4) 

13(7)(^) 


Models 


[145] 
[244] 
[246] 


10 

8 

29 


23 



Table 14: b baryons. Only statistical errors and, where applicable, systematical errors 
due to scale setting except for the quenched scale ambiguity and fitting are shown. 
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aMo 






aMps{Ks) 


Mps{Ks)[GeV] 


aq* = 1 


aq* = TV 


1.71 


1.73(10) 


0.513(3) 


2.24(10) 


4.46(22) 


0.169(4) 


0.172(4) 


2.0 


2.02(9) 


0.515(4) 


2.54(9) 


5.05(21) 


0.169(3) 


0.174(3) 


4.0 


4.07(5) 


0.528(4) 


5.00(5) 


9.95(22) 


0.169(3) 


0.187(4) 


8.0 


7.63(16) 


0.527(6) 


8.16(16) 


16.24(45) 


0.167(5) 


0.195(5) 



Table 15: Parameters and results for strange heavy-light mesons from [79]. The last two 
columns show renormalized decay matrix elements a^^'^f\/M for matching scales q*. 





Nf = 2 


NRQCD 


fniMeV] 
/sjMeV] 


191(4)(27) 
220(4)(31) 


204(8)(29) 
242(9) (34) 


relativistic 


/B[MeV] 
/fijMeV] 
/z?[MeV] 
/dJMcV] 


188(3)(26) 
220(2)(31)(1«) 
218(2)(39) 
250(1) (45) (+[;^ 


208(10)(29) 
250(10)(35)(i^) 
225(14)(40) 
267(13) (48) (1°) 



Table 16: Heavy-light decay constants from [137,174]. The scale is set using nip. 



Matrix el. 




Ai [GeV3/2] 


A2 [GeV5/2] 






1.95 


1.236(3) 


-1.977(26) 


2.03(5) 


f^^VMp 


2.1 


0.912(30) 


-0.72(31) 


-0.08(79) 


1.95 


-0.0134(21) 


-0.402(17) 


0.329(35) 




2.1 


-0.0086(19) 


-0.422(19) 


0.395(50) 



Table 17: Fit parameters of the unrenormalized matrix elements as a function of I/Mq 
according to Eq. ()230|) . 
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Ratios of decay constants 

Ref. 




[261] (iV/^ 
[210] {Nf 

[137] (iV^: 

[261] (iV;: 



0) 
0) 

2) 
2) 



0.88(1) 

0.891(12)(^0) 

0.81(6) 

0.94(6) 

0.922(13)(i) 



Table 18: Ratios of decay constants. The first error is statistical, the second the system- 
atical errors given by the authors added in quadrature, where applicable. 
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Ref. 



scale 



/B[MeV] /BjMeV] 



Lattice 



Nf = 



[259] 




nip 


173(4)(13) 


199(3)(14) 


[260] 




Ik 


164(lf)(8) 


185(^3) (9) 


[101] 




nip 


147(11)(13) 


175(8)(16) 


[178] 




Va = 427MeV 


170(5)(15) 


191(4)(17) 


[257] 




M* , , 


173(13)(i^) 


196(11)(^2) 


[89] 




u 


195(6)(i) 


220(6)(i) 


[258] 




fx 


177(17)(22) 


204(12) (i) 


[179] 




nip 




187(4)(15) 


[137] 




nip 


188(3)(26) 


220(2)(i) 


[174] 




nip 






[103] 






174(22)(8) 


204(15)(8) 


[261] 




fn 


173(6)(16) 


199(5)(i) 


[262] 




To 


170(11)(23) 


192(9)(25) 


[263] 




ro 




192(6)(4) 


[210] 




ro 




205(12) 


[202] 




ro 




225(9) (34) 


averaj 


^el 




175(7)(f[''^) 


198(5)(tJ^]) 


averaj 


ie2 






201 (6) (ig'^'^') 








= 2 




[124] 




nip 


186(5)(25) 


215(3)(i) 


[137] 




nip 


208(10) (29) 


250(10)(i) 


[174] 




nip 


204(8)(29) 


242(9) (34) 


[261] 






mmfr) 


217(6)(i) 


[122] 




nip 


mml'2) 


215(9)(1^) 


[267] 




ro = 0.49fm 


181(7)(i^) 




averaj 


lel 




190(10)(i5''') 


226(15)(ft^l) 


averaf 


ie2 






226(15)(ft^l) 



2 + 1 



[264,265] T(25-l^),ri 216(9)(20) 260(7)(28 



Table 19: fs and from the lattice. Statistical errors and systematical errors given by 
the authors are included, adding the systematical errors in quadrature. The method to 
set the scale is indicated in the second column. The first error on the averages is due to 
the statistical and systematical errors of the individual results, while the second error is 
from chiral extrapolation uncertainties and scale ambiguities as explained in the text. 
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Ref. 


/^[MeV] 


/B.[MeV] 


Sum rules 


[277] 


203(23) 


236(30) 


[278] 


206(20) 




[279] 


210(19) 


244(21) 


[280] 


180 - 190(30) 




Potential models 


[281] 


178(15) 


196(20) 



Table 20: Jb and Jb^ from sum rule and potential model calculations. 
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